• 제목/요약/키워드: linear approximation

검색결과 748건 처리시간 0.02초

삼중대각행렬 시스템 풀이의 빠른 GPU 구현 (Fast GPU Implementation for the Solution of Tridiagonal Matrix Systems)

  • 김영희;이성기
    • 한국정보과학회논문지:시스템및이론
    • /
    • 제32권11_12호
    • /
    • pp.692-704
    • /
    • 2005
  • 컴퓨터 하드웨어의 급속한 발전으로 그래픽 프로세서 유닛(Graphics Processor Units : GPUs)은 굉장한 메모리 대역폭과 산술 능역을 보유하게 되어 범용 계산에 많이 활용되고 있으며, 특히 계산 집약적인 물리 기반 시뮬레이션(physics based simulation)의 GPU 구현이 활발하게 연구되고 있다. 물리 기반 시뮬레이션의 기본이 되는 미분방정식 풀이 과정에서 삼중대각행렬(tridiagonal matrix) 시스템은 유한차분(finite-difference) 근사에 의해서 자주 나타나는 선형시스템으로 물리 기반 시뮬레이션 관점에서 삼중대각행렬 시스템의 빠른 풀이는 중요한 연구 분야이다. 본 논문에서는 GPU에서 삼중대각행렬 시스템 풀이를 빠르게 구현할 수 있는 방법을 제안한다. 벡터 프로세서(vector processor) 계산에서 삼중대각행렬 시스템 풀이 방법으로 널리 사용되는 cyclic reduction 또는 odd-even reduction 알고리즘을 GPU에서 구현하였다. 본 논문에서 제안한 방법을 삼중대각행렬 시스템 풀이 방법으로 잘 알려져 있는 Thomas 방법과 GPU를 이용한 선형시스템 풀이에서 좋은 성과를 보이고 있는 conjugate gradient 방법과 비교할 때 상당한 성능 향상을 얻을 수 있었다. 또한, 열전도(heat conduction) 방정식, 이류 확산(advection-diffusion) 방정식, 얕은 물(shallow water) 방정식에 의한 물리 기반 시뮬레이션의 GPU 구현에 본 논문에서 제안한 방법을 사용하여 1024x1024 격자의 계산 영역에서 초당 35프레임 이상의 놀라운 성능을 보여주었다.

Direct numerical simulations of viscoelastic turbulent channel flows at high drag reduction

  • Housiadas Kostas D.;Beris Antony N.
    • Korea-Australia Rheology Journal
    • /
    • 제17권3호
    • /
    • pp.131-140
    • /
    • 2005
  • In this work we show the results of our most recent Direct Numerical Simulations (DNS) of turbulent viscoelastic channel flow using spectral spatial approximations and a stabilizing artificial diffusion in the viscoelastic constitutive model. The Finite-Elasticity Non-Linear Elastic Dumbbell model with the Peterlin approximation (FENE-P) is used to represent the effect of polymer molecules in solution, The corresponding rheological parameters are chosen so that to get closer to the conditions corresponding to maximum drag reduction: A high extensibility parameter (60) and a moderate solvent viscosity ratio (0.8) are used with two different friction Weissenberg numbers (50 and 100). We then first find that the corresponding achieved drag reduction, in the range of friction Reynolds numbers used in this work (180-590), is insensitive to the Reynolds number (in accordance to previous work). The obtained drag reduction is at the level of $49\%\;and\;63\%$, for the friction Weissenberg numbers 50 and 100, respectively. The largest value is substantially higher than any of our previous simulations, performed at more moderate levels of viscoelasticity (i.e. higher viscosity ratio and smaller extensibility parameter values). Therefore, the maximum extensional viscosity exhibited by the modeled system and the friction Weissenberg number can still be considered as the dominant factors determining the levels of drag reduction. These can reach high values, even for of dilute polymer solution (the system modeled by the FENE-P model), provided the flow viscoelasticity is high, corresponding to a high polymer molecular weight (which translates to a high extensibility parameter) and a high friction Weissenberg number. Based on that and the changes observed in the turbulent structure and in the most prevalent statistics, as presented in this work, we can still rationalize for an increasing extensional resistance-based drag reduction mechanism as the most prevalent mechanism for drag reduction, the same one evidenced in our previous work: As the polymer elasticity increases, so does the resistance offered to extensional deformation. That, in turn, changes the structure of the most energy-containing turbulent eddies (they become wider, more well correlated, and weaker in intensity) so that they become less efficient in transferring momentum, thus leading to drag reduction. Such a continuum, rheology-based, mechanism has first been proposed in the early 70s independently by Metzner and Lamley and is to be contrasted against any molecularly based explanations.

