• 한국해양공학회지
• /
• 제33권1호
• /
• pp.17-25
• /
• 2019

A numerical method is suggested for evaluating the singular integral of curved panels in the higher-order boundary element method. Two-step mapping procedures that are significantly related to the physical properties of singular behaviors were developed and illustrated. As a result, the singular behaviors were significantly alleviated, and the efficiency and robustness of the present method for tangentially and axially deformed elements were proven. However, inaccuracies and numerical instabilities of twisted elements were discovered as a result of nonlinearities.

• 대한기계학회논문집
• /
• 제14권6호
• /
• pp.1408-1416
• /
• 1990

본 연구에서는 과도 연성 열탄성문제의 해를 구할때 사용되는 직접시간적분방 법과 Laplace 변환방법은 상호 장 단점을 가지고 있다. 각 방법들의 장단점은 서로 배타적이므로 서로의 장점을 살리는 수치방법이 필요하다. 그런데, 대부분의 과도 열탄성문제는 급격한 온도변화로 인한 물체의 변형에 관심이 있기때문에 이 형태의 문 제를 효율적으로 다루는 데 주안점을 두고 본 연구를 수행하였다. 유도된 유한요소 방정식은 결국 열탄성 지배 방정식 중 열전달방정식인 에너지보존식은 Gurtin의 범함 수로부터 유도된 원래의 형태를 사용하나 수치적 안정성(numerical stability)을 보장 하기 위하여 운동방정식은 시간에 대한 2차미분 형태로 수정하였다. 에너지보존식은 시간에 대한 합성적분(convolution)형태로 표현되므로 온도의 시간미분항이 소거되므 로 경계에서의 급격한 온도변화로 인한 수치 해석적 문제점은 간단히 해결된다. 그 러므로, 제안된 수치해법은 직접시간적분방법의 일종이나 결과식인 유한요소방정식은 기존의 문헌들과 상당한 차이가 있다. 과도 연성 열탄성해석을 위한 새로운 근사수 치해법의 장점을 이론적으로 설명하기보다 수치계산면에서의 안전성, 정확성 및 효율 성이 있음을 증명하기 위하여 이미 발표된 문헌들에서 다룬 예제를 선정하여 해석결과 를 비교하였다.

• 자원환경지질
• /
• 제27권2호
• /
• pp.191-199
• /
• 1994

층상대지속에 있는 3차원 물체에 의한 지자기 지전류 (MT)응답을 계산하기 위하여 적분방정식을 이용한 수치모델링법을 개발하였다. 이 방법에서는 3차원 이상체를 몇 개의 세포로 분할하여 펄스기조함수로 근사할 수 있는 전류분포로 치환한다. 층상대지를 표현하는 전기 텐서그린함수를 쓰면 행렬방정식을 유도할 수 있으며 이를 각 세포의 벡터전류에 대하여 푼다. 결국 미지의 산란장은 산란전류에 관한 전기 및 자기 텐서그린함수를 적분함으로써 얻어진다. 지표면 근처에 3차원 전도성물체가 존재할 때 2차원의 TE모드 모델링으로는 깊은 곳에 가짜의 저 비저항을 가정해야 한다. 이는 TE모드 모델링에서는 경계면 전하의 영향을 고려할 수 없기 때문이다. 그러나 긴 3차원 직방체의 가운데를 가로지르는 단연은 2차원 TM모드 아르고리즘으로 정확히 모델화할 수 있으며, 이는 정식화과정에서 경계면 전하가 고려되어 있기 때문이다. 다중 주파수에 관한 수치계산 결과 겉보기 비저항과 위상은 상호보완적인 변수라는 것이 밝혀졌다. 따라서 이들 변수는 주파수영역 MT 해석시 함께 취급되는 것이 바람직하다.

• 한국전산구조공학회논문집
• /
• 제23권3호
• /
• pp.331-340
• /
• 2010

