• 한국게임학회 논문지
• /
• 제12권5호
• /
• pp.67-78
• /
• 2012

합산 영역 테이블은 모든 픽셀에 대해 임의의 크기 사각영역의 이미지 필터링 처리를 일정시간 안에 가능케 한다. 이러한 특성은 각각의 픽셀에 대해서 주변 픽셀의 밝기의 합 혹은 평균을 필요로 하는 이미지 처리 적용 분야에 유용하게 쓰일 수 있다. 합산 영역 테이블의 생성은 단지 행 혹은 열 단위의 합만을 구하는 메모리 바운드 작업임에도 불구하고 기존 연구들은 이미 존재하는 데이터 병렬성만을 활용하기 위하여 대기 시간이 긴 전역 메모리에 과도한 접근을 하여야만 했다. 본 논문에서는 입력 데이터를 정방의 서브 이미지로 분할하고 매개 데이터를 이들 간에 파급시킴으로써 GPGPU 환경 적합한 알고리즘을 제안하고자 한다. 이를 통하여 기존 방법 대비 전역 메모리 접근 량을 거의 반으로 줄임으로써 주어진 메모리 대역폭을 효율적으로 사용한다. 결과에서도 성능이 대폭 향상되었다.

• 대한기계학회논문집A
• /
• 제26권3호
• /
• pp.576-585
• /
• 2002

A three-dimensional Lagrangian explicit time-integration finite element code for analyzing the dynamic impact phenomena was developed. It uses four node tetrahedral elements. In order to consider the effects of strain rate hardening, strain hardening and thermal softening, which are frequently observed in high-velocity deformation phenomena, Johnson-Cook model is used as constitutive model. For more accurate and robust contact force computation, the defense node contact algorithm was adopted and implemented. In order to evaluate the performance of the newly developed three-dimensional hydrocode NET3D, numerical simulations of the oblique impact of mild steel plate by mild steel sphere were carried out. Ballistic limit about various oblique angle between 0 degree and 80 degree was estimated through a series of simulations with different initial velocities of sphere. Element eroding by equivalent plastic strain was applied to mild steel spheres and targets. Ballistic limits and fracture characteristics obtained from simulation were compared with experimental results conducted by Finnegan et al. From numerical studies, the following conclusions were reached. (1) Simulations could successfully reproduce the key features observed in experiment such as tensile failure termed "disking"at normal impacts and outwards bending of partially formed plus segments termed "hinge-mode"at oblique impacts. (2) Simulation results fur 60 degrees oblique impact at 0.70 km/s and 0.91 km/s were compared with experimental results and Eulerian hydrocode CTH simulation results. The Lagrangian code NET3D is superior to Eulerian code CTH in the computational accuracy. Agreement with the experimentally obtained final deformed cross-sections of the projectile is excellent. (3) Agreement with the experimental ballistic limit data, particularly at the high-obliquity impacts, is reasonably good. (4) The simulation result is not very sensitive to eroding condition but slightly influenced by friction coefficient.

• 한국측량학회지
• /
• 제27권1호
• /
• pp.693-700
• /
• 2009

장기선 네트워크 RTK(Real-Time Kinematic) 측량, 정밀단독측위(precise point positioning) 및 전리층/대류권 지연 모니터링등 GPS를 이용한 실시간 데이터 처리를 위해서는 IGS(International GNSS Service)에서 제공하는 정밀궤도 수준의 정확도가 시간지연 없이 확보되어야 한다. 본 연구는 준 실시간 위성궤도결정을 위한 선행연구로서 일반적인 위성궤도 결정 방법, 특히 동역학적 방법에 대한 이론적 고찰과 가속도 적분을 위한 지구기준/관성좌표계 변환 방법에 대한 테스트를 수행하였다. IAU 1976/1980 세차/장동모델은 IAU 2000A 모델과 0.05mas 수준의 허용범위 내 차이를 보였다. IAU 2000A 모델은 기본적으로 복잡한 장동모델로 인해 변환행렬 계산에 많은 시간이 소요된다. IAU 2000A 모델에 기존 변환 방법을 사용하는 경우가 NRO(non-rotating origin) 방법에 비해 2배정도 빠른 결과를 보인 반면 회전행렬의 실질적인 차이는 없는 것으로 나타났다.

