• 제목/요약/키워드: 시간 영역 유한 차분 방정식

검색결과 54건 처리시간 0.023초

질량 원천항을 이용한 내부조파 (Internal Wave Generation with the Mass Source Function)

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

  • PDF

광대역 전자파를 이용한 역산란 해석 연구 (Analysis of Microwave Inverse Scattering Using the Broadband Electromagnetic Waves)

  • 이정훈;정용식;소준호;김준연;장원
    • 한국전자파학회논문지
    • /
    • 제17권2호
    • /
    • pp.158-164
    • /
    • 2006
  • 본 논문에서는 시간 영역 유한차분법(FDTD: Finite-Difference Time-Domain Method)과 설계 민감도법(design sensitivity analysis)을 이용하여 미지의 유전체 산란체(dielectric scatterer) 복원을 위한 역산란 문제(inverse scattering)의 새로운 해석 기법을 제안하였다. 본 연구에서는 빠른 수렴을 위하여 목적 함수의 도함수를 이용한 설계 민감도법을 도입하였고, 시간 영역 유한차분법으로부터 직접 설계 민감도 수식을 도출하였다. 해석의 효율성을 위하여 보조 변수법(adjoint variable method)을 도입하여 보조 변수 방정식을 도출하고 최적화 알고리듬으로 최대 경사도법을 이용하여 반복적인 추정을 통하여 미지의 유전체를 복원하였다. 본 연구의 타당성의 보이기 위하여 2차원 $TM^z$에서의 유전체 복원 사례를 제시한다.

시추공간 레이다 측정에서 전자기 도파관 효과의 수치모델링 (Numerical modelling of electromagnetic waveguide effects on crosshole radar measurements)

  • 장한누리;박미경;김희준
    • 지구물리와물리탐사
    • /
    • 제10권1호
    • /
    • pp.69-76
    • /
    • 2007
  • 시추공 지중레이다 고주파수 전자기파 전파와 관련된 현상은 복잡하다. 그 물리적 과정을 보다 잘 이해하기 위해서 본 연구에서는 원통좌표계에서 맥스웰방정식의 유한차분 시간영역 해를 이용하였다. 이 방법은 시추공간 레이다탐사에서 전자기파의 전파형을 모델링할 수 있다. 그리고 원통좌표계는 계산 효율이 높고, 파동장의 3 차원적인 기하학적 확산을 정확하게 계산할 수 있으며 또 작은 크기의 시추공을 이산화할 때 효과적이다. 수치계산 결과를 통해 물로 채워진 시추공은 전파 파형에 영향을 주는 강한 도파관 효과를 일으킬 수 있으며, 그 도파관 효과는 시추공의 크기와 안테나의 길이에 의존한다.

시간영역 유한차분법을 이용한 휴대용 전화기의 모노폴 안테나 특성해석 (Analysis of the monopole antenna characteristcs of handy phone using Finite Difference Time Domain(FDTD) Method)

  • 손영수;윤현보
    • 한국전자파학회지:전자파기술
    • /
    • 제6권3호
    • /
    • pp.3-14
    • /
    • 1995
  • 800MHz대 휴대용 전화기에서, 모노폴 안테나의 입력임피던스, 입력전력 및 복사패턴을 시간영역 유한차분 볍(FDTD)을 이용하여 계산하였다 .. FDTD로 모노폴 안테나의 주파수특성 해석을 위해, 휴대용 전화기를 모노 폴 안테나가 도체 박스에 접속된 형태로 모델링하고, 도션구조의 FDTD 해석을 위해, Yee 알고리즘을 가는 도 선 근사법 및 백스웰 적분방정식을 이용하여 변형한 알고리즘[U]을 적용하였다. 또한, 도체 박스에 의한 모노폴 안테나의 복사특성 영향을 관찰하기위해, 모노폴 안테나 및 도체 박스주위에 분포하는 주변자계로 부터, 전 류분포를 직접 구하여 복사패턴을 계산하는 방볍을 제시하고, 기존의 FDTD 복사패턴 계산방법에 의한 결과와 비교하였다. 모든 FDTD 계산결과는 휴대용 전화기의 전체 길이를 $\lambda$/2가 되도록 제작, 측정결과와 비교, 잘 열 치하고 있음을 확인하였다.

  • PDF

