• 제목/요약/키워드: Newton-Raphson

검색결과 588건 처리시간 0.03초

연소 반응을 가지는 후류 유동의 불안정성 (Stability Analysis of Wakes with Chemical Reaction)

  • 신동신;홍성제
    • 한국추진공학회:학술대회논문집
    • /
    • 한국추진공학회 1998년도 제10회 학술강연회논문집
    • /
    • pp.17-17
    • /
    • 1998
  • 연소반응을 가지는 후류(wake)는 가스터빈 연소실의 flame holder 등에서 발생한다. 후류유동의 안정성 혹은 불안정성은 이러한 유동에 있어서 많은 영향을 끼치므로 상당히 중요하다. 본 연구는 위와 같이 연소 반응에 의한 밀도구배를 가지는 후류유동에 대하여 불안정성 해석을 수행하였다. 교란에 대한 지배방정식은 Navier-Stokes 방정식에서 점성항을 제외한 Euler 방정식을 고려하였다. 충류유동의 압력은 일정하다고 가정하였다. 교란 방정식을 유도하기 위하여 충류 유동이 평행하여 유동 방향에 수직한 방향의 구배만이 존재한다고 가정하였다. 모든 변수들은 충류 유동의 값과 움직이는 파장의 형태를 가지는 작은 교란의 합으로 생각하여 압력에 대한 교란방정식을 구하며, shooting법과 Newton-Raphson법에 근거를 두는 반복법을 사용하여 풀었다. 불안정성 해석을 수행하는 기본 유동의 속도장 및 온도장은 불안정성 해석을 수행하는 기본 유동의 속도장 및 온도장은 비압축성의 경우 우선 Gaussian Profile 가정함과 동시에 연소반응을 포함하는 유동장을 보다 정확히 구하기 위하여 Navier-Stokes 방정식으로부터 구한 결과를 사용하였다. 연소반응을 포함하는 유동장을 구할 때에는 해석상 편의를 위해 예혼합물질은 이상기체로, 화학반응은 1단계의 비가역반응으로 가정하였으며, 반응열로 인한 부력의 효과는 무시하였다. 위와 같은 유동장을 가지고 불안정성 해석을 수행한 결과 후류유동은 두 개의 변곡점을 가지며 sinuous 모드와 varicose모드의 두 개의 불안정성 모드가 있음을 보였다. 밀도 변화가 있는 경우나 밀도 변화가 없는 경우 모두 sinuous 모드의 가장 불안정한 모드가 varicose 모드의 가장 불안정한 모드보다 더 불안정함을 보여주어 후류 유동은 자유 유동에 가까운 위상 속도를 가지는 sinuous 모드에 의해 지배될 것임을 예측할 수 있다. 연소반응이 완전연소에 가까울수록 그리고 압축성 효과가 클수록 유동내부의 온도가 증가하고 점성 또한 증가하여 후류유동은 안정됨을 알 수 있었다 유동변수들의 contour로부터 유동의 특성을 예측한 결과 baroclinic 항이 dilatational 항보다 상대적으로 크며, 중심선 상하에 생기는 vortex를 더욱 성장시킬 것으로 생각된다.

  • PDF

선형혼합모형을 이용한 유전체 자료분석방안에 대한 연구 (Efficient strategy for the genetic analysis of related samples with a linear mixed model)

  • 임정민;성주헌;원성호
    • Journal of the Korean Data and Information Science Society
    • /
    • 제25권5호
    • /
    • pp.1025-1038
    • /
    • 2014
  • 가족 자료를 활용한 연속형 표현형의 전장유전체분석 (genome-wide association analysis)은 주로 선형혼합모형을 이용하며, 분산공분산행렬은 가족 구성원간의 유전적 거리를 고려하여 결정된다. 그러나 가족 구성원들의 표현형의 유사성은 유전적 요인과 환경적 요인에 의하여 발생함에도 불구하고, 표현형의 유사성은 단지 유전적 요인에 의해서 발생한다고 가정한다. 예를 들어 키의 경우 부부 사이에 양의 상관관계가 존재하나 유전적 요인만 고려하여 독립으로 가정한다. 선형혼합 모형에서 분산공분산 구조를 잘못 가정하는 경우, 검정통계량의 1종 혹은 2종의 오류를 적절히 관리할 수 없다. 본 논문에서는 다양한 유형의 분산공분산구조를 가정할 수 있는 선형혼합모형과 이를 기반으로 한 검정통계량을 제안하였다. 모의실험을 통하여 제안한 방법이 기존의 모형보다 통계적 검정력이 우수함을 확인하였다. 또한 체질량지수 (body mass index; BMI)의 전장유전체 분석에 적용하여 기존에 알려지지 않은 새로운 원인 유전자를 규명하였다.

