• Title/Summary/Keyword: 근사해석법

Search Result 523, Processing Time 0.023 seconds

BGK 수치기법을 이용한 로켓 노즐 내의 유동장 해석

  • 신동신;이재성
    • Proceedings of the Korean Society of Propulsion Engineers Conference
    • /
    • 2000.04a
    • /
    • pp.16-16
    • /
    • 2000
  • 충격파를 포함하는 초음속 유동을 해석하는 수치해법 중에서 많이 사용되어진 것은 엄밀 및 근사 리만 해법과 플럭스 분할 기법들로서 이들은 Euler 방정식에 기반을 두고 선형 또는 비선형파의 상호작용을 풍상 차분법으로 기술하는 방법들이다. 이러한 수치기법들은 과거 광범위하게 사용되어 왔으나 최근 여러 가지 단점이 발견되었다. 이와 같은 문제점을 극복하고자 입자의 통계적인 운동을 기술하는 기체 운동론에 근거하여 BGK 수치기법이 제시되었다. 이는 비충돌 볼츠만 방정식으로부터 입자의 수준에서 플럭스 분할 기법 형태의 풍상차분법을 구현하는 것으로 볼츠만 방정식의 충돌항을 BGK 모델로 대치하고 이것의 적분해로부터 수치 플럭스를 구한다. 이 수치기법은 기존의 리만해법에 비하여 수치적으로나 물리적으로 매우 타당한 성질인 강건성, 정확성, 엔트로피 조건, 양수보존성 등을 가지고 있음이 밝혀졌다. 이와 같은 수치기법을 사용하여 로켓 노즐 내의 아음속, 천이음속, 초음속에서의 유동장 해석을 위한 프로그램을 작성하였다. 시간 적분에 대하여는 정상 상태의 계산을 위하여 내재적 시간 적분 방법을 사용하였으며, 공간 이산화 방법으로는 임의의 제어체적에 대하여 적분형 보존 방정식을 적용하는 유한 체적법을 사용하였다. 초음속 입구 유동과 출구에서 초음속과 저음속 유동의 두가지 경우를 고려하여 얻은 결과를 기존의 연구 결과와 비교하여 본 결과 잘 일치하였다. 입구 유동이 저음속이고 출구 유동이 초음속인 경우에 대하여도 해석결과가 실험결과와 잘 일치하였다. 상대적으로 낮은 온도, 압력 조건과 높은 온도, 압력 조건을 가지는 고체 로켓 모터 노즐 내의 유동을 해석하였다. 이들 해석 결과를 전압, 전온도로 표준화시킨 결과 서로 일치하였으며, 파라서 저온, 저압에서 얻은 결과도 표준화시킬 경우, 고온, 고압에서도 사용될 수 있음을 알 수 있었다.의 영향에 초점을 맞추었다.다고 판단되며 배기 가스 자체에 대기 공기중에 함유되어 있던 습기가 얼어붙는(Icing화) 문제가 발생하기 때문에 배기가스의 Icing을 방지하기 위하여 압축기 끝단에서 공기를 추출하여 배기부분에 송출할 필요성이 있는 것으로 판단되었다. 출구가스의 기체 유동속도가 매우 빠르므로 (100-l10m.sec) 이를 완화하기 위한 디퓨저의 설계가 요구된다고 판단된다. 또 연소기 후방에 물을 주입하는 경우 열교환기 및 기타 부분품에 발생할 수 있는 부식 및 열교환 효율 저하도 간과할 수 없는 문제로 파악되었다. 이러한 기술적 문제가 적절히 해결되는 경우 비활성 가스 제너레이터는 민수용으로는 대형 빌딩, 산림, 유조선 등의 화재에 매우 적절히 사용되어 질 수 있을 뿐 아니라 군사적으로도 군사작전 중 및 공군 기지의 화재 그리고 지하벙커에 설치되어 있는 고급 첨단 군사 장비 등의 화재 뿐 아니라 대간첩작전 등에 효과적으로 활용될 수 있을 것으로 판단된다.가 작으며, 본 연소관에 충전된 RDX/AP계 추진제의 경우 추진제의 습기투과에 의한 추진제 물성 변화는 미미한 것으로 나타났다.의 향상으로, 음성개선에 효과적이라고 사료되었으며, 이 방법이 편측 성대마비 환자의 효과적인 음성개선의 치료방법의 하나로 응용될 수 있으리라 생각된다..7%), 혈액투석, 식도부분절제술 및 위루술·위회장문합술을 시행한 경우가 각 1례(2.9%)씩이었다. 13) 심각한 합병증은 9례(26.5%)에서 보였는데 그중 식도협착증이 6례(17.6%), 급성신부전증 1례(2.9%), 종격동기흉과 폐염이 병발한 경우와 폐염이 각 1례(2.9%)였다. 14) 식도경 시행회수는 1회가 17례(54.8%), 2회가 9례(29.0%), 3회 이상이 5례(16.1%)였다.EX>$IC_{50}$/ 값이 210 $\mu\textrm{g}$<

  • PDF