충격하중을 받는 구조부재의 탄소성 파괴해석 프로그램 개발 (Development of Elastic-Plastic Fracture Analysis Program for Structural Elements under an Impact Loadings)

  • 김경수;박준범
    • 대한조선학회논문집
    • /
    • 제35권1호
    • /
    • pp.61-71
    • /
    • 1998
  • 본 논문에서는 2차원 공간에서의 탄성 또는 탄소성응력파를 받는 구조부재의 동적 파괴거동을 다룬다. 이러한 문제에 대한 지배방정식은 운동방정식과 탄소성 구성방정식에 대한 증감식으로 구성된 쌍곡선 편미분 방정식으로 나타나고, 이를 풀기 위해 유한차분법을 기초로 한 Zwas방법이 도입된다. 또한 탄소성문제의 동적거동을 나타내기 위해 응력공간내 탄소성 loading path가 소성항복 현상을 모델링하는데 제안된다. 이러한 계산결과를 바탕으로 탄성체의 균열선단의 동적응력확대계수가 계산되어지고, 탄소성체에 대한 소성영역의 형상의 시간이력을 보여준다.

  • PDF

코앤다 효과를 이용한 순환 제어 날개의 수치적 연구 (A Numerical Study on a Circulation Control Foil using Coanda Effect)

  • 박제준;이승희
    • 대한조선학회논문집
    • /
    • 제37권2호
    • /
    • pp.70-76
    • /
    • 2000
  • 2차원 순환 제어 날개 주위의 점성 유동을 수치적으로 해석하여 선박 해양 분야에 적용하기 위한 기초적 연구를 수행하였다. 지배 방정식으로는 2차원 비 압축성 RANS 방정식과 연속 방정식을 사용하였으며, 유한 차분법과 MAC법을 이용하여 이를 차분화 하였다. 또한, 압력 Poisson 방정식은 SOR 법을 이용하여 계산하였고, 모든 물리량이 격자점에 위치하도록 하였으며, 슬롯의 형상을 잘 표현하기 위해 격자계가 엇갈린 O형 비교차 중첩 격자계를 도입하였다. 코앤다 효과는 주로 난류 영역에서 나타나므로 수정된 Baldwin-Lomax 난류 모델을 도입하였다. 본 연구에서는 이러한 수치 계산 조직을 활용하여 뒷날이 원형으로 변형된 20% 두께의 타원 날개에 대한 코앤다 효과를 조사하기 위하여 분사 모멘텀을 달리 하며 시간에 따른 양 항력 특성의 변화를 계산하고, 그 결과를 실험치와 비교하여 잘 일치함을 보였다.

  • PDF

시간영역에서 가우스뉴튼법을 이용한 탄성파 파형역산 (Time Domain Seismic Waveform Inversion based on Gauss Newton method)

  • 신동훈;박창업
    • 한국지구물리탐사학회:학술대회논문집
    • /
    • 한국지구물리탐사학회 2006년도 공동학술대회 논문집
    • /
    • pp.131-135
    • /
    • 2006
  • 본 논문에서는 가우스 뉴튼법을 이용한 중합전 탄성파 자료의 파형역산에 관한 연구를 수행하였다. 탄성파 파형역산에 가우스 뉴튼법을 적용하는 방법은 80년대에 제시되었으나 최근 들어서야 활발히 연구가 진행되고 있는데 이는 연산 능력과 기억용량의 한계에 기인한 것이다. 이를 극복하기 위해 본 연구에서는, 파동 전파 수치모의와 역산과정에서 각각 다른 크기의 격자간격을 사용하고, 필요한 시간영역의 파동전파 모사와 가상 진원의 근사를 통해 편미분 파형을 계산하였으며, 효과적으로 슈퍼컴퓨터를 활용하기 위해 병렬처리 기법을 사용하였다. 수치모의를 통해, 가우스 뉴튼법을 이용한 파형 역산의 수렴속도가 빠르고 정확한 것을 알 수 있었으며, 이를 통해 본 연구에서 제시한 방법의 실제 탄성파 자료를 이용한 역산에의 적용가능성을 확인하였다.

  • PDF

반도체 공정을 위한 OpenMP와 MPI 기반의 FDTD 시뮬레이션 연산 환경 구축 (Implementation of FDTD simulation using OpenMP and MPI for semiconductor process)

  • 이승일;이철훈
    • 한국콘텐츠학회:학술대회논문집
    • /
    • 한국콘텐츠학회 2015년도 춘계 종합학술대회 논문집
    • /
    • pp.59-60
    • /
    • 2015
  • 반도체 공정에서는 소자 내부의 물리량 계산을 통해 결점를 검출하는 시뮬레이션을 수행하게 된다. 이를 위해 3차원적 형상을 표현하여 물리적 미분 방정식을 계산하는 유한 차분 시간 영역 알고리즘(Finite-Difference Time-Domain, 이하 FDTD)과 같은 수치해석 기법이 사용된다. 반도체 테스트 범위의 확장으로 인해 시뮬레이션 사이즈 또한 커지고 있는 추세이다. 이에 따라 하나의 프로세서에서 수행할 수 없는 문제가 발생하기도 한다. 이를 해결하기 위해 본 논문에서는 openMP와 MPI를 이용한 하이브리드 컴퓨팅 기법을 바탕으로 다중 GPU 제어를 통해 시뮬레이션 환경을 구축하였으며 정상 동작함을 확인하였다.

  • PDF

