• 한국추진공학회지
• /
• 제6권2호
• /
• pp.75-84
• /
• 2002

고체 추진 로켓 내부 연소실의 비정상 유동을 수치적으로 해석하였다. 완전 보존식을 이용하여 2 차원 축-대칭 연소실 안의 연소 불안정을 해석하기 위한 수치 기법을 구성하였는데 비정상 유동을 해석하기 위한 수정된 $\kappa$-$\varepsilon$ 난류 모델이 사용되었다. 이산화한 지배 방정식은 연관된 경계 조건을 포함하여 dual time-stepping 방법으로 시적분 하였다. 정상 상태의 계산을 기반으로 연소실 내의 천이 압력파의 비정상 상태를 수치적으로 모사하기 위하여 압력 펄스 및 압력 변동을 연소실 상단에 부과하였다. 로켓 모터 연소실 내의 다양한 정상 상태 및 비정상 상태의 특성을 계산 및 해석하였다.

• Nuclear Engineering and Technology
• /
• 제19권3호
• /
• pp.198-204
• /
• 1987

발전소 과도현상과 비냉각재 상실사고를 모의할 수 있는 가압경수로발전소 모의코드 MCSIM을 개발하였다. 원자로 냉각재계통은 에너지 방정식과 운동량 방정식을 분리 취급하면서 Drift Flux 2상 유동모델, 적분 운동량 방정식 등을 사용하여 모델링하였다. 증기발생기의 모사는 Pot Boiler 모델을 사용하였고, 2차계통을 위해서는 분리 취급된 정상상태 에너지 방정식과 운동량방정식을 핵출력 계산을 위해서는 점 동특성 방정식을 사용하였다. 현재의 코드성능을 시험하기 위해 완전 냉각재 유동상실사고와 제어봉 집합체 인출 사고를 계산하여 그 결과를 원자력 5/6호기 최종 안전 보고서의 결과와 비교하였다.

• 한국지진공학회논문집
• /
• 제3권2호
• /
• pp.87-96
• /
• 1999

시간영역에서의 구조물의 해석은 직접적분법과 모드중첩법에 의하여 구해질 수 있다 그 중에서도 모드중첩법에 의한 해석방법은 몇가지 저차 진동모드를 사용하여 비교적 정확한 해를 구할 수 있기 때문에 동적해석에 널리 사용되고 있다, 그러나 비선형해석에서는 각 부재들의 상태에 따라 강성이 달라지므로 고유 진동모드를 정의할수 없거나 변화되는 강성에 따라 고유진동 모드를 지속적으로 다시 구하여야 하는 불편 있으므로 모드 중첩법을 이용한 비선형해석은 완전탄소성모델 등 극히 제한된 조건에서만 실행이 가능하였다 본논문에서는 강성행렬을 각 부재별로 분리시키고 비선형복원력과 초기선형복언력과 초기선형복원력의 차이를 하중항에 반영시킴으로써 모드중첩법을 이용하여 비선형 해석은완전탄소성모델 등 극히 제한된 조건에서만 실행이 가능하였다 본 논문에서는 강성행렬을 각 부재별로 분리시키고 비선형 복원력과 초기선형복원력의 차이를 하중항에 반영시킴으로써 모드중첩법을 이용하여 비선형해석이 가능한 방법을 제시하고자 한다. 특히 각 부재 강성을 각 부재의 상대변위의 함수로 나타냄으로써 연속적인 계산이 가능하게 하였다 본 논문에서 제시된 방법은 전단보모델에 적용하였으며 모드 개수를 변화시켜 지진하중에 의한 최대변위를 계산하여 이를 직접적분버에 의한 결과와 비교하였다.

• 한국지반공학회논문집
• /
• 제18권4호
• /
• pp.215-228
• /
• 2002

동반논문(오세붕, 2002)에서는 전체 변형도 영역의 거동을 모델할 수 있는 GUX( 구성모델을 ABAQUS 코드(Hibbit 등, 2001)에 구현하였다. 먼저 정확도 해석을 통하여 동반논문에서 구현한 내재적 응력적분 기법이 해의 정확도를 확보함을 알 수 있었다. GUX 모델은 전응력 개념에 의거한 비등방경화 탄소성 구성모델이므로 완전 배수 또는 비배수 삼축시힘시 응력-변형도 관계를 입력하여 대상문제를 배수 및 비배수조건에 따라 해석할 수 있다. 본 연구에서는 연약지반 및 풍화토 지반상 성토시 안정성 문제와 지반-말뚝계의 축하중 재하시 거동에 대한 예제 해석을 수행하였다. 그리고 대변형 해석을 통하여 기하학적 비선형성을 고려하였고 전체 변형도 거동을 합리적으로 모델하는 GUX 모델을 이용한 결과를 Mises 모델 결과와 비교하였다.

• 한국방재학회:학술대회논문집
• /
• 한국방재학회 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) 기법을 사용하였다. 수치모형실험이 전체적으로 수리모형실험에서 관측한 파랑의 처오름 및 처내림 현상을 잘 재현하고 있는 것으로 나타났으며, 정량적인 비교를 통해 수치모형의 성능을 검증하였다.