200 Gbps WDM 시스템에서 Mid-Span Spectral Inversion 기법을 이용한 채널 왜곡의 보상 특성 (Compensation Characteristics of Distorted Channels in 200 Gbps WDM Systems using Mid-Span Spectral Inversion Method)

  • 이성렬
    • 한국전자파학회논문지
    • /
    • 제14권8호
    • /
    • pp.845-854
    • /
    • 2003
  • 광섬유의 색 분산과 비선형 효과에 의해 왜곡된 채널 신호를 보상하기 위한 방법으로 고비선형 분산 천이 광섬유의 광 위상 공액기를 이용한 경로 평균 강도 근사 MSSI(Mid-Span Spectral Inversion)를 채택한 채널 전송률 40 Gbps의 5-채널 WDM시스템의 시뮬레이션을 통하여 채널 각각의 MSSI 보상 특성을 다양한 광섬유 분산 계수에 따라 분석해 보았다. 채널별 보상 특성의 분석은 1 dB 기준 눈 열림 패널티, 수신단에서의 비트 에러율 특성, $10^{-9}$ BER에서의 채널 간 파워 패널티 등을 이용했다. 시스템 전송 길이, 광섬유의 분산 계수, 광 위상 공액기 펌프 광 파장, 광 위상 공액기에서의 WDM신호 파장에 따른 공액파 변환 효율 등에 관계하여 광 위상 공액기를 중심으로 한 첫 번째 전송 링크에서의 신호 평균 전력과 두 번째 전송 링크에서의 공액파 평균 전력을 동일하게 만드는 펌프 광 전력을 적절히 선택하면 고비선형 분산 천이 광섬유가 MSSI 보상을 통한 광대역 WDM 시스템을 위한 광 위상 공액기에서의 비선형 매질로 매우 적합하다는 것을 확인할 수 있었다.

Skeletal stability following mandibular advancement: is it influenced by the magnitude of advancement or changes of the mandibular plane angle?

  • Tabrizi, Reza;Nili, Mahsa;Aliabadi, Ehsan;Pourdanesh, Fereydoun
    • Journal of the Korean Association of Oral and Maxillofacial Surgeons
    • /
    • 제43권3호
    • /
    • pp.152-159
    • /
    • 2017
  • Objectives: The aim of this study was to investigate the effects of advancement magnitude and changes in mandibular plane angle on the stability of mandibular advancement. Materials and Methods: This retrospective cohort study evaluated the postoperative stability of mandibular advancement in class II skeletal subjects who underwent bilateral sagittal split osteotomy. Radiographs taken preoperatively, immediately postoperatively and 1 year postoperatively were traced and analyzed using linear and angular measurements. To determine horizontal and vertical relapse, an X-Y coordinate system was established in which the X-axis was constructed by rotating S-N downward by $7^{\circ}$ (approximation of the Frankfort horizontal plane) and the Y-axis was defined as a line perpendicular to the X-axis and passing through the point Sella. For certain reference points including point A, point B, pogonion and menton, the perpendicular distance between each point and both axes was determined and cephalometric variables were recorded as X and Y coordinates. Results: Twenty-five subjects were studied. A significant correlation between the amount of mandibular advancement and relapse in the B point (vertical and horizontal) and the pogonion point was observed (vertical and horizontal, P<0.001). Evaluation of data demonstrated a positive correlation between the mandibular plane angle (SN/ML) change and vertical relapse in the B point (P<0.05). A simple regression model demonstrated that 74% of horizontal relapse and 42.3% of vertical relapse in the B point was related to the amount of mandibular advancement. The receiver operating characteristic test showed that 8.5 mm mandibular advancement is related to a relapse rate of 1 mm or more in the pogonion, vertically or horizontally. Conclusion: The magnitude of mandibular advancement is a stronger surgical predictor for horizontal rather than vertical relapse at the B point. Changes in mandibular plane angle (SN/ML) during surgery affect vertical, but not horizontal relapse at the B point.