급전선과 직교된 전자기결합 마이크로스트립 다이폴 배열안테나의 해석 (Analysis of the Transversely fed EMC Microstrip Dipole Array Antenna)

  • 손영수;윤현보
    • 한국전자파학회논문지
    • /
    • 제7권2호
    • /
    • pp.105-116
    • /
    • 1996
  • 급전선과 직교된 전자기결합 마이크로스트립 다이폴안테나는 주로 주파수 영역의 적분방정식 및 모멘트법에 의해, 다이폴의 전류분포를 구하므로써 설계 및 해석이 이루어졌으나, 본 연구에서는 FDTD 방법을 이용하여 각 다이폴의 전류분포를 직접 계산하여, 전자기결합 다이폴 배열안테나의 설계 및 해석 가능성을 제시하였다. 이 경우, 맥스웰방정식 적분형을 유한차분형태로 표현한 수식을 적용하였으며, 계산된 각 다이폴의 전류분포에 의해 복사전계를 구하므로써, 표면전류 및 자류밀도를 계산하여 복사패턴을 구하는 방법에 비해 계산과정과 시간이 크게 단축됨을 보였다.

  • PDF

내부조파를 활용한 파랑 변형 수치모의 (Numerical simulation of wave transformation with the internal wave maker)

  • 하태민;조용식
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2011년도 학술발표회
    • /
    • pp.43-43
    • /
    • 2011
  • 최근 컴퓨터 기술의 발전으로 수심방향으로 완전한 운동방정식인 Navier-Stokes 방정식을 이용하는 3차원 수치모형에 의한 연구가 기초과학 분야인 수학에서부터 공학은 물론 의학 분야까지 광범위하게 진행되고 있다. 3차원 수치모형을 이용한 연구는 이론적으로 매우 우수하긴 하나 정확도 높은 결과를 얻기 위해서는 매우 조밀한 격자를 필요로 하기 때문에 아직까지 막대한 계산시간이 필요하다는 단점이 있으나 근래의 컴퓨터 기술이 엄청난 속도로 발전하고 있는 점을 감안할 때 적용 가능성은 계속 높아지고 있다. 이러한 흐름에 따라 해안공학 분야에서도 3차원 수치모형을 이용한 다양한 연구가 시도되고 있다. 복잡한 해안구조물로 인한 파랑변형 및 구조물의 안전성 등을 검토할 때 기존의 평면 2차원 수치모형이나 연직 2차원 수치모형으로는 재현이 힘들어 상대적으로 수치모형을 활용한 연구의 효용성이 낮았던 것에 비해 3차원 수치모형을 활용할 경우 복잡한 지형 및 구조물의 형상을 재현할 수 있기 때문에 좀 더 유용하게 사용될 수 있다. 한편, 파랑변형을 다루는 수치모형실험을 수행할 때 외부조파를 사용할 경우 구조물이나 지형에 의해 반사되어 나온 파랑이 조파지점에 도달할 때 실험영역으로 재 반사되는 문제가 발생한다. 이를 해결하기 위해 내부조파기법의 개발에 대한 연구가 필수적이었으며, 2차원 수치모형을 중심으로 그 연구가 매우 활발하게 진행되어 왔으나 Navier-Stokes 방정식 모형의 경우 상대적으로 연구가 미흡하였다. 본 연구에서는 3차원 Navier-Stokes 방정식 모형에서 내부조파기법을 도입하여 목표 파랑을 조파하였다. 수치모형은 Navier-Stokes 방정식을 엇갈림 격자체계에서 계산하는 유한차분모형으로 자유수면 추적에는 2차 정확도의 VOF(volume-of-fluid) 기법을 사용한다. 수치모형실험 결과를 일정한 수심을 전파하는 정현파의 해석해와 비교하였으며, 수치모형실험 결과가 비교적 정확하게 목표 파랑을 조파할 수 있음을 검증하였다.

  • PDF