• 제목/요약/키워드: 수치적 근사

검색결과 464건 처리시간 0.035초

다중 경사 MUSCL을 이용한 천수방정식의 수치해석 (A Numerical Analysis of the Shallow Water Equations Using the Multi-slope MUSCL)

  • 황승용;이삼희
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2011년도 학술발표회
    • /
    • pp.158-158
    • /
    • 2011
  • 천수방정식과 같은 쌍곡선형 미분방정식의 불연속 해에 대한 Riemann 해법은, 1950년대 말 공기동역학 분야에서 S. K. Godunov의 선구적인 시도 이후, 다양한 영역에서 성공적으로 적용되고 있다. 당초 제안된 해법은 공간에 대해 1차 정도였으나, 2차의 정도를 얻을 수 있는 기법이 1970년대 말 B. van Leer에 의해 제안되었으며, MUSCL로 불린다. 서로 인접한 격자의 보존변수가 고려된 경사가 도입되어 두 격자에 의해 공유되는 변의 좌 우에서 선형으로 보존변수가 재구축되는 MUSCL은 제한자와 함께 이용될 때, 구조 격자 체계에서 비교적 단순하면서도 효과적인 적용성이 입증되었다. 그런데, 이 기법을 2차원의 비구조 격자 체계에 적용하는 경우, 인접한 모든 격자의 보존변수를 고려한 평면의 경사를 결정해야 하는 어려움이 따른다. 특히, 삼각형 비구조 격자에 적용할 경우 최적의 평면을 결정하기 위해 Green-Gauss 적분식이나 최소-자승법 등을 이용하게 된다. 이에 비해, 2010년 T. Buffard와 S. Clain이 제안한 다중경사 기법은 격자의 각 변에서 경사가 각각 결정되는 방법으로 계산량이 많은 Green-Gauss 적분식이나 최소자승법을 피할 수 있는 장점이 있는 것으로 알려져 있다. 정확해가 알려진 두 경우에 대해 몇 가지 제한자를 적용한 결과를 1차 정도의 해와 함께 비교하였으며, superbee 제한자에 의한 결과가 우수하였으나, 희유파와 충격파가 맞닿는 곳에서 수치 분산이 나타났다. minmod 제한자의 결과가 대체로 무난하였으며, 이를 2차원 댐 붕괴 문제에 적용하여 1차 정도의 해와 비교하였다. 마찰이 없고 초기 수심이 댐 상류에서 10 m, 하류에서 5 m로서 물이 차 있는 경우, 1차 정도의 해에서 나타나는 수치 소산이 2차 정도에서는 발생되지 않았다. 댐 하류에서 초기에 수심이 영으로 바닥이 드러난 경우에서 마찰의 영향을 검토하였다. 마찰이 있는 경우, 마찰 경사 항의 Manning 계수를 0.04로 두었으며, 마찰에 의한 영향이 잘 드러났다. 수심이 50 mm 보다 작은 경우에는 마찰을 적용하지 않았다. 이 연구는 환경부 '차세대 핵심환경기술개발 사업'의 지원에 의한 것이다.

  • PDF