• 한국컴퓨터그래픽스학회논문지
• /
• 제5권1호
• /
• pp.27-33
• /
• 1999

본 논문은 유연한 비정형 물체의 애니메이션을 위한 효율적인 기법을 제안한다. 비정형 물체를 표현하기 위해 질량-스프링 모텔이 사용되었다. 지금까지 많은 기법들이 부드러운 객체의 사실적인 애니메이션을 생성하기 위해 질량-스프링 모델을 사용하였다. 질량-스프링 모텔의 애니메이션을 수행하기 위한 가장 손쉬운 접근법은 명시적 오일러 방법 (explicit Euler method)인데, 이 방법은 '불안정성 문제'라는 잘 알려진 문제가 발생한다는 단점을 가지고 있다. 이 불안정성 문제를 해결하기 위한 해법으로 암시적 적분법이 사용될 수 있다. 그러나, 이 암시적 방법의 가장 결정적인 약점은 대규모의 선행 시스템을 풀어야 한다는 것이다. 본 논문은 암시적 방법의 근사(approximation)를 이용하여 질량-스프링 모델을 빠른 시간에 애니메이션 할 수 있는 기법을 제시한다. 제안된 기법은 n 개의 질점이 O(n) 개의 스프링으로 연결되어 있을 때, 각 질점의 상태를 O(n) 시간에 안정적으로 갱신할 수 있다. 본 논문은 사실적인 결과를 위해 비정형 물체와 공기와의 상호 작용도 고려하였다.

• 한국항공우주학회지
• /
• 제36권1호
• /
• pp.72-78
• /
• 2008

본 연구에서는 블레이드 구조 변형 효과를 고려하여 스테이터-로터 상호간섭 케스케이드 모델의 성능평가를 위한 유체-구조 연계해석 시스템을 개발하였다. 고정된 스테이터와 회전하는 로터는 상호간섭 영향이 유동해석에 고려되었으며, 레이놀즈-평균화 난류 방정식인 Spalart-Allmaras 모델과 k-ω SST 난류 모델이 압축성 유동박리 효과를 고려한 유동하중을 예측하기 위해 적용되었다. 정적인 유체-구조 연계해석과 수렴율 증진을 효과적으로 수행하기 위하여 큰 인공 감쇠를 가지는 연계 Newmark 시간적분 기법을 적용하였다. 수치실험을 통해 탄성축 위치에 따른 구조변형 효과가 케스케이드 성능에 미치는 영향을 파악하였다. 구조변형 효과가 고려된 경우 일반적인 강체 블레이드 모델에 대한 성능예측 결과와 다소 차이가 유발될 수 있음을 보였으며 공력탄성학적 영향을 고찰하였다.

• 한국방재학회:학술대회논문집
• /
• 한국방재학회 2011년도 정기 학술발표대회
• /
• pp.59-59
• /
• 2011