산화물 반도체식 가스센서의 입출력 고차 캘리브레이션 방정식 도출 (Induction of the High Order Calibration Equation of Metal Oxide Semiconductor Gas Sensors)

  • 박규태;김강민;이형기;윤명섭
    • 한국가스학회지
    • /
    • 제24권2호
    • /
    • pp.44-49
    • /
    • 2020
  • 본 논문에서는 산화주석 기반 반도체식 가스센서의 제조사양을 분석하여 측정회로를 설계하였다. 또한, 실제 가스안전관리 산업현장에서 가스센서를 사용하기 위하여 응답특성시험을 통한 최적 입출력 방정식을 도출하였다. 산화물 반도체 가스센서의 응답특성은 6종의 농도를 가진 표준가스를 제조하여 주입하며 출력전압을 측정을 통해 분석되었다. 각각 농도가 다른 가스의 주입에 대한 센서의 출력 전압의 관계를 룩업 테이블로 생성하였다. 생성된 룩업테이블의 데이터가 등간격이면 근사함수법으로 다항식을 도출할 수 있다. 그래서 5차의 다항식을 정의하고, 그 계수를 최소자승법으로 구하여 센서의 입출력 특성을 대표하는 5차 다항식을 완성하였다. 제안된 5차 다항식이 가스검지기 등에 실제로 적용되기 위해서는 다항식의 역변환이나 사이 값이 추가된 배열형태로 프로그래밍이 다시 생성되어져야 한다. 본 연구에서는 사이 값을 1/5로 세밀하게 구분 및 보간을 계산하여 배열형태로 구현되었다. 제안한 5차 캘리브레이션 교정식의성능은 표준가스를 주입 및 그 농도를 측정했을 때 일차방정식보다 오차가 감소됨을 확인되었다.

밀도 변화를 포함하는 3차원 연직함수 전개모형 (A Three-Dimensional Galerkin-FEM Model with Density Variation)

  • 이호진;정경태;소재귀;강관수;정종율
    • 한국해안해양공학회지
    • /
    • 제8권2호
    • /
    • pp.123-136
    • /
    • 1996
  • 본 연구에서는 시 공간적으로 변화하는 밀도구조를 고려한 3차원 연직함수 전개모형을 제안한다. 정수압을 가정하며 열보존식과 상태방정식의 도입으로 밀도변화가 고려된다 수평방향으로의 변화는 유한 차분 격자상에서 계산되며 유속 및 수온의 연직구조는 선형보간함수를 사용하여 계산된다. 구성된 행렬방정식은 유사변환 개념을 도입하여 시간적분된다. 개발된 모델의 테스트를 위해 간단한 이상해역과 황해역에서 대기와 해양간의 열교환에 따른 수온구조 변동을 실험하였다. 이류효과는 열수송방정식에서만 고려하였으며 연직 와동점성계수와 와동확산계수는 시 공간적으로 일정한 값을 사용하였다. 이상해역에서의 수치실험결과 모델영역의 수심의 차이에 따른 열저장의 차이로 인해 수온의 수평적 구배가 발생하였다. 결과적으로 전향력을 고려하지 않을 경우에는 상층에서는 수심이 증가하는 방향으로 흐름이 발생하고 하층에는 반대방향의 흐름이 유도된 반면 전향력을 고려할 경우에는 수온차에 의한 압력구배력과 전향력이 균형을 이루면서 지형류가 뚜렷하게 나타났다. 황해역에서는 복잡한 흐름이 나타났지만 전체적으로는 지형류의 특성이 우세하게 나타났다.

  • PDF