Charge Trap Flash 메모리 소자 프로그램 동작 시 전하수송 메커니즘

  • 유주태;김동훈;김태환
    • 한국진공학회:학술대회논문집
    • /
    • 한국진공학회 2011년도 제41회 하계 정기 학술대회 초록집
    • /
    • pp.273-273
    • /
    • 2011
  • 현재 사용되고 있는 플로팅 게이트를 이용한 플래시 메모리 소자는 비례축소에 의해 발생하는 단 채널 효과, 펀치스루 효과 및 소자간 커플링 현상과 같은 문제로 소자의 크기를 줄이는데 한계가 있다. 이러한 문제를 해결하기 위하여 silicon nitride와 같은 절연체를 전자의 트랩층으로 사용하는 charge trap flash (CTF) 메모리 소자에 대한 연구가 활발히 진행되고 있다. CTF 메모리 소자의 전기적 특성에 대한 연구는 활발히 진행 되었지만, 수치 해석 모델을 사용하여 메모리 소자의 전하수송 메커니즘을 분석한 연구는 매우 적다. 본 연구에서는 수치 해석 모델을 적용하여 개발한 시뮬레이터를 사용하여 CTF 메모리 소자의 프로그램 동작 시 전하수송 메커니즘에 대한 연구를 하였다. 시뮬레이터에 사용된 모델은 연속방정식, 포아송 방정식과 Shockley-Read-Hall 재결합 모델을 수치해석적 방법으로 계산하였다. 또한 CTF 소자 프로그램 동작 시 트랩 층으로 주입되는 전자의 양은 Wentzel-Kramers-Brillouin 근사 법을 이용하여 계산하였다. 트랩 층에 트랩 되었던 전자의 방출 모델은 이온화 과정을 사용하였다. 게이트와 트랩 층 사이의 터널링은 Fowler-Nordheim (FN) tunneling 모델, Direct tunneling 모델, Modified FN tunneling 모델을 적용하였다. FN tunneling 만을 적용했을때 보다 세가지 모델을 적용했을 때가 더 실험치와의 오차가 적었다. 그 이유는 시뮬레이션 결과를 통해 인가된 전계에 의해 Bottom Oxide 층의 에너지 밴드 구조가 변화하여 세가지 tunneling 모델의 구역이 발생하는 것을 확인 할 수 있었다. 계산된 결과의 전류-전압 곡선을 통해 CTF 메모리 소자의 프로그램 동작 특성을 관찰하였다. 트랩 층의 전도대역과 트랩 층 내부에 분포하는 전자의 양을 시간에 따라 계산하여 트랩 밀도가 시간이 지남에 따라 일정 값에 수렴하고 많은 전하가 트랩 될 수록 전하 주입이 줄어듬을 관찰 하였다. 이와 같은 시뮬레이션 결과를 통해 CTF 메모리의 트랩층에서 전하의 이동에 대해 더 많이 이해하여 CTF 소자가 가진 문제점 해결에 도움을 줄 것이다.

  • PDF

유도전동기 회전자 제작시 압입작업 평가 II (Evaluations of Swaging Process for Rotor Core of Induction Motors II)

  • 박상철
    • 한국산학기술학회논문지
    • /
    • 제18권6호
    • /
    • pp.464-469
    • /
    • 2017
  • 유도전동기 회전자는 rotor core를 일정한 높이로 쌓은 후 rotor core slot에 Cu bar를 열박음(shrinkage fit)작업으로 고정한 후 Cu bar 표면을 punch로 길이방향, 원주방향으로 순차적으로 압입(swaging)하여 제작된다. 본 연구에서는 Cu bar 표면을 압입할 때 Cu bar 표면에 발생하는 압입깊이, 잔류응력 분포 그리고 rotor core slot과 Cu bar 접촉면에 작용하는 접촉력(contact force)을 평가하고자 하였다. 한편 실제의 압입작업은 복잡한 3차원 문제이나 수치해석적인 접근을 위하여 단순화된 2차원 cyclic symmetry 해석모델을 개발하였으며 이러한 해석모델을 사용하여 일정한 압입압력에서 4가지 종류의 Cu bar 형상에 대하여 압입량과 접촉력을 평가하였다. 탄소성 수치해석 결과 1) sample Cu bar block model을 사용하여 단순화된 2차원 수치해석 모델을 개발하였고 2) rotor core 압입작업시 일정한 압입압력에서 Cu bar 표면의 압입량은 Cu bar 형상에 크게 관계없이 거의 일정하며 3) 압입작업시 punch에 가해지는 압입압력으로부터 rotor core slot과 Cu bar 접촉면에 작용하는 총 접촉력을 추정할 수 있는 근사식을 제안하였다.