Meta-Heuristic Algorithms를 이용한 확률분포의 매개변수 추정 (Parameters Estimation of Probability Distributions Using Meta-Heuristic Algorithms)

  • 윤석민;이태삼;강명국;정창삼
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2012년도 학술발표회
    • /
    • pp.464-464
    • /
    • 2012
  • 수문분야에 있어서 빈도해석의 목적은 특정 재현기간에 대한 발생 가능한 수문량의 규모를 파악하는데 있으며, 빈도해석의 정확도는 적합한 확률분포모형의 선택과 매개변수 추정방법에 의존하게 된다. 일반적으로 각 확률분포모형의 특성을 대표하는 매개변수를 추정하기 위해서는 모멘트 방법, 확률가중 모멘트 방법, 최대우도법 등을 이용하게 된다. 모멘트 방법에 의한 매개변수 추정은 해를 구하기 위한 과정이 단순한 반면, 비대칭형의 왜곡된 분포를 갖는 자료들에 대해서는 부정확한 결과를 나타내게 된다. 확률가중 모멘트 방법은 표본의 크기가 작거나 왜곡된 자료일 경우에도 비교적 안정적인 결과를 제공하는 반면, 확률 가중치가 정수로만 제한되는 단점을 갖고 있다. 그리고 대수 우도함수를 이용하여 매개변수를 추정하게 되는 최우도법은 가장 효율적인 매개변수 추정치를 얻을 수 있는 것으로 알려져 있으나, 비선형 연립방정식으로 표현되는 해를 구하기 위해서는 Newton-Raphson 방법을 사용하는 등 절차가 복잡하며, 때로는 수렴이 되지 않아 해룰 구하지 못하는 경우가 발생되게 된다. 이에 반해, 최근의 Genetic Algorithm, Ant Colony Optimization 및 Simulated Annealing과 같은 Meta-Heuristic Algorithm들은 복잡합 공학적 최적화 문제 있어서 효율적인 대안으로 주목받고 있으며, Hassanzadeh et al.(2011)에 의해 수문학적 빈도해석을 위한 매개변수 추정에 있어서도 그 적용성이 검증된바 있다. 본 연구의 목적은 연 최대강수 자료의 빈도해석에 적용되는 확률분포모형들의 매개변수 추정을 위해 Meta-Heuristic Algorithm을 적용하고자 함에 있다. 따라서 본 연구에서는 매개변수 추정을 위한 방법으로 Genetic Algorithm 및 Harmony Search를 적용하였고, 그 결과를 최우도법에 의한 결과와 비교하였다. GEV 분포를 이용하여 Simulation Test를 수행한 결과 Genetic Algorithm을 이용하여 추정된 매개변수들은 최우도법에 의한 결과들과 비교적 유사한 분포를 나타내었으나 과도한 계산시간이 요구되는 것으로 나타났다. 하지만 Harmony Search를 이용하여 추정된 매개변수들은 최우도법에 의한 결과들과 유사한 분포를 나타내었을 뿐만 아니라 계산시간 또한 매우 짧은 것으로 나타났다. 또한 국내 74개소의 강우관측소 자료와 Gamma, Log-normal, GEV 및 Gumbel 분포를 이용한 실증연구에 있어서도 Harmony Search를 이용한 매개변수 추정은 효율적인 매개 변수 추정치를 제공하는 것으로 나타났다.

  • PDF