Stress Analysis of the Corner Part of Ship Structures Using the New Equivalent Curved Beam Theory (신(新) 등가(等價) 곡선(曲線)보 이론(理論)에 의한 선체(船體) Corner부(部)의 응력(應力) 해석(解析))

  • Chang-Doo Jang;Seung-Soo Na
    • Journal of the Society of Naval Architects of Korea
    • /
    • v.29 no.4
    • /
    • pp.193-201
    • /
    • 1992
  • A new equivalent curved beam theory is developed for the analysis of the corner part of ship structures, in which effects of distributed loads and asymmetricity with two or three connected parts are considered. Equivalent loads are obtained from equilibrium conditions between the distributed loads and the member forces and moments at the ends of curved beam. And an equivalent curved beam for the asymmetric structure is obtained by superposing the equivalent symmetric parts which have equivalent stiffness. From the sample calculation, it is found that the results of the new equivalent curved beam theory are well agreed with those of finite element method using membrane elements with little computing time and sufficient accuracy.

  • PDF

Analysis of the thermal fluid flow between the gas torch and the steel plate for the application of the line heating (선상 가열을 위한 가스 토치와 강판 사이의 열유동 해석)

  • Jong-Hun Woo;Jong-Gye Shin
    • Journal of the Society of Naval Architects of Korea
    • /
    • v.39 no.2
    • /
    • pp.52-60
    • /
    • 2002
  • Line heating is a forming process which makes the curved surface with the residual strain created by applying heat source of high temperature to steel plate. in order to control the residual strain, it is necessary to understand not only conductive heat transfer between heat source and steel plate, but also temperature distribution of steel plate. In this paper we attempted to analyze is temperature distribution of steel plate by simplifying a line heating process to collision-effusive flux of high temperature and high velocity, and conductive heat transfer phenomenon. To analyze this, combustion in the torch is simplified to collision effusive phenomenon before analyzing turbulent heat flux. The distribution of temperature field between the torch and steel plate is computed through turbulent heat flux analysis, and the convective heat transfer coefficient between effusive flux and steel plate is calculated using approximate empirical Nusselt formula. The velocity of heat flux into steel plate is computed using the temperature distribution and convective heat transfer coefficient, and temperature field in the steel plate is obtained through conductive heat transfer analysis in which the traction is induced by velocity of heat flux. In this study, Finite Element Method is used to accomplish turbulent heat flux analysis and conductive heat transfer analysis. FEA results are compared with empirical data to verify results.

Analysis of Electromagnetic Wave Scattering From a Perfectly Conducting One Dimensional Fractal Surface Using the Monte-Carlo Moment Method (몬테칼로 모멘트 방법을 이용한 1차원 프랙탈 완전도체 표면에서의 전자파 산란 해석)

  • 최동묵;김채영
    • Journal of the Institute of Electronics Engineers of Korea TC
    • /
    • v.39 no.12
    • /
    • pp.566-574
    • /
    • 2002
  • In this paper, the scattered field from a perfectly conducting fractal surface by the Monte-Carlo moment method was computed. An one-dimensional fractal surface was generated by using the fractional Brownian motion model. Back scattering coefficients are calculated with different values of the spectral parameter(S$\_$0/), and fractal dimension(D) which determine characteristics of the fractal surface. The number of surface realization for the computed field, the point number, and the width of surface realization are set to be 80, 2048, and 64L, respectively. In order to verify the computed results these results are compared with those of small perturbation methods, which show good agreement between them.

Evaluation of the Cochlear Electrode Behavior in the 3D Human Cochlea Model by FEM Analysis (3차원 달팽이관 모델에서의 인공와우 전극의 거동에 대한 유한요소해석)

  • 임윤섭;박세익;김용협;오승하;김성준
    • Journal of Biomedical Engineering Research
    • /
    • v.25 no.3
    • /
    • pp.207-215
    • /
    • 2004
  • A finite element analysis is used in this study to model 6 kinds of electrode by a genetic program in order to evaluate the mechanical effect on the 3D human cochlea model and the behavior of electrode. Human cochlea is modeled by the spiral-approximation method and the shape of scala tympani is extracted from the mid section of the human cochlea. Contact pressure at the tip and the insertion force are found to be highest when the wires stack horizontally. Axial rotation of electrode is minimal comparing with the stimulating current spread. The results indicate that the electrode stiffness is important to minimize the trauma.