일정진폭하중을 받는 유한 길이 봉의 유한요소해석 (Finite Element Analysis in Finite Length Bar under Constant Amplitude Loading)

  • 황은하
    • 한국전산구조공학회논문집
    • /
    • 제23권5호
    • /
    • pp.525-533
    • /
    • 2010
  • Newmark방법과 같은 직접시간적분법은 시간증분 구간 사이에서 하중이 변하더라도 하중값을 그 시간 구간에서 일정한 하중으로 사용하기 때문에 일정진폭하중과 같은 연속적인 하중함수를 불연속적인 하중함수로 가정하고 수치계산을 수행한다. 따라서 이러한 하중함수의 근사에 따른 오차로 인하여 정확한 수치결과를 계산할 수 없다. 이에 반해, Gurtin의 변분식 에 기초한 유한요소방정식은 하중함수를 시간이력에 대하여 합성적분하여 계산한다. 따라서 시간증분 구간에서 하중이 변하더라도 연속적인 하중함수의 곡선을 따라 가면서 계산하기 때문에 신뢰할 수 있는 수치결과를 구할 수 있다. 본 논문에서는 1차원 막대의 자유단에서 일정진폭하중을 받는 문제를 수치해석하여 Gurtin방법이 Newmark방법 보다 일정진폭하중을 받는 문제에 더 적합한 방법임을 보인다. 또한, Gurtin방법이 일정한 하중을 받는 문제보다 일정진폭하중을 받는 문제에 더 효과적인 방법임을 보인다. Gurtin방법을 FORTRAN으로 프로그래밍하여 해석한 수치결과와 해석용 소프트웨어인 ADINA의 Newmark방법에 의한 수치결과를 비교하여 제시된 수치해의 정확성과 타당성을 검증한다.

일반적인 연결선 구조의 해석을 위한 효율적인 행렬-벡터 곱 알고리즘 (An Efficient Matrix-Vector Product Algorithm for the Analysis of General Interconnect Structures)

  • 정승호;백종흠;김준희;김석윤
    • 대한전자공학회논문지SD
    • /
    • 제38권12호
    • /
    • pp.56-65
    • /
    • 2001
  • 본 논문은 이상적인 균일한 무손실 유전체를 갖는 일반적인 3차원 연결선 구조에서의 커패시턴스 추출 시, 널리 사용되는 일차 대조법(First-order collocation) 외에 고차 구적법을 결합하여 사용함으로써 정확성을 제고하고, 반복적 행렬-벡터의 곱을 효율적으로 수행하기 위한 알고리즘을 제안한다. 제안된 기법은 연결선에서 전기적 성질이 집중되어 있는 코너나 비아를 포함한 경우에 일차 대조법 대신에 구적법을 이용하여 고차로 근사함으로써 정확성을 보장한다. 또한, 이 기법은 경계 요소 기법에서 행렬의 대부분이 수치적으로 저차 계수(low rank)를 이룬다는 회로상의 전자기적 성질을 이용하여 모형차수를 축소함으로써 효율성을 증진한다. 이 기법은 SVD(Singular Value Decomposition)에 기반한 저차 계수 행렬 축소 기법과 신속한 행렬의 곱셈 연산을 위한 Krylov-subspace 차수 축소 기법인 Gram-Schmidt 알고리즘을 도입함으로써 효율적인 연산을 수행할 수 있다. 제안된 방법은 허용 오차 범위 내에서 효율적으로 행렬-벡터의 곱셈을 수행하며, 이를 기존의 연구에서 제시된 기법과의 성능 평가를 통하여 보인다.

  • PDF