폐합형 수계에 대한 수리학적 계산모형 개발 (Development of Computational Model for Looped Network Channel)

  • 구강민;전경수
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2015년도 학술발표회
    • /
    • pp.387-387
    • /
    • 2015
  • 최근 기후변화에 따른 국지성 집중호우 및 돌발홍수 증가로 도심지역에 많은 침수 피해가 발생하고 있다. 본 연구는 이러한 내외수에 영향을 받는 도시 중소하천의 침수 저감을 위하여 침수 예측모형에 적용 가능한 수리학적 계산모형을 개발하는 것이다. 일반적으로 자연하천은 수지형 수계가 대부분이며, 도시 우수관로와 관계배수 시스템 등은 폐합형 수계에 포함되어 수지형 계산모형으로는 도시 중소하천의 유출 모의를 할 수 없다. 폐합형 수계의 계산모형은 수지형에 비해 복잡하지만, 적용대상 수계가 폐합형이 아닌 경우에도 합류점의 유입량을 처리하기가 편리하고, 역방향의 월류 흐름이 존재하는 감조하천에서의 월류 흐름 모의가 가능한 장점을 갖고 있어 도시 내배수 시스템은 물론 자연하천에도 적용 할 수 있다. 본 모형은 절점, 수로 및 계산점으로 구성되는데 동력학적 방법인 1차원 Saint-Venant의 연속 방정식과 운동량 방정식에 수치해법을 이용하여 구하고자 하는 시간과 지점의 수위와 유량을 계산할 수 있게 구성하였다. 수치해법으로는 가장 보편적으로 사용되는 유한차분법 중 안정성과 정확성이 우수한 것으로 평가된 Preissmann의 4점 음해법으로 차분방정식에 Newton-Raphson 방법을 사용하여 유량과 수위 보정치에 Double Sweep 알고리즘을 적용하였다. 유한차분법은 안정성 문제를 수반할 수 있기에 시간 증분이 작을수록 Courant 조건을 만족할 수 있다. 모형 비교 검증을 위하여 동력학적 방법을 적용한 대표적인 폐합형 수계모형인 EPA SWMM을 지원하는 CHI사의 PC-SWMM 프로그램을 이용하여 가상의 폐합 수계를 구축하였다. 일반적으로 상류지점의 경계조건은 하류로 추적될 입력 자료로써, 상류지점과 합류지점은 유입되는 유량값을 그리고 하류지점은 유출되는 수위값을 경계조건으로 입력하였다. 운동량방정식의 에너지 경사와 마찰경사 항에 포함된 조도계수는 변화량에 중요한 물리적 요소이지만, 고정 상수값인 0.03을 적용하여 검증에 용의하도록 하였다. 구축된 모형과 PC-SWWM을 통해 산출된 계산점별 수위와 유량에 RMS 오차를 비교한 결과 만족할만한 결과를 얻을 수 있었다. 따라서 향후 내외수를 연계에 침수예측모형에 적용 가능할 것으로 사료된다.

  • PDF

발전용 가스터빈 동적 거동 시뮬레이션 Tool 개발 및 해석 (Development of Transient Behavior Simulation Tool and Analysis of Gas Turbines)

  • 김정호;김동섭
    • 플랜트 저널
    • /
    • 제13권4호
    • /
    • pp.48-50
    • /
    • 2017
  • 산업용 가스터빈의 동적 거동 해석을 위하여 시뮬레이션 Tool을 개발하였다. 시뮬레이션 Tool의 확장성을 향상시키기 위해 모든 가스 터빈 부품(압축기 및 연소기, 터빈, 덕트)을 모듈화하였다. 우리는 이 목적을 위해 객체 지향 프로그래밍을 사용했다. 질량 및 에너지 평형식은 다변수 뉴튼렙슨법을 사용하여 수치적인 해를 구했다. 가스터빈의 성능을 예측하기 위하여 압축기와 터빈의 성능선도를 사용하였다. 연소는 완전연소로 가정하였다. 가스터빈의 회전수와 터빈 배기 가스온도를 일정하게 유지하기 위해서 PID 제어를 사용하여 연료량과 압축기 입구 안내깃을 동시에 제어하였다. 가스터빈을 안정적으로 제어할 수 있었고 매우 빠른 부하 변화에도 대응이 가능함을 확인하였다.

  • PDF