파랑의 전파와 변형에 대한 연구에는 수심방향으로 적분한 2차원방정식인 완경사방정식과 Boussinesq 방정식을 기반으로 한 수치모형을 이용한 연구가 최근까지 가장 활발하게 진행되어 오고 있다. 그러나 실제 구조물의 설계에는 2차원 수치모형에서 고려할 수 없는 수심방향 유속에 기인한 정확도의 문제로 인해 구조물의 형상과 재원을 설계하기 위한 정교한 수치모형실험이 어려워 주로 수리모형실험에 의존해 왔다. 수리모형실험은 실제 현상을 가장 잘 재현해낼 수 있어 신뢰성이 매우 높지만 다양한 실험을 수행하기가 어렵고 많은 시간과 비용이 소요되는 단점이 있다. 이에 따라 최근 수심방향으로 완전한 운동방정식인 Navier-Stokes 방정식을 푸는 3차원 수치모형에 대한 연구가 활발히 진행되고 있다. 이론적으로 매우 우수한 모형이긴 하나 정확도 높은 결과를 얻기 위해서는 매우 조밀한 격자를 필요로 하기 때문에 아직까지 막대한 계산시간이 필요하다는 단점이 있으나 컴퓨터 기술이 급격한 속도로 발전하고 있어 Navier-Stokes 방정식 모형의 적용 가능성은 계속 높아지고 있다. 파랑변형을 다루는 수치모형실험을 수행할 때 외부조파를 사용할 경우 구조물이나 지형에 의해 반사되어 나온 파랑이 조파지점에 도달할 때 실험영역으로 재 반사되는 문제가 발생한다. 이를 해결하기 위해 내부조파기법의 개발에 대한 연구가 필수적이었으며, 자유수면변위를 변수로 사용하는 모형의 경우 그 연구가 매우 활발하게 진행되어 왔다. 한편, Navier-Stokes 방정식 모형의 경우 자유수면변위를 변수로 사용하는 2차원 모형에 비해 상대적으로 연구가 미흡하였다. 본 연구에서는 기존의 연직 2차원 Navier-Stokes 방정식 모형에 사용된 연속방정식에 질량 원천항을 추가하는 내부조파기법을 도입하여 3차원 수치모형에서 고립파를 내부조파하고, 급경사에서의 고립파의 처오름 및 처내림 현상을 수리모형 실험결과와 비교 및 분석하였다. 수치모형은 Navier-Stokes 방정식을 엇갈림 격자체계에서 계산하는 동수압 모형으로서, Two-step projection 기법을 사용하는 유한차분모형을 사용하였다. 본 수치모형은 난류의 해석을 위해서 상대적으로 큰 에디(eddy)만을 고려하는 SANS(spatially averaged Navier-Stokes) 방정식을 계산하는 LES(large-eddy-simulation) 기반의 수치모형으로, 난류 모델링을 위해 Smagorinsky LES 모형을 사용한다. 또한, 압력장의 계산을 위해 Bi-CGSTAB 기법을 이용하여 Poisson 방정식의 해를 구하였으며, 자유수면 추적을 위하여 2차 정확도의 VOF(volume-of-fluid) 기법을 사용하였다. 수치모형실험이 전체적으로 수리모형실험에서 관측한 파랑의 처오름 및 처내림 현상을 잘 재현하고 있는 것으로 나타났으며, 정량적인 비교를 통해 수치모형의 성능을 검증하였다.

• Journal of Astronomy and Space Sciences
• /
• 제2권1호
• /
• pp.19-33
• /
• 1985

3체 문제에 있어서 Newtonian 운동방정식의 특이접은 수치적분에 의한 특수해를 구하는 과정에서 정확도를 떨어뜨리고 computer 사용시간을 증가시킨다. 이러한 특이점은 Newtonian 운동방정식의 독립변수와 좌표축을 변환하는 정칙화(regularization)과정을 통해 제거할 수 있다.이 논문에서는 정칙화된 Newtonian 운동방정식을 $5^{th}$ -order Runge-Kutta 방법으로 특수해를 구하기 위해 BASIC 언어로 작성한 computer program 으로 적절한 초기위치와 초기속도를 가정하여 시간에 따른 3체의 위치와 속도를 계산하였다. 그 결과, 이러한 3중성계의 진화는 결국 분열되어 3체 중 1개는 쌍곡선 궤도를 그리면서 계를 탈출하고 나머지 2개는 연성계를 형성하게 되었다. 이는 연성계의 기원을 설명할 수 있는 하나의 방버이 되지 않을까 생각한다.

• 한국전산구조공학회논문집
• /
• 제16권3호
• /
• pp.311-319
• /
• 2003

