• 제목/요약/키워드: 적분방정식

검색결과 667건 처리시간 0.034초

불확실성을 고려한 비선형 수문현상의 해석 (The Analysis of Nonlinear Hydrologic Phenomenon with Uncertainty)

  • 장수형;김상단;윤용남
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2004년도 학술발표회
    • /
    • pp.660-665
    • /
    • 2004
  • 본 연구에서는 비선형적인 메커니즘을 갖는 수문현상의 불확실성을 해석하고자 하는 목적으로 새로운 개념의 지배방정식이 유도된다. 제안된 모형의 불확실성은 토양 특성치의 공간적 변동성에 기인하고 있는 것으로 가정하여, 유도된 방정식은 Fokker-Planckl 방정식의 형태를 가지고 있다. 실제 유역단위에서 토양 내 수분 흐름의 연직방향 흐름을 모의하기 위해 미소단위에서 유도된 Richards 방정식은 토양의 공간적 변동성으로 말미암아 불확실한 매개변수를 갖는 비선형 추계학적 편미분방정식의 형태를 갖게 된다. 이는 먼지 수직 방향적분을 통하여 단순화된 비선형 추계학적 상미분방정식으로 전환되고, 이렇게 전환된 비선형 추계학적 상미분방정식은 다시 추계학적 Liouville 방정식을 이용하여 선형 추계학적 편미분방정식으로 전환되어진다. 최종적으로 cumulant 급수방법을 이용하여 상기 방정식을 선형 결정론적 편미분방정식으로 전환시킴으로써, 강우 시 토양 내 수분 침투현상을 모형화할 경우 유역단위에서 토양의 공간적 변동성을 설명할 수 있는 지배방정식을 유도할 수 있다.

  • PDF

비등방 다공성 2차원 천수모형을 적용한 도시홍수 모델링 (Urban Flood Modeling with Anisotropic Porous Shallow Water Model)

  • 김병현;김현일;한건연
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2020년도 학술발표회
    • /
    • pp.414-414
    • /
    • 2020
  • 고전 천수방정식을 적용한 2차원 도시홍수 모델링에서는 지형의 정확한 반영을 위해 고해상도 격자가 요구되며 이는 많은 계산시간과 노력을 필요로 한다. 최근에는 다공성 천수방정식을 적용한 도시홍수해석으로 많은 계산 노력이 요구되는 도시홍수모델링의 한계를 극복한 연구가 많이 이루어지지고 있다. 이러한 연구는 도시 홍수에서 흐름이 공간적으로 변화할 때 불균일 공극이 존재하므로 격자의 크기를 다르게 하여 이러한 불균일성을 해결하고자 하는 등방성 천수모형의 적용에서 시작되었다. 하지만, 등방성 공극을 고려한 도시홍수 해석모형은 대표요소체적(REV)보다 더 큰 격자의 적용을 해야 하는 제한성을 가진다. 반면, 비등방성 공극은 대표요소체적의 적용이 필요하지 않아 불균일 공극의 크기에 관계없이 이론상으로는 동일한 해상도의 격자가 사용가능하긴 하지만, 실제 도시홍수 해석에서 중요하면서도 도전적인 연구이다. 본 연구에서는 도시홍수의 효율적 계산을 위해 비등방성 공극을 고려한 적분형 다공성 천수방정식을 기반으로 하는 2차원 도시홍수 해석모형을 개발하였다. 모형의 개발을 위해, 적용 격자내에서 도시지역의 건물이 차지하는 길이 및 면적을 산정하고 그 값을 2차원 천수방정식에 적용 가능하도록 체적공극(𝜙j)와 면적공극(𝜓k)을 2차원 고전 천수방정식에 추가하였다. 개발모형은 고전 천수방정식, 등방성 공극 고려(미분형 다공성) 천수방정식 및 비등방성 공극 고려(적분형 다공성) 천수방정식의 적용이 가능하여, 각 모형에 적합한 2차원 격자 생성, 각 모형의 매개변수를 보정 그리고 정확성, 효율성, 적용성이 비교 가능하다. 각 모형의 정확성과 효율성 비교를 위해 3가지의 오차 비교 (구조적 오차, 격자크기 오차, 공극 모형 오차), 계산시간 비교, 공간 변동성 검증을 위한 수심 종단형상 비교하였다.

  • PDF