적응적 hp-세분화 기법과 그 기법의 효과적인 구성방법을 포함한 새로운 적응적 유한요소 알고리즘의 기초이론 및 적용이 이 연구를 통해 제시되었다. 적응적 hp-세분화 기초의 유한요소기법은 적분형 르장드르 형상함수와 요소별로 불균등한차수의 분배 및 비정형적인 절점연결과 관련된 연속조건을 만족시킬 수 있는 제약조건을 필요로 한다. 따라서 요소간의 접합부분에서 적응적 hp-유한요소망의 연속성이 중요한 문제로 대두된다. 이러한 문제를 요소경계에 연속성 제약조건을 절점연결 사상행렬을 적용하여 해결하였다. 또한, 적분형 르장드르 형상함수의 계층성질을 이용하여 제시된 알고리즘의 효율적 정식화 방안을 제시하였다. 간단한 캔틸레버문제가 h-세분화, p-세분화 그리고 hp-세분화 방법에 의해 계산되었다. hp-세분화의 결과는 다른 방식의 세분화에 비해 보다 빠른 수렴성을 보여 주는 것이 확인되었다. 그러므로 제시된 hp-세분화 알고리즘은 실제문제에 효율적으로 적용될 수 있을 것으로 생각된다.

• 대한조선학회논문집
• /
• 제36권2호
• /
• pp.22-28
• /
• 1999

복합원형 실린더군의 미소진폭 동위상 진동유동에 대해 수치해석적 및 실험적으로 연구하였다. 유동장은 내부와 외부로 분리하여 해석되어진다. 1차 비선형해는 스톡스 진동해로서 해석적으로 구해지며, 경계층 외부의 2차 정상유동은 경계적분법(Panel Method)와 유한체적법(Finite Volume Code)에 의해 구하였다. 수치해석해와 정성적 정량적으로 비교하기 위하여, 전유동장 측정기법의 하나인 입자영상유속계(PIV)기법을 사용하였으며, 계산시간 단축을 위해 영상압축 및 코드화기법에서 개발된 삼단계추적 알고리즘을 도입하였다. TLP등 저 댐핑 시스템에 있어서 작지만 0 아닌 항력계수의 정확한 시스템 거동 및 공진에 의한 파괴문제에 있어서 중요한 과제이다. 본 논문에서는 실린더의 갯수와 배열, 간격의 변화에 따른 항력계수 값들을 계산하여 도표화하였다.

• 한국해안해양공학회지
• /
• 제6권4호
• /
• pp.413-420
• /
• 1994

하천수 플룸(plume)의 퍼짐을 다루는 역학은 플룸의 경계면이 시간과 공간에 따라서 변하기 때문에 자유경계조건의 문제(free boundary problem)로 다루어야 하는 대단히 복잡한 비선형 문제이다. 더욱이 플룸경계를 통한 주변수의 혼합까지 고려할 경우 그 복잡성은 한층 더해진다. 이러한 비선형성과 복잡성을 피하는 기법의 하나가 적분해석법인 바, 본 논문에서는 하천수 플룸의 흐름축에 수직한 횡방향 및 수심방향에 대하여 기본방정식들을 적분함으로 3차원 문제를 1차원 문제로 치환하는 적분해석법을 사용하였다. 다만 이 일이 가능하기 위해서는 유동변수들(유속, 밀도 등)의 횡방향 및 수심방향의 분포함수가 알려져 있음이 전제되어야 하는데 유속의 축방향성분 및 플룸과 주변수 간의 밀도차가 상기 두 방향에 대해서 가우스(Gauss)분포를 갖는다는 잘 검증된 가정을 활용하였다. 그리고 이 가정에서 플룸의 횡방향 유속을 도출해낸 본 연구자들의 기발표된 논문의 결과도 활용하였다. 결과로 얻어진 연행(entrainment)효과까지 포함한 방정식들을 Runge-Kutta 수치해석법을 사용하여 풀었다. 그리하여 하천수 풀룸의 3차원적 해석을 쉽게 수행할 수 있는 수치해석기법을 얻어냈다.

• 한국전산구조공학회논문집
• /
• 제32권4호
• /
• pp.215-222
• /
• 2019

복합재료의 열전달 문제는 일반적으로 만족시켜야 하는 보존방정식과 경계조건 외에 추가적으로 만족시켜야 하는 계면경계조건의 존재로 인해 새로운 수치기법의 개발에 어려움이 있다. 계면경계조건이 미분방정식의 해에 불연속성을 유발시키기 때문에 이것을 적절하게 처리할 수 있는 특수한 함수의 도입이 필요하며, 이산화를 통한 계 방정식의 구성도 쉽지 않다. 본 논문에서는 계면경계의 불연속성을 모사하는 특수함수를 포함하면서 계면경계조건을 항상 만족시킬 수 있도록 계면경계식 자체를 매입한 미분근사식을 제안하고, 불연속 재료상수를 갖는 열전달 문제를 무요소 강형식으로 이산화한 이동최소제곱 차분법을 제시한다. 개발된 수치기법은 기존의 수치기법들과 달리 수치적분과 계면경계조건을 만족시키기 위한 별도의 구속 방정식이 필요없으며, 빠르고 정확하게 이종재료 열전달 문제의 수치해를 구해준다. 개발된 수치기법으로 다양한 복합재료 열전달 문제를 해석하고 오차의 수렴률을 조사한 결과, 높은 정확성과 계산 효율성을 갖는다는 것을 확인할 수 있었으며, 특히, 계면경계가 기하학적 특이성을 나타내는 문제에서도 우수한 성능을 발휘하는 것을 보였다.

• 한국추진공학회지
• /
• 제22권1호
• /
• pp.36-44
• /
• 2018

고체 로켓 추진기관 노즐의 내열재로 사용되는 탄소/페놀릭 복합재료의 열반응 수치해석을 수행하였다. 본 논문에서 탄소/페놀릭 재료의 열반응 해석은 (1) 로켓 노즐벽에서 대류열전달계수를 구하기 위한 연소가스의 경계층 적분방정식 수치해석과 (2) 삭마두께, 숯깊이 및 온도를 계산하기 위한 탄소/페놀릭의 열반응(열분해, 삭마)을 고려한 1차원 열전도 해석으로 구성된다. 시험결과와 해석결과를 비교 분석하였으며, 목삽입재 좌우 인접 부위를 제외하고 잘 일치하는 것을 확인 할 수 있었다.

• 한국해양공학회지
• /
• 제12권1호
• /
• pp.85-98
• /
• 1998

본 논문의 수치해법은 경계치문제를 풀기 위하여 코시이론(Cauchy's theorem)을 사용하였다. 경계치문제는 완전한 물체표면조건과 자유표면조건을 만족시키는 초기치문제로 귀결된다. 현 수치해법에서 무한영역은 수치계산 영역인 비선형 영역과 선형 자유표면조건을 만족하는 선형영역으로 나누어진다. 선형영역의 해는 과도 그린(Green)함수를 사용하여 정합조건을 부과함으로써, 수치계산은 비선형 영역에서만 수행된다. 본 논문에서 저자는 수치계산 영역에서 코시이론을 사용하여 적분방정식을 도출하였고, 무한영역의 해는 정합면에서 과도 그린함수를 사용하여 표현하였다. 본 수치계산에서 자유표면에 요소 재분배법을 적용함으로써 쇄파현상에 대해서도 안정적인 수치해석을 할 수 있었다. 본 논문에서 개발된 수치방법을 적용한 문제는 다음과 같다. 첫째는 자유표면에서 실린더가 강제동요하는 경우에 자유표면형상과 힘을 계산하여 이전의 실험치 및 계산치와 비교하였다. 두번째로는 실린더가 자유수면하에서 일정한 속도로 항주하는 경우에는 조파저항과 양력을 계산하여 고차 스펙트럴법과 비교하였다.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2004년도 학술발표회
• /
• pp.1376-1380
• /
• 2004

본 연구에서는 홍수범람지도 제작을 위해 사용되는 수치모형의 검증을 수행하고자 한다. 검증할 수치모형은 스위스의 Beffa에 의해 개발된 FLUMEN(FLUvial Modeling ENgine)으로서 수심 적분된 2차원 비선형 천수방정식에 불규칙 심각망을 이용한 유한체적법(finite volume method)이 적용된 수치모형이며 스위스, 독일, 오스트리아 등에서 홍수범람해석에 사용된 바 있는 모델이다. FLUMEN 모형의 검증을 위하여 범람해석시 가장 중요한 문제인 이동경계조건(moving boundary condition)을 포함하고 있는 원형섬에서의 고립파에 처오름높이를 계산하여 수리모형실험 결과와 비교한다. 수리모형실험은 미국 육군 공병단 산하의 해안공학연구소(CERC, Coastal Engineering Research Center)에서 수행되었으며(Liu 등, 1995) 수조의 중앙에 높이 0.625m, 해저지름 7.2m, 경사각 $14.04^{\circ}$의 원형섬이 위치한다. 본 연구 결과, 수치해석으로 계산된 섬에서의 실제 처오름높이와 입사파의 파고의 비(R/A)는 수리모형실험의 결과와 어느 정도 일치하였다.