• 한국전산구조공학회논문집
• /
• 제18권4호통권70호
• /
• pp.395-404
• /
• 2005

본 논문에서는 동적하중 재하 시 공간 강뼈대 구조물의 비선형 시간이력 해석에 대한 간편하고 효율적인 수치해석 기법을 제시하였다. 뼈대의 강성행렬에 안정함수를 적용함으로서 기하학적 비선형성 즉 P-$\delta$ 및 P-$\Delta$ 효과를 고려하였다. 접선계수 개념 및 New Orbison 완전 소성면에 기초한 연화소성힌지를 적용하여 축방항력에 의한 부재의 점진적 항복 및 휨에 의한 부재의 부분 소성화를 고려하였다. 증분형태로 정식화된 뼈대 구조의 동적평형방정식의 해를 수치적으로 구하기 위하여 적분형태의 평균가속도법을 사용하여 컴퓨터 프로그램을 개발하였다. 검증 예제에 대하여 보요소 모델을 사용한 ABAQUS 및 개발된 프로그램에 의한 해석결과를 비교함으로서 개발된 프로그램의 정밀도와 해석시간에 대한 효율성을 입증하였다.

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

최근 우리나라에서는 기상이변과 기후변화에 의한 국지성 집중호우의 발생으로 인해 인명 및 재산 피해가 증가하는 추세이다. 따라서 이러한 기상현상을 좀 더 정확하게 예측하고 이를 대응하고자 악기상 모형의 개발과 구축 및 활용에 대한 연구들이 활발하게 진행 중에 있다. GCM이 제공하고 있는 많은 유용한 정보에도 불구하고 대부분의 모델이 시 공간 분해능과 물리 과정의 한계점으로 인해 지역적인 기후 특성이나 변화를 예측하기에는 많은 문제점들이 나타나고 있다. GCM의 한계점을 극복하기 위한 방법으로 세밀한 규모의 기후 정보를 얻기 위해 복잡한 지형과 해안선, 호수, 식생, 지표특성과 같은 아격자 규모의 강제 효과를 반영할 수 있는 고해상도 지역 기후 모델(Regional Climate Model, RCM)의 필요성이 제기되었다. 본 연구에서는 전지구 20km 격자자료를 입력장으로 하여 8km 격자로 한반도를 포함하는 도메인에 대해 비정역학 완전 압축성 중규모 모델인 WRF를 이용하여 상세예측자료를 생산하고자 하였다. 강수 예측의 경우 돌발적으로 발생하는 경우가 많아, 이를 예측하기 위해서는 상세한 강수량 정보를 빠른 시간 내에 정확히 제공할 수 있는 모델을 사용하여야 한다. 강수의 경우 온도와는 달리 공간적 편차가 매우 커 지역적으로 정확한 강수량을 예측 하는데 어려움이 있다. 상세강수 예측을 위해 미세 격자 규모의 비 정역학 모형을 사용할 경우 계산양이 매우 늘어나기 때문에 장시간의 모형 적분 시간뿐 아니라, 상당한 컴퓨터 자원을 필요로 하므로 이에 대한 대안으로 지형효과를 포함한 강수량 진단 모형인 QPM(Quantitative Precipitation Model)을 사용하였다. 최종적으로 한반도의 복잡한 지형적 영향을 반영하기 위해 1 km의 수평해상도를 가지는 고해상도 강수량 진단 모형(QPM)과 상세한 지리적, 공간적 분석을 할 수 있는 ARCGIS를 이용하여 한반도 도별 상세 강수자료를 생산하고자 한다.

• 에너지공학
• /
• 제16권1호
• /
• pp.46-52
• /
• 2007