Domain decomposition technique to simulate crack in nonlinear analysis of initially imperfect laminates

  • Ghannadpour, S. Amir M.;Karimi, Mona
    • Structural Engineering and Mechanics
    • /
    • 제68권5호
    • /
    • pp.603-619
    • /
    • 2018
  • In this research, an effective computational technique is carried out for nonlinear and post-buckling analyses of cracked imperfect composite plates. The laminated plates are assumed to be moderately thick so that the analysis can be carried out based on the first-order shear deformation theory. Geometric non-linearity is introduced in the way of von-Karman assumptions for the strain-displacement equations. The Ritz technique is applied using Legendre polynomials for the primary variable approximations. The crack is modeled by partitioning the entire domain of the plates into several sub-plates and therefore the plate decomposition technique is implemented in this research. The penalty technique is used for imposing the interface continuity between the sub-plates. Different out-of-plane essential boundary conditions such as clamp, simply support or free conditions will be assumed in this research by defining the relevant displacement functions. For in-plane boundary conditions, lateral expansions of the unloaded edges are completely free while the loaded edges are assumed to move straight but restricted to move laterally. With the formulation presented here, the plates can be subjected to biaxial compressive loads, therefore a sensitivity analysis is performed with respect to the applied load direction, along the parallel or perpendicular to the crack axis. The integrals of potential energy are numerically computed using Gauss-Lobatto quadrature formulas to get adequate accuracy. Then, the obtained non-linear system of equations is solved by the Newton-Raphson method. Finally, the results are presented to show the influence of crack length, various locations of crack, load direction, boundary conditions and different values of initial imperfection on nonlinear and post-buckling behavior of laminates.

Ultimate strength estimation of composite plates under combined in-plane and lateral pressure loads using two different numerical methods

  • Ghannadpour, S.A.M.;Shakeri, M.;Barvaj, A. Kurkaani
    • Steel and Composite Structures
    • /
    • 제29권6호
    • /
    • pp.785-802
    • /
    • 2018
  • In this paper, two different computational methods, called Rayleigh-Ritz and collocation are developed to estimate the ultimate strength of composite plates. Progressive damage behavior of moderately thick composite laminated plates is studied under in-plane compressive load and uniform lateral pressure. The formulations of both methods are based on the concept of the principle of minimum potential energy. First order shear deformation theory and the assumption of large deflections are used to develop the equilibrium equations of laminated plates. Therefore, Newton-Raphson technique will be used to solve the obtained system of nonlinear algebraic equations. In Rayleigh-Ritz method, two degradation models called complete and region degradation models are used to estimate the degradation zone around the failure location. In the second method, a new energy based collocation technique is introduced in which the domain of the plate is discretized into the Legendre-Gauss-Lobatto points. In this new method, in addition to the two previous models, the new model named node degradation model will also be used in which the material properties of the area just around the failed node are reduced. To predict the failure location, Hashin failure criteria have been used and the corresponding material properties of the failed zone are reduced instantaneously. Approximation of the displacement fields is performed by suitable harmonic functions in the Rayleigh-Ritz method and by Legendre basis functions (LBFs) in the second method. Finally, the results will be calculated and discussions will be conducted on the methods.

IoT 장비에 있어서 실시간 데이터 압축 전송을 위한 BL-beta 유니버설 코드의 경량화, 고속화 연구 (The study on Lightness and Performance Improvement of Universal Code (BL-beta code) for Real-time Compressed Data Transferring in IoT Device)

  • 김정훈
    • 한국정보전자통신기술학회논문지
    • /
    • 제15권6호
    • /
    • pp.492-505
    • /
    • 2022
  • 본 연구는 IoT 센싱 데이터의 무 손실 실시간 전송에 활용 가능한 BL-beta 코드의 엔코딩 및 디코딩 성능 개선을 통해 효과적으로 압축 데이터를 실시간으로 전송하고, 해독할 수 있도록 로직을 개선한 결과에 대한 연구이다. BL-beta 코드의 엔코딩 과정에는 비교적 연산 부담이 큰 로그 함수와 지수 함수, 나눗셈 및 제곱근 연산 등이 포함되어 있는데 이를 개선하여 비트 연산과 이진수 패턴 분석 그리고 비트 패턴을 이용한 뉴턴-랩슨 방법의 초기 값 설정을 통해 빠르게 데이터를 BL-beta 코드로 엔코딩 및 디코딩 할 수 있는 새로운 규칙성을 발견하였으며 이를 적용하여 기존 연구와 비교하여 알고리즘의 엔코딩 속도를 평균 24.8%, 디코딩 속도를 평균 5.3% 개선하였다.