차분 간섭도로부터 지표변위의 시계열 관측을 위한 개선된 Small Baseline Subset (SBAS) 알고리즘 (Improvement of Small Baseline Subset (SBAS) Algorithm for Measuring Time-series Surface Deformations from Differential SAR Interferograms)

  • 정형섭;이창욱;박정원;김기동;원중선
    • 대한원격탐사학회지
    • /
    • 제24권2호
    • /
    • pp.165-177
    • /
    • 2008
  • 최근 spatial decorrelation을 극복하기 위하여 비교적 짧은 기선(baseline)을 지니는 여러 시기의 차분간섭도(differential interferogram)로부터 시계열 지표변위(time-series surface deformation)를 관측할 수 있는 small baseline subset(SBAS) 알고리즘이 개발되었다. 이 SBAS 알고리즘은 singluar value decomposition(SVD)을 이용하여 시간별로 완벽하게 연결되지 않는 차분간섭도로부터 시계열 지표변위를 관측하였을 뿐 아니라 공간적으로 저주파 필터와 시간적으로 고주파 필터를 이용하여 대기효과를 보정하였다. 그러나 이 알고리즘은 초기 관측시 시계열 지표변위를 선형으로 가정하였기 때문에 각 차분간섭도의 phase unwrapping 오차를 정확하게 보정하기 어려웠을 뿐 아니라 시계열의 지표변위에 존재하는 노이즈 성분을 완화시키지 못했다. 이와 같은 단점을 보완하기 위하여 이 연구에서는 기존의 SBAS 알고리즘을 개선하였다. 이 개선된 SBAS 알고리즘은 각 차분간섭도의 phase unwrapping 오차를 최소화하기 위하여 반복적으로 시계열 지표면 변위를 개선하였고, 시계열 지표변위의 노이즈를 제거하기 위하여 유한차분근사법(finite difference approximation)을 이용하였다. 서로 다른 지역의 26개의 ERS-12자료와 21개의 RADARSAT-1 fine beam (F5) 자료를 이용하여 개선된 SBAS 알고리즘을 실험하고 분석하였다. ERS-1/2자료에서는 LOS(line-of-sight) 지표변위가 약 13년 동안 최대 -40cm가 관측되었고, RADARSAT-1 fine beam 자료에서는 약 2년 동안 최대 -3cm의 LOS 지표변위가 관측되었다.

응력근사해법(應力近似解法)을 이용한 평면(平面)트러스구조물(構造物)의 형상최적화(形狀最適化)에 관한 연구(研究) (Optimization of the Truss Structures Using Member Stress Approximate method)

  • 이규원;유희중
    • 대한토목학회논문집
    • /
    • 제13권2호
    • /
    • pp.73-84
    • /
    • 1993
  • 본(本) 연구(研究)에서는 분할기법(分割技法)을 이용하여 평면(平面)트러스구조물(構造物)의 형상최적화(形狀最適化)를 시도(試圖)하였다. 본(本) 연구(研究)의 제(第)1단계(段階)(Level 1)에서는 다른 연구(研究)와 달리 응력제약(應力制約)을 감도해석(感度解析)에 효율적(效率的)이라고 알려진 설계공간법(設計空間法)에 의해서 부재응력근사화(部材應力近似化)를 하므로서 비선형최적화문제(非線形最適化問題)가 선형계획문제(線形計劃問題)로 변환(變換)되어 해(解)를 효율적(效率的)으로 구할 수 있고 또한 감도해석(感度解析)을 위한 구조해석수(構造解析數)를 줄일 수 있다. 목적함수(目的凾數)는 구조물(構造物)의 중량(重量)이 최소(最小)가 되도록 중량함수(重量凾數)를 택하였다. 제약조건식(制約條件式)으로는 허용응력(許容應力), 좌굴응력(挫屈應力), 변위제약(變位制約) 및 설계변수(設計變數) 상하한치제약(上下限値制約)을 부과(附課)하였고 다(多) 재하조건(載荷條件)을 고려(考慮)하여 최적화문제(最適化問題)를 형성(形成)하였다. 제(第)2단계(段階)(Level 2)에서는 설계변수(設計變數) 및 조정변수(調整變數)를 절점좌표(節點座標)로 하고 목적함수(目的凾數)로는 중량함수(重量凾數)로 하여 최적화문제(最適化問題)를 형성(形成)하였다. 절점좌표(節點座標)만을 설계변수(設計變數)로 하므로서 무제약최적화문제(無制約最適化問題)로 형성(形成)되므로 최적화(最適化) 과정(過程)이 용이(容易)하다. 본(本) 연구(研究)의 제(第)1단계(段階)에서는 부재응력(部材應力)을 근사화(近似化)하여 단면(斷面)을 최적화(最適化)하고 제(第)2단계(段階)에서는 형상(形狀)만 최적화(最適化)하는 분할기법(分割技法)을 트러스구조물(構造物)에 적용(適用)한 결과 본(本) 연구(研究)는 트러스구조물(構造物)의 형태(形態), 제약조건식(制約條件式)에 구애받지 않고 최적해(最適解)에 부재응력근사화(部材應力近似化)로 인하여 효율적(效率的)으로 수렴(收斂)하였고 또한 타(他)의 연구(研究)와 거의 동일(同一)한 연구(研究) 결과(結果)를 얻었으며 형상최적화(形狀最適化)로 트러스구조물(構造物)의 중량(重量)을 5.4% - 15.4% 까지 감소(減少)시켰다.

  • PDF