특정한 방향에 대해 방향중성자속(angular neutron flux)을 정의하는 방향차분 방정식(discrete-ordinates or $S_{N}$ equation)과 달리 방향변수를 구분된 방향영역에 대하여 적분한 값을 사용하고, 해당 방향영역 내에서 방향중성자속이 일정하다고 가정하는 영역상수법(piecewise-constant method)을 개발하였다. 기존 방향차분법과 본 연구에서 개발된 영역상수법을 1계 수송방정식(1'st-order Boltzmann transport equation)과 2계 우성 방정식(even-parity equation)에 적용하여 방향차분 방정식인 $S_{N}$ 방정식과 유사 방향차분방정식($S_{N}$-like equation)인 $PC_{N}$ 방정식을 유도하였다. 우성 방정식에 영역상수법을 적용한 경우 기존 방향차분법의 단점인 광첨두 현상(ray effect)이 현저히 감소함을 확인하였는데 이는 우성 방정식의 혼합 미분항의 기여도가 작아지기 때문인 것으로 판단된다. 이러한 이론은 우성 방정식에서 혼합 미분항이 제거된 단순우성 방정식(simplified even-parity equation)을 사용하는 경우 광첨두 현상이 완전 제거 또는 극단적으로 감소되었던 이전의 결과를 이론적으로 설명한다.

• Composites Research
• /
• 제16권6호
• /
• pp.23-32
• /
• 2003

복합 신소재의 건설 분야 적용은 확대되고 있는 추세이고, 높은 비강도 비를 지닌 복합 신소재 교량 바닥 판의 개발은 국내외에서 활발히 진행되고 있다. 본 연구의 목적은 이동 하중을 받는 비등방성 복합재료 적층판의 동적 응답을 정식화하고, 유한요소법을 사용하여 적층에 따른 동적거동 특성을 분석하는 것이다. 수치 해석 모델에 대해서 이동하중의 속도를 증가시키면서 동적 확대 계수를 계산하였다. 또한 적층형식 및 순서, 섬유 보강 각도 등의 변화에 따른 동적 거동 특성을 분석하였다. 본 연구 프로그램의 타당성을 확보하기 위해 휨과 자유 진동 해석에 관한 기존 문헌 결과와 비교하여 검증하였다. 또한 이동 하중에 의한 동적 해석에 대해 모우드 중첩법과 Newmark 직접 적분법을 사용하였다. 이러한 이동 하중과 적층 수, 적층 순서 및 섬유 보강 각도에 따른 수치 해석 결과는 완전 복합 신소재 교량 바닥 판을 개발하는데 있어서 중요한 기초 자료로 제시될 수 있을 것이다.

• 자원환경지질
• /
• 제48권6호
• /
• pp.451-465
• /
• 2015

주변 해역을 포함한 한반도 일원에서 측정된 중력자료로부터 상층중력의 고도이상(free-air anomaly)을 계산하였다. 주변 영역에서는 인접국가가 발표한 중력자료가 있는 경우 발표된 자료를 이용하였으며, 없는 경우 EGM2008(Earth Gravitational Model 2008)로부터 계산한 고도이상을 이용하였다. 중력의 상향연속은 Dragomir가 제안한 방법으로 계산하였다. 상층중력 고도이상 계산의 정확성과 계산 속도를 고려하여 적분반경은 계산 고도의 10배로 하였다. 적분에 필요한 측지선의 거리는 Bowring이 개발한 공식을 사용하였다. 위도 $33^{\circ}N{\sim}43^{\circ}N$, 경도 $124^{\circ}E{\sim}131^{\circ}E$에서 계산된 고도이상은 고도 1 km에서 -41.315에서 189.327 mgal까지 변화하고 표준 편차는 22.612 mgal이다. 고도 3 km에서는 -36.478에서 156.209 mgal까지 변화하고 표준 편차는 20.641 mgal이다. 고도 1,000 km에서는 3.170에서 5.864 mgal까지 변화하고 표준 편차는 0.670 mgal이다. 3 km 고도에서 계산된 고도이상을 같은 높이에서 측정한 항공 중력 고도이상과 비교하였다. 이들의 rms 오차는 3.88 mgal로 나타났다. 항공 중력 측정 교차점오차가 2.2 mgal 임을 고려하면 이들 오차에 의미를 부여할 수 없으며, 원인으로는 이번 연구에서 발생한 계산상 오차와 함께/또는 발표된 항공중력의 보정오차에 기인하는 것으로 사료된다. 상층중력 고도이상에 완전식으로 계산한 지구타원체 외부의 정규중력을 더하여 상층중력을 예측하였다. 이번 연구에서 국내 최초로 계산한 고도에 따른 상층중력 고도이상은 한반도 일원의 상층중력장을 잘 표현하고 있는 것으로 보이며, 상층중력장은 관성항법장치의 정확도 향상 등에 이용될 수 있을 것이다.