• 한국해안·해양공학회논문집
• /
• 제20권5호
• /
• pp.461-471
• /
• 2008

WKB 방법에 의한 일차원 완경사 파랑식의 해석해를 유도하였고 이는 Porter(2003)의 해와 유사한 형태를 갖는다. 적용적 측면에서 본 해석해는 관련 수치해에 견줄 만큼 해석상의 일반성을 갖는다. 유도과정에서 해면함수로 표현된 굴절 방정식의 해도 얻었다. Bremmer 방법을 이용한 본 해석해에 대한 수치계산 결과를 제시하였고 이들은 기존 결과와 일치한다.

• 대한기계학회논문집B
• /
• 제35권7호
• /
• pp.723-733
• /
• 2011

본 연구에서는 비정렬격자계와 체적포착법을 사용하여 표면장력이 지배적인 다상유동의 수치해석 방법을 제시하였다. 먼저 표면장력에 대한 CSF(Continuum Surface Force) 모델을 비정렬격자계에 적용할 수 있도록 수치해석 방법을 확립시켜 Myong(2009)이 개발한 비정렬격자계와 체적포착법을 사용한 수치 해석코드에 삽입하였다. 테스트 문제로 오직 표면장력만이 존재하는 평형상태의 정적(static) 액적 및 비평형상태의 동적(dynamic) 액적 문제에 적용하여, 이 해석방법의 유용성과 정확도를 평가하였다. 연구결과, 매끄러운 곡률 계산을 위해 필요한 필터로 본 연구에서 제안한 Laplacian 필터와 함께 CSF 모델로는 밀도보정(density-scaled)한 CSF 모델이 예측성능이 우수한 것으로 나타났다. 또한 표면장력 계산을 위한 이 모델을 채용한 본 수치해석방법은 표면장력이 지배적인 다상유동인 평형상태의 정적 액적 및 비평형상태의 동적 액적 문제 모두에 대해 정확성과 유용성이 입증되었다.

• 터널과지하공간
• /
• 제21권3호
• /
• pp.216-224
• /
• 2011

본 연구에서는 Class II 거동에 대해 입자결합모델을 이용하여 수치해석적인 방법으로 살펴보았으며, 횡방향 변형률 제어 압축시험을 수치해석적으로 모사할 수 있는 방법을 제시하였다. 수치해석에서 사용된 미시변수는 스웨덴 Aspo Hard Rock Laboratory에서 수행한 일축압축시험을 이용하여 결정하였다. 제시된 수치해석 기법을 이용해 Aspo 암석의 Class II 거동을 효과적으로 모사할 수 있었으며, 수치해석의 결과는 실험실 시험 결과와 잘 일치하였다.

• 한국지반공학회논문집
• /
• 제17권6호
• /
• pp.245-255
• /
• 2001

기존 개탁식터널은 굴착경사와 관계없이 터널상부에는 되메움지반의 자중이 그리고, 터널벽체에는 정지토압이 작용하는 것으로 가정하여 설계되고 있다. 그러나 개착식터널에 실제 작용하는 토압은 개착규모에 따라 다르며 이에 따라 터널라이닝의 구조적거동 또한 다르게 나타난다. 본 연구에서는 수치해석적방법을 사용하여 굴착경사가 터널라이닝에 미치는 영향을 다음과 같이 조사하였다. 첫째로, 지반범용수치해석프로그램인 FLAC2D를 사용하여 개착공법으로 시공된 기존 터널의 변형거동을 모사할 수 있는 수치해석방법을 활용하여 현장의 계측결과와 비교하므로써 수치해석방법의 타당성을 검증하였고, 둘째로, 동일한 수치해석기법을 적용하여 굴착경사가 $30^{\circ}\;, 456{\circ},\; 60^{\circ},\; 75^{\circ}%인 개착식터널에 작용하는 토압 및 변위와 터널라이닝에 발생하는 단면력을 조사하였다. 수치해석에 사용된 개착식터널은 2차로 도로터널이며, 터널라이닝과 되메움지반 사이의 마찰력을 고려하기 위하여 접합면 요소를 사용하였다. 수치 해석결과, 되메움지반과 터널라이닝 및 굴착사면 사이의 마찰력에 의한 되메움지반 내부의 아칭현상으로 인하여 굴착 경사가 클수록 터널벽체에 상대적으로 작은 토압이 작용하였으며, 이에 따라 터널라이닝의 변형, 모멘트, 전단력이 증가하였음을 알 수 있었다.