Preliminary Design Conditions for a Thermally Actuated Refrigerator Based on the Vuilleumier Cycle (Vuilleumier 사이클로 작동되는 열구동 냉동기의 예비설계조건)

  • 유호선;강병하
    • Transactions of the Korean Society of Mechanical Engineers
    • /
    • v.16 no.12
    • /
    • pp.2358-2367
    • /
    • 1992
  • This paper deals with preliminary design conditions for a thermally actuated Vuilleumier refrigerator/heat pump. The previously reported approximate adiabatic analysis which is based on the 8-volume model makes it possible to evaluate exchanged heats per cycle as well as cyclic pressure, temperature and mass variation of each working volume. Calculated results reveal not only there exists an optimum value for the phase angle and the swept volume ratio maximizing a specific thermal output, but also design parameters can be determined independently of each other. Under a given combination of operating temperature levels, the optimum conditions for refrigeration are different from those for heat pumping and the differences between two operating modes become larger with decreasing the dead volume ratio. Both the optimum phase angle and the optimum swept volume ratio are increased asymptotically toward 0.5 pi and 1.0 respectively, as the dead volume ratio approaches to unity. When a VM machine is used for cooling and heating simultaneously, the design parameters should be carefully determined to reach the best performance.

Novel Polarization Condition of Optical Polarized Beamsplitter with Misaligned Input Interface and Guiding Loss (불연속 입력 경계면과 전송 손실이 있는 광 편파기 설계를 위한 새로운 편파조건)

  • Ho, Kwang-Chun
    • Journal of IKEEE
    • /
    • v.11 no.1 s.20
    • /
    • pp.61-68
    • /
    • 2007
  • Using modal transmission-line theory (MTLT), itexplores rigorously the polarization characteristics of TE-TM mode optical polarized beamsplitter with misaligned input interface. To verify the validity of the approach proposed in this paper, a novel polarization condition is defined, and it compares and evaluates with the polarization condition of coupled-mode theory (CMT) used widely. The numerical results reveal that the discontinuity of misaligned input interface has nothing to do with the polarization condition of beamsplitter. Furthermore, the condition of CMT is no longer effective at the asymmetric guiding structures with strong power coupling.

  • PDF

A Study on Stiffness-based Optimal Design of Tall Plane Frameworks using Composite Member (합성부재를 이용한 고층평면골조의 강성최적설계에 관한 연구)

  • Kim, Ho-Soo;Lee, Han-Joo
    • Journal of Korean Association for Spatial Structures
    • /
    • v.4 no.3 s.13
    • /
    • pp.77-84
    • /
    • 2004
  • This study presents an effective stiffness-based optimal technique to control quantitatively lateral drift for tall frameworks using composit member subject to lateral loads. To this end, displacement sensitivity depending on behavior characteristics of tall frameworks is established and approximation concept that preserves the generality of the mathematical programming and can efficiently solve large scale problems is introduced. Specifically, under the 'constant-shape' assumption, resizing techniqe of composite member is developed. Two types of 50 story frameworks are presented to illustrate the features of the quantitative lateral drift control technique proposed in this study.

  • PDF

Molecular Orbital Analyses on Hammett Substituent Constant (I) (Hammett 치환기 상수에 대한 분자궤도론적 해석 (I))

  • Byung-Kak Park;Gab-Yong Lee
    • Journal of the Korean Chemical Society
    • /
    • v.30 no.2
    • /
    • pp.172-180
    • /
    • 1986
  • The Hammett's substituent constants were analyzed for the substituted-benzoic acid, -benzene and -naphthalene by LCAO MO method. Hammett's ${\sigma}$-values have been theoretically proved to treat as the sum of independent contribution of inductive effect and resonance effect. That is, Is (static inductive effect) and Id (dynamic inductive effect) corresponding to inductive effect are chosen as net charge and approximate self atom polarizability indices respectively, and ${\Delta}E^{HOMO}$ which is the difference in HOMO energy level between substituted molecule and unsubstituted molecule as resonance effect. In conclusion, it has been found that the observed ${\sigma}$-values depend on the sum of Id, Is and ${\Delta}E^{HOMO}$.

  • PDF

Automatic Interpretation of the Borehole Normal Resistivity Data by Using a Personal Computer (퍼스널 컴퓨터를 이용한 비저항 물리검층자료의 자동해석)

  • Kim, Jin-Hu
    • Journal of Ocean Engineering and Technology
    • /
    • v.2 no.2
    • /
    • pp.51-60
    • /
    • 1988
  • A data transform is performed by using a point-electrode focusing method in order to obtain accurate and objective interpretation of the borehole normal resistivity data. Two new synthetic curves can be generated through the data transform. The one is an approximate apparent resistivity curve, which would be used to predict the true resistivity of the formation. The other one is a bed boundary coefficient curve, which would be used to distinguish bed boundaries. The accuracy of the normal data interpretation can be improved and this method takes much less computational time than a linear inversion technique. Moreover, this method does not require an initial guess model and limitation of number of unknown parameters. Since this algorithm can be run on a personal computer, an immediate interpretation would be possible at the field work site. If an additional set of electrodes(a=125cm)is attached to a normal resistivity tool which is being used (a=25cm, 50cm, 100cm), the apparent resistivity for the point-electrode focusing device can be calculated, and it would maximize the use of short and long normal resistivity data and promote the accuracy of the interpretation.

  • PDF