논문 : 헬리콥터 비행 시뮬레이션을 위한 로터운동방정식 유도 (Papers : Implicit Formulation of Rotor Aeromechanic Equations for Helicopter Flight Simulation)

  • 김창주
    • 한국항공우주학회지
    • /
    • 제30권3호
    • /
    • pp.8-16
    • /
    • 2002
  • 헬리콥터 비행 시뮬레이션을 위한 로터 운동방정식을 implicit formulation 형태로 유도하였다. 좌표계 사이의 상대운동을 고려한 일반화된 벡터 kinematics 를 유도하고 이를 적용하여 브레이드 임의 위치 에서 관성속도 및 관성가속도를 구하였다. 유도된 속도 및 가속도 벡터를 이용하여 플래핑, 리드래그 및 토오크 방정식 등을 implicit form으로 유도하였다. 브레이드 스팬에 따른 공간 적분 방법을 살펴보고, 다양한 힌지형상 및 힌지배열 순서에 관계없이 응용영역을 확장할 수 있음을 밝혔다. DAE(Differential Algebraic Equation) 형태를 갖는 본 연구의 결과식을 이용하여 동특성 계산을 위한 시간적분법을 검토하였다.

안전 구조물의 퍼포먼스 측정시 나타나는 Euler 방정식의 수치해석적 안정성 (Numerical Experiments on the Stability of Euler Equations of the Performance Test of Safety Structures)

  • 고만기;우광성
    • 전산구조공학
    • /
    • 제7권3호
    • /
    • pp.83-93
    • /
    • 1994
  • 도로 안전 구조물을 설계하고 동적 퍼포먼스를 측정하기 위하여 시행되는 충돌시험시 충돌차량의 운동량, 특히 각가속도를 Euler방정식의 수치적분을 통하여 구하는 기법을 설명하였다. 수치적으로 가장 안정적인 9-array 시스템안에 내지된 여러 형태의 7,8-array 서브시스템들의 시스템 미분방정식 및 이들 방정식의 수치 적분시 안정성을 실험데이터를 이용하여 파악하였다. 기본적인 9-array안에 있는 6개의 8-array시스템들은 모두 수치적으로 안정성을 보였고 12개의 7-array시스템들중 6개는 안정성을 보이고 6개는 불안정하였다. 안정성을 보이는 내재된 서브시스템을 활용하면 기본 9-array시스템을 구성하는 센서의 일부 고장시 효과적으로 각가속도를 측정할 수 있는 fail-safe 시스템을 구성할 수 있다.

  • PDF

유전적분형 물성방정식에 근거한 선형 점탄성문제의 시간영역 유한요소해석 (Time-domain Finite Element Formulation for Linear Viscoelastic Analysis Based on a Hereditary Type Constitutive Law)

  • 심우진;이호섭
    • 대한기계학회논문집
    • /
    • 제16권8호
    • /
    • pp.1429-1437
    • /
    • 1992
  • 본 연구에서는 이완형 물성방정식을 바탕으로 하며 프와송 비가 일정하다는 가정을 하지 않는다. 또한 점탄성 지배방정식에 변분원리를 적용하고 유도되어진 식 에 유한요소해법을 사용하여 시스템 기본해석을 위한 연립방정식을 유도한다. 이와 함께 점탄성 물성함수의 유도 및 응력계산을 위한 공식화 과정도 설명한다. 제시된 방법론의 타당성 및 정확성을 보이기 위해서 평면응력 및 평면변형 문제의 변위 및 응력을 수치해석하여 이론해와 비교 검토하며, 아울러 시간증분의 변화와 Gauss poi- nts수가 수치정확도에 끼치는 영향을 조사한다.

저주파 필터를 이용한 수치등각사상에 관한 연구 (A Study on Numerical Conformal Mapping by Low Frequency pass Filter)

  • 송은지
    • 한국정보통신학회:학술대회논문집
    • /
    • 한국해양정보통신학회 2009년도 추계학술대회
    • /
    • pp.821-824
    • /
    • 2009
  • 등각사상은 함수론의 기본적인 문제의 하나로 물리학, 공학 등에 응용되고 있다. 예를 들어 선박이나 비행기의 모양을 고려할 때 나타나는 유체의 문제등에 유용하게 쓰이고 있다. 본 논문에서는 단위원 내부로부터 Jordan 영역 내부에로의 수치등각사상을 구하는 문제를 다루는데 이것은 비선형 적분방정식인 Theodorsen방정식을 푸는 것으로 귀착된다. Theodorsen방정식을 구하는 해법 중 Wegmann의 방법은 계산량에 있어서 가장 유효하다고 알려져 있다. 그러나 수치실험을 통하여 난이도가 높은 문제에 있어 수렴하지 않는 문제점을 발견하게 되었다. 본 논문에서는 수렴하지 않는 불안정의 원인을 이론적으로 분석하고 저주파필터에 의하여 개선한 방법을 제안한다. 개선한 방법에 의한 수치 실험결과를 보고한다.

  • PDF