• 한국방재학회 논문집
• /
• 제8권1호
• /
• pp.133-138
• /
• 2008

외해에서 내습하는 파랑 자료를 분석하는 일은 연안에서 발생하는 문제를 해결함에 있어 기본이 되기 때문에 매우 중요하다. 파랑을 해석하는 방법에는 크게 수치 모델을 이용하는 방법과 해석 해를 이용하는 방법이 있다. 수치 모델의 경우, 다양한 지형과 파랑 조건에 대해 적용할 수 있다는 장점이 있지만 수치 오차를 고려해야 하는 번거로움이 있다. 반면, 해석 해의 경우 수치 오차 없이 빠르고 정확하게 해를 구할 수 있다는 장점이 있지만 특정한 지형 및 파랑 조건에서만 성립한다는 단점이 있다. 본 연구에서는 수치적인 기법과 해석적인 접근을 혼합하여 수치 오차를 최소화시키면서 다양한 조건에 적용이 가능한 완경사 방정식의 해를 유도하였다. 유도된 해를 기존의 수치 해와 비교한 결과 매우 잘 일치한다는 알 수 있었다.

• 대한기계학회논문집B
• /
• 제39권4호
• /
• pp.359-370
• /
• 2015

본 연구는 다양한 조건하에서 모세관력 불균형에 의해 구동되는 수평 표면 위의 액적 거동을 수치해석적으로 연구한 것이다. 액적 거동은 자체개발 코드(PowerCFD)를 사용하여 수치해석하였다. 수치해석에 사용된 코드는 보존적인 압력기반 유한체적방법에 기초한 비정렬 셀 중심 방법 및 VOF 방법에 체적포착법인 CICSAM을 채용하고 있다. 상세한 액적 거동이 다양한 초기 액적형상, 접촉각 및 표면장력(또는 Bond 수)의 조건하에서 얻어졌다. 또한 액적 이송 메커니즘이 액적 형상에 대한 수치해석 결과로부터 검토되었다.

• 한국음향학회지
• /
• 제25권5호
• /
• pp.229-238
• /
• 2006

본 논문에서는 수치 영역의 포물선 지배 방정식의 근사 차수와 수치 영역 경계의 비국소적 경계 조건의 근사 차수가 서로 다를 때 음파 해에 미치는 영향을 해석적으로 보였다. 우선 평면파 분석법을 이용해 비국소적 경계 조건을 반 무한 매질 영역으로 변환했다. 그리고 실제 수치 영역과 반 무한 매질 영역의 경계에서 해석적 반사 오차를 유도했다. 지배 방정식과 비국소적 경계 조건의 해석적 오차가 간단한 대수 식으로 표현 가능한 경우에 대해서는 대수적인 오차식을 유도하고 그 경향을 고찰했다. 지배 방정식이 일반적인 고차 포물선 방정식일 때는 대수적인 오차 식은 보다 복잡하게 표현되며 수치적 방법을 이용해 그 특성을 고찰했다. 최종적으로 지배 방정식의 차수에 따른 비국소적 경계 조건의 정밀도를 유도하고 해석적 반사 오차의 전반적인 특성에 대해 논의했다. 본 연구의 핵심 공헌은 포물선 방정식과 비국소적 경계 조건의 근사 차수가 다를 때 해석적 오차 추정 방법과 사용한계를 제시했다는데 있다.

• 한국수자원학회논문집
• /
• 제36권5호
• /
• pp.751-760
• /
• 2003