레일 파단시 장대레일 개구부에서의 궤도-차량 동적상호작용 (Dynamic Interaction of Track and Train System on Open Gap by Rail Breaks)

  • 강윤석;강영종;양신추
    • 대한토목학회논문집
    • /
    • 제28권6D호
    • /
    • pp.895-904
    • /
    • 2008
  • 장대레일은 동절기에 중위온도 이하로 떨어져 레일의 인장파단에 의해 개구부가 발생할 수 있다. 레일 파단을 미리 감지 못하는 경우에 개구부에서 동적하중에 의해서 차륜과 레일의 파손이 발생할 수 있으며, 열차주행 안전성이 감소될 수 있다. 본 연구에서는 개구부에서의 궤도와 차량의 동적거동 분석을 위해 레일 파단시 개구부에서의 궤도와 차량의 상호작용을 고려한 동적해석모델을 제안하였다. 궤도와 차량의 선형시스템은 비선형 헤르찌안 접촉 스프링에 의해 연성되었으며 전체 궤도-차량 시스템 운동방정식을 정식화하였다. 그리고 개구부에서의 상호작용 현상을 고려하여 궤도 불균일부의 함수를 정의하고 개구부에서의 전후방 레일 사이의 개구량을 고려하였다. 비선형 방정식을 풀기위해 동적해석은 $Newmark-{\beta}$ 방법에 의해 수행되었다. 그리고 차량속도, 개구량, 레일지지강성에 따른 매개변수에 대한 궤도-차량의 동적상호작용해석을 수행하였고 변수에 대한 영향을 분석하였다.

평면(平面) 트러스 구조물(構造物)의 형상최적화(形狀最適化)에 관한 구연(究研) (A Study on Shape Optimization of Plane Truss Structures)

  • 이규원;변근주;황학주
    • 대한토목학회논문집
    • /
    • 제5권3호
    • /
    • pp.49-59
    • /
    • 1985
  • 탄성(彈性) 이론(理論)에 의하여 트러스의 형상최적화(形狀最適化) 문제(問題)를 형성(形成)하게 되면 부재(部材)의 단면적(斷面積)과 절점(節點)의 좌표(座標)를 동시에 고려(考慮)해야 하는 복잡(複雜)한 비선형(非線型) 계획문제(計劃問題)가 된다. 이런 비선형(非線形) 계획문제(計劃問題)를 해석(解析)할 수 있도록 제시(提示)된 기법(技法)이 별로 없고 현재 사용(使用)하고 있는 기법(技法)들도 실제(實際)의 적용(適用)에 제한(制限)을 받는 경우가 많다. 그러므로 트러스의 형태(形態), 재하조건(載荷條件) 등에 구애됨이 없이 트러스의 형상(形狀)을 최적화(最適化)할 수 있는 일반(一般) 해석기법(解析技法)이 필요(必要)한 것이다. 이에 본연구(本硏究)에서는 전(全) 해석과정(解析過程) two-phases로 나누어 phase 1 에서는 단면(斷面)을 최적화(最適化)하고 phase 2 에서는 트러스의 절점좌표(節點座標)를 변수(變數)로 하여 형상(形狀)을 최적화(最適化)하는 알고리즘을 개발(開發)한 것이다. 이 알고리즘의 phase 1 에서 유도(誘導)된 비선형(非線型) 계획문제(計劃問題)를 SUMT 문제(問題)로 변환(變換)시켜 Modified Newton-Raphson Method에 의한 SUMT 법(法)을 채택(採擇)하고 phase 2 에서는 Rosenbrock Method의 일방향(一方向) 탐사기법(探査技法)에 의해 목적함수(目的凾數)만이 최소(最小)가 되도록 하는 기법(技法)을 도입(導入)하여 최적화(最適化) 알고리즘 개발(開發)하였다. 개발(開發)된 알고리즘을 트러스의 형태(形態), 설계제약조건(設計制約條件), 재하조건(載河條件) 등을 변화(變化)시켜 가면서 수종(數種)의 트러스에 적용(適用)하여 수치계산(數値計算)을 실시(實施)하고 그 결과(結果)를 다른 알고리즘의 결과(結果)와 정교(正較)하므로서 개발(開發)된 알고리즘의 타당성(妥當性) 안정성(安定性) 적용성(適用性)을 검토(檢討)하였다. 연구(硏究) 결과(結果) 개발(開發)된 이 two-phases 알고리즘은 트러스의 설계조건(設計條件)에 구애받지 않고 트러스의 형상최적화(形狀最適化)에 적용(適用)할 수 있으며 안정성(安定性)있게 빠른 속도(速度)로 최적해(最適解)에 수렴(收斂)한다는 사실(事實)이 확인(確認)되었다. 이에 본(本) 알고리즘을 트러스의 형상최적화(形狀最適化) 알고리즘으로 새로이 제안(提案)하고 본(本) 알고리즘이트러스의 경제적(經劑的)인 설계(設計)에 도움을 줄 수 있을 것으로 사료(思料)된다.

  • PDF