PSC-I 거더교의 지진취약도 평가를 위한 HAZUS 방법의 국내 적용성 연구 (Modified HAZUS Method for Seismic Fragility Assessment of Domestic PSC-I Girder Bridges)

  • 서형열;이진학;김두기;송종걸
    • 한국구조물진단유지관리공학회 논문집
    • /
    • 제14권2호
    • /
    • pp.161-170
    • /
    • 2010
  • 교량구조물의 지진재해평가기술을 개발하기 위해서는 전체 교량에 대한 분류체계별 표준교량을 선정한 후, 간편한 정식화 방법을 통해 지진취약도를 분석하고, 이로부터 교량의 근사적인 지진위험도를 평가하는 것이 현실적이라 할 수 있다. 본 연구에서는 HAZUS에서 사용 중인 간편하면서도 실용적인 방법을 미국과 국내의 내진설계수준 등을 고려하여 국내 PSC-I 거더교에 적용 가능하도록 '수정된 HAZUS 방법에 의한 한국형 지진취약도함수'를 개발하였다. 이를 위해 국내 표준교량 형식 중 하나인 PSC-I 거더교에 대해 수치해석적 방법을 적용하여 해석적 지진취약도함수를 구하고, 수정된 HAZUS 방법을 적용하여 지진취약도함수를 구한 후 그 결과를 비교 분석하였다. 수정된 HAZUS 방법에서의 주요 계수는 수치해석적 방법에 의한 지진취약도함수와 가장 유사한 경향을 보이도록 하는 계수를 선택하는 방법으로 결정하였으며, 강도감소계수의 경우 HAZUS에서 제시한 값의 70% 수준을 사용할 때 해석적 결과와 유사한 지진취약도 함수를 구할 수 있음을 알 수 있었다.

민감도분석을 통한 지진하중을 받는 CFRD 정상부 침하량 예측 (Estimation of Settlement on the Crest of CFRD Subjected to Earthquake Loading Using Sensitivity Analysis)

  • 하익수
    • 한국지반공학회논문집
    • /
    • 제23권1호
    • /
    • pp.39-49
    • /
    • 2007
  • 본 연구의 목적은 수치해석을 이용한 정량적 민감도분석으로부터 지진하중을 받는 CFRD 정상부 침하량 변화의 주 영향인자를 찾고, 주 영향인자 변화에 따른 침하량의 변화양상을 고찰하고, 이로부터 지진하중을 받는 CFRD 정상부 침하량을 근사적으로 추정하는 방법을 제안하는 데 있다. 대형삼축압축시험으로부터 얻어진 사력재료 물성을 통계 분석하여 사력재료 물성 통계특성값을 산정하고, 산정된 통계특성값을 이용하여 작성된 27개 해석단면에 대해 해석조건을 달리한(2개 지진파, 각 지진파별 2가지 가속도 크기) 총 108개 CFRD 해석단면에 대한 동적수치해석을 수행, 그 결과를 이용한 정량적 민감도분석을 수행하였다. 민감도분석 결과, 지진하중 작용 시 CFRD 정상부 침하량은 입력 물성 중 전적으로 사력재료의 전단탄성계수에 의존하는 것으로 나타났고 점착력과 마찰각의 영향은 미미하였다. 민감도분석 결과와 전단탄성계수 변화에 따른CFRD 동적해석 결과로부터, 사력재료의 전단탄성계수와 간단한 댐 정보를 이용하여 지진하중 작용시 CFRD 정상부 침하량을 근사적으로 추정하는 방법을 제안하였다.