본 연구에서는 지수형 함수로 표현된 비선형 충격력을 받는 시스템의 동적 응답 정밀도를 향상시키고자 선형 충격력의 경우와 비교ㆍ검토하였다. Newmark 계열의 중앙법 알고리듬이 사다리꼴법과 같이 충격력이 없는 경우 모멘텀과 에너지 보존을 성립하도록 유도되었다. 중앙법, 사다리꼴법, 시간구간 종점 평가법((n+1)점 방법)을 선형 충격력에 적용하면 적분간격의 크기에 상관없이 보존성질을 만족하나, 비선형 충격력의 경우 모멘텀과 에너지 보존 상수값이 과소 또는 과대평가 되어졌다. 이러한 오차를 제거하고자 시간간격을 늘리면서 평가함수의 개수를 최소로 하는 다단계 방법중의 하나인 Simpson 1/3법을 사용하여 보존상수값의 정밀도를 향상시켰다. 아울러 유한회전을 포함한 유한운동을 해석할 때에도 제안된 알고리듬이 확정ㆍ적용될 가능성을 보여주고 있다.

• 한국지반공학회논문집
• /
• 제29권8호
• /
• pp.65-73
• /
• 2013

지중 열교환 시스템은 지속적인 에너지 효율의 개선으로 공간 냉난방을 위한 친환경적 에너지 기술로 주목받고 있다. 지중에 매설된 파이프는 내부 유체 순환을 통하여 인접한 지반과 열적 상호작용으로부터 직접적인 열에너지 교환을 수행한다. 하지만, 파이프의 수치모델링에서 열-수리가 연관된 난류해석과 파이프의 긴 세장비에 의한 메쉬사이즈의 부적합성은 열교환 시스템의 적절한 수치해석을 어렵게 하고 있다. 본 논문에서는 파이프 내부 유체흐름에 대한 에너지 보존의 법칙을 적용하여 지배방정식을 유도하였으며, Galerkin수식화와 시간적분을 통하여 열-수리 연동일차원 파이프 요소를 개발하였다. 그리고 제안된 파이프 요소를 기 개발된 다공질 재료를 위한 열-수리-역학(Thermo-Hydro-Mechanical) 해석을 위한 유한요소 프로그램과 결합하였다. 개발된 요소를 이용한 수치해석 결과는 열응답 시험(Thermal Response Test) 결과로부터 주위지반의 유효 열전도도를 평가하기 위하여 사용하는 선형 열원 모델이 인접 파이프간의 열적상호작용과 파이프의 단부효과에 의하여 지반의 열전도도를 과다 평가하는 것으로 보여주었다. 따라서 열응답 시험 해석 결과에 대한 역해석을 적용하여 최적의 수렴성을 보여주는 변환행렬을 제시하였다.

• 한국전자파학회논문지
• /
• 제10권4호
• /
• pp.550-563
• /
• 1999

본 논문은 시변 맥스웰 방정식에 기초한 웨이브릿-갤러킨 설계를 제안하였다. 두 개의 모멘트 함수가 0 이 되는 Daubechies 웨이브릿 함수를 기저함수로 전개하고 Yee가 제안한 Leap-frog 접근법을 적용하였다. D Daubechies 웨이브릿의 변위된 보간 특성을 이용하여 적분이나 매체 연산자에 대한 부가적인 행렬이 필요없 는 방정식을 유도하였다 안정화 조건을 유도하고 분산특성을 분석한 후 유한차분 시간영역법과 다해상도 시 간영역법의 결과와 비교하였고. 분산특성의 분석을 통해 기저함수의 정규성(Regularity)과 받침폭(Support width) 사이의 균형을 확인했다. 기저함수가 단 2개의 0이 되는 웨이브릿 모멘트 함수를 가지지만. 이는 수치 해석 상에서 무시할 수 있는 분산 오류를 수반하였고, 컴팩트 받침(Compact support)에 의해 노드 당 적은 수의 계수만이 고려되었다. 제안된 설계의 저장계수의 효율, 실행 시간의 감소와 정확도를 균일 공진기와 비 균일 공진기의 공진주파수 해석을 통해 검증하였다.