준설토 유동해석을 위한 유한요소 수식화 (Numerical Formulation for Flow Analysis of Dredged Soil)

  • 신호성
    • 한국지반환경공학회 논문집
    • /
    • 제15권3호
    • /
    • pp.41-48
    • /
    • 2014
  • 준설토에 대한 연구는 주로 준설토의 1차원 침강 및 자중압밀 특성을 파악하는 실험적 연구가 진행되었다. 하지만 양질의 준설지반 확보를 위한 효과적인 투기장의 설계와 배출수에 의한 환경오염을 최소화하기 위해서는 준설토의 투기에 의한 유동특성의 체계적인 연구가 필요하다. 본 연구에서는 준설토 투기장의 펌핑에 의한 토사의 유동 형상을 모사하기 위하여 준설토사를 단일상으로 가정하고 연속 방정식을 유도하여 좌표축에 따른 힘 평형 방정식을 유도하였다. 준설토장의 3차원 거동 해석을 위한 컴퓨터 연산 부하와 모델링 소요시간을 최적화하기 위하여, 토체의 깊이 방향으로 적분을 수행하는 깊이 적분 방법을 지배 방정식에 적용하여, 3차원적 지형조건을 고려할 수 있도록 하였다. 지배 방정식의 보간함수를 이용한 공간분할에서 Petrov-Galerkin 수식화 기법을 적용하였다. 일반화된 사다리꼴 법칙으로 시간적분을 수행하고 Newton의 반복과정을 이용할 수 있도록 근사화시켰다. 가중행렬은 DG과 CDG 기법을 적용하였으며, 준설토 유동해석에서 가중행렬에 따른 수치적인 안정성을 평가하기 위하여 사각형 기둥 슬럼프 시뮬레이션을 수행하였다. 수치기법에 대한 비교 분석 결과는 DG 기법을 적용한 SU/PG 수식화가 유사진동을 최소화시키는 가장 안정적인 수치해석결과를 보여주는 것으로 나타났다.

스펙트럴 요소법을 이용한 SH파 전파의 외연적 시간이력해석 (Explicit Transient Simulation of SH-waves Using a Spectral Element Method)

  • 윤승욱;강준원
    • 한국전산구조공학회논문집
    • /
    • 제31권2호
    • /
    • pp.87-95
    • /
    • 2018
  • 이 논문에서는 스펙트럴 요소법과 외연적 시간적분법을 이용해 SH파의 전파 거동을 계산하는 수치해석 기법을 제시한다. 2차원 영역에서의 탄성파 해석을 위해 해석영역을 유한 영역으로 한정하고 파동이 반사되지 않도록 수치적 파동흡수 경계조건인 perfectly matched layer(PML)를 도입하였다. PML이 포함된 시간영역 파동방정식의 유한요소해법을 위해 스펙트럴 요소법을 적용하였고 Legendre- Gauss-Lobatto 수치적분법을 사용하여 질량행렬을 대각화하였다. 2차 미분방정식 시스템의 파동방정식을 1차 미분방정식 시스템으로 변환하였고 병렬화를 통한 탄성파 해석 성능의 최적화를 위해 외연적 시간적분법인 4차 Runge-Kutta 방법을 이용해 해석영역에서의 변위응답을 계산하였다. 2차원 해석영역에서 SH파의 전파 거동을 계산하는 수치예제를 통해 제시한 외연적 스펙트럴 요소법의 정확성을 검증하였고 PML로 인한 반사파의 감쇠효과를 확인하였다. 외연적 시간적분법을 통한 탄성파 해석 기법은 3차원 영역과 같은 대규모 문제에서의 탄성파 수치해석을 효율적으로 수행할 수 있을 것으로 기대된다.