방사성 폐기물 저장시설에서의 방사선 작업 중 피폭선량 최소화를 위한 최적 작업경로 계획 (Planning of Optimal Work Path for Minimizing Exposure Dose During Radiation Work in Radwaste Storage)

  • 박원만;김경수;황주호
    • Journal of Radiation Protection and Research
    • /
    • 제30권1호
    • /
    • pp.17-25
    • /
    • 2005
  • 원자력 발전의 안정성이 사회적 문제로 제기된 이후 작업종사자에 대한 작업 중 방사선 피폭량에 대한 관심이 높아지고 있다. 현재의 방법에 의하면 방사선 작업 계획의 수립 시 작업 공간 내 선량률이 일정하다는 가정 하에 피폭선량을 예측하므로 작업 경로에 따른 피폭선량의 변화에 대한 고려가 이루어지지 않고 있다. 본 연구에서는 작업자와 선원과의 거리가 고려되 수정된 방사선 피폭량 계산식을 이용하여 방사성 폐기물 저장시설에서의 작업 중 작업경로 변화에 따른 방사선 피폭량을 계산하였다. 이 계산식을 이용하여 주어진 작업 공간과 선원 조건하에서 작업 중 방사선 피폭량을 최소로 하는 최적 작업경로를 탐색할 수 있는 수치해석 알고리즘을 제안하였다. 이를 위하여 2차원 작업공간에서 무한개의 작업경로를 유한개의 경로로 근사화하고 근사화된 모든 작업경로 중 피폭선량이 최소가 되는 작업경로를 탐색하였다. 또한, 3차원 그래픽 기술과 Java 프로그래밍을 이용한 가상작업 시뮬레이션 프로그램을 개발하고 작업 공간의 선량률 가시화 및 가상 작업 시뮬레이션을 수행하여 방사선 작업 계획의 수립을 위한 도구로서의 가능성을 검토하였다. 수치해석 계산과 시뮬레이션 과정을 통하여 최적 작업경로는 작업자와 선원과의 거리를 증가시키고 작업 시간을 단축시킬 수 있는 경로로 제시되었고, 이를 바탕으로 방사선 피폭량은 작업시간뿐만 아니라 작업자와 선원간의 거리에 영향을 받음과 최적화된 방사선 방호를 위해 작업경로가 고려되어야 함을 확인하였다.

$CO_2$에 의한 층류 확산화염의 소화기구에 관한 정량적 분석 (A Quantitative Analysis on the Extinguishing Mechanism of Laminar Diffusion Flames by $CO_2$)

  • 문선여;황철홍;오창보;박원희
    • 한국화재소방학회:학술대회논문집
    • /
    • 한국화재소방학회 2011년도 추계학술논문발표회 논문집
    • /
    • pp.378-381
    • /
    • 2011
  • 대표적 가스계 소화약제인 $CO_2$의 소화기구(Extinguishing mechanism)를 재조명하기 위하여, $CH_4$/air 저 신장율 대향류 확산화염을 대상으로 $CO_2$가 소화에 미치는 영향에 대한 수치해석이 수행되었다. 부력이 지배적인 화염의 소화현상은 복사 열손실에 의해 큰 영향을 받기 때문에, 소화농도 예측에 대한 복사모델의 성능평가가 우선적으로 이루어졌다. 주요 결과로서, 공기류에 첨가된 $CO_2$의 소화농도는 복사모델이 고려되지 않은 경우 과다 예측되는 반면에, 간략화된 광학적으로 얇은 근사(Optically thin approximation) 모델과 비교적 높은 정확도를 갖는 좁은 밴드(Statistical narrow band) 모델은 실험의 오차범위 내에서 큰 차이를 보이지 않았다. $CO_2$가 소화에 미치는 순수 희석효과, 희석에 의한 복사효과, 화학적 효과 및 열적효과를 정량적으로 분석하기 위하여 가상의 소화약제의 개념을 도입하였다. 이를 통해 화염의 총괄신장율에 따른 $CO_2$ 소화효과에 대한 구체적인 이해가 시도되었다.

  • PDF

디이젤엔진내의 복사열전달 효과에 관한 수치해석적 연구 (Numerical simulations of radiative and convective heat transfer in the cylinder of a diesel engine)

  • 임승욱;김동우;이준식
    • 오토저널
    • /
    • 제14권2호
    • /
    • pp.54-64
    • /
    • 1992
  • During combustion process in a diesel engine radiation heat transfer is the same order of magnitude as the convection heat transfer. An approximation of heat and momentum source distributions is applied at a level consistent with those used in modelling the soot distribution and the turbulence instead of modelling the fuel spray and the chemical kinetics. This paper illustrates a use of the third order spherical harmonics approximation to the radiative transfer equation and delta-Eddington approximation to the scattering phase function for droplets in the flow. Results are obtained numerically by a time marching finite difference scheme. This study aims to compare heat transfer with convection heat transfer and to investigate the importance of scattering by fuel droplets and of accounting for spatial variations in the extinction coefficient on the radiative heat flux distributions at the walls of a disc shaped diesel engine.

  • PDF