시간 및 공간가중치를 고정하지 않는 Muskingum-Cunge 홍수추적방법에 대한 오차해석을 수행하였다. 오차해석 결과 시간가중치와 공간가중치의 합이 1.0이상인 경우에는 홍수파가 진행하면서 증폭되어 수치해가 발산하였다. 시간가중치와 공간가중치의 합이 작을수록 수치확산이 크게 발생하였다. 격자의 해상도가 낮을수록 수치확산 및 수치진동이 크게 발생하였다. 수치실험과 자연하천에 대한 적용 결과, 공간가중치를 고정하지 않는 경우에는 공간가중치를 0.5로 고정하는 전통적인 Muskingum-Cunge방법보다 첨두의 감쇄가 큰 홍수파 모의에 효과적임을 알 수 있었다.

• 물과 미래
• /
• 제21권3호
• /
• pp.291-291
• /
• 1988

2차원 흐름이 원주주위를 지날 때 발생하는 흐름의 변화가 기본방정식인 연속방정식과 운동량방정식에 의하여 수치적으로 해석된다. 수치해석 과정은 기본방정식에 유함수, 와도 및 흐름의 특성을 나타내는 무차원 매개변수를 도입하여 무차원 유함수-와소수송식을 유도한 후, successive over relaxation scheme과 alternating direct implicit scheme으로 수행된다. 수치실험은 레이놀즈수 125-275를 기존의 수치해석에서는 주로 수치실험 결과와 비교한다. 원주표면의 압력을 구하는 방법에 있어서 기존의 수치해석에서는 주로 방사 운동량방정식만을 사용하였으나, 본 논문에서는 기존의 방법외에 방사 운동량방정식 및 접선 운동량방정식에 의해 압력을 계한하고, 두 값을 비교하여 레이놀즈수에 따른 압력을 구하는 방법을 제시한다. 또한 와도의 분포를 도시하여 원주에 의한 후류의 영향을 받지 않는 외부경계의 한계를 새로이 설정한다.

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

본 논문에서는 진동중인 교량 거더에 작용하는 풍하중을 산정하고 그에 따른 플러터 발생풍속을 예측하기 위하여 분산형 전산환경을 활용한 수치해석 연구를 수행하였다. 분산형 전산환경은 웹 포탈을 기반으로 수치해석 환경을 제공하는 수치풍동 시스템으로서, 전산유체역학(CFD : Computational Fluid Dynamics)에 대한 전문지식이 부족한 사용자들도 격자생성, 수치해석자를 이용한 계산, 가시화 등의 전 과정을 편리하게 수행할 수 있는 차세대 토목분야 연구 환경이다. 본 시스템은 그리드스피어(GfidSphere)를 기반으로 구성되었으며, 기본적으로 사용자 관리, 세션 관리, 그룹 관리, 레이아웃 관리 등을 제공하여 사용자가 포탈을 통해서 다양한 서비스를 쉽게 사용할 수 있는 환경을 구축하도록 도와준다. 수치해석을 위한 유체 지배방정식은 2차원 비정상 비압축성 RANS(Reynolds-Averaged Navier-Stokes) 방정식이며, pseudo compressibility 방법을 적용하였다. 비정상 유동장을 해석하기 위하여 이중시간 전진법(dual time stepping)을 사용하였으며, 수렴가속화를 위해 Multi-grid 기법을 적용하였다. 또한 난류 유동장 해석을 위해서 $k-{\omega}$ SST 난류 모델을 사용하였으며, 난류 천이 과정에서의 유동을 모사하기 위하여 Total stress limitation 방법을 적용하였다. 교량 거더의 연직과 회전방향의 2자유도 움직임을 모사하기 위하여 동적격자 기법을 도입하였다. 교량 거더 주변의 비정상 유동해석 결과를 통해, 거더 표면에서 떨어져나가는 크고 작은 와류의 영향으로 양력 및 모멘트 계수 그래프가 중첩된 진폭과 주기를 갖고 주기적으로 나타나는 것을 확인할 수 있었다. 또한 계산된 비정상 공기력을 적용한 2자유도 플러터 방정식을 통하여 플러터 발생풍속을 산정하였다. 최종적으로 본 연구에서 계산된 결과의 타당성을 검증하기 위하여 수치적으로 구한 플러터 발생풍속과 기존의 실험 및 수치해석 결과를 비교하였으며, 결과는 잘 일치하였다.