• Journal of the Korean Society of Propulsion Engineers
• /
• v.3 no.3
• /
• pp.31-39
• /
• 1999

In this paper, numerical experiment is attempted to analyze and compare the combustion efficiency of the burning sprays due to OFO, FOF triplet / FOOF split doublet injectors. Preconditioned Wavier-Stokes equation system with low Reynolds number $\kappa$-$\varepsilon$ model for turbulence closure, is LU-SGS time-integrated. Spray processes are modeled by DSF analysis with experimentally determined injection characteristics. n-heptane/air global reaction model approximates the combustion for simplicity, and the influence of turbulence on the chemical reaction is included using eddy dissipation model. The results showed the FOF triplet injector of highest combustion efficiency, whereas the OFO type of poet performance. It was also observed that the droplet mean diameter and the average gas temperature due to the mixing efficiency, are the representative parameters for the performance design of combustion.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.2
• /
• pp.70-76
• /
• 2000

A numerical study on the viscous flow around a 2-dimensional circulation control foil is carried out for application on the field of naval architecture and ocean engineering. The governing equations are the RANS and the continuity equations. The equations are discretized by finite difference method and MAC method and the pressure poisson equation is calculate by a SOR method and an O-type non-staggered boundary fitted coordinate system which is overlapped near the slot is used to improve the numerical accuracy. Turbulence is approximated by a modified Baldwin-Lomax turbulence model. In the present paper, the Coanda effect on a 2-dimensional foil of a 20% thickness ellipse with modified rounded trailing edge has been numerically studied. The change in drag and lift of the foil with various jet momentums are calculated and compared to the experimental results to show good agreements.

• The Korean Journal of Rheology
• /
• v.11 no.1
• /
• pp.18-27
• /
• 1999

Predicting the deformation behaviors of sheets in thermoforming processes has been a daunting challenge due to the strong nonlinearities arising from very large deformations, mold-polymer contact condition and hyperelasticity constitutive equations. Nonlinear numerical analysis is always required to face this challenge especially for realistic processing conditions. In this study a 3-D algorithm and the membrane approximation are developed for thermoforming processes. The constitutive equation is expressed in terms of the 2nd Piola-Kirchhoff stress tensor and the Cauchy-Green deformation tensor. The 2-term Mooney-Rivlin model is used for the material model equation. The algorithm is established by the finite element formulation employing the total Lagrangian coordinate. The deformation behavior and the stress distribution results of 3-D algorithm with various point boundary conditions are compared to those of the membrane approximation algorithm. Also, the slip boundary condition and the no-slip boundary condition are applied for the systems that have molds. Finally, the effect of sheet temperatures on the final thickness distribution is investigated for the ABS material.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2006.11a
• /
• pp.363-366
• /
• 2006

The effectiveness of two passive control techniques for alleviating the pressure oscillation generated in a supersonic cavity flow is investigated numerically. The passive devices suggested in the present research include a triangular bump and a sub-cavity installed near the upstream edge of a rectangular cavity. The supersonic cavity flow characteristics are examined by using the three-dimensional, unsteady Wavier-Stokes computation based on a finite volume scheme. Large eddy simulation (LES) is carried out to properly predict the turbulent features of cavity flow. The results show that the pressure oscillation near the downstream edge dominates overall time-dependent cavity pressure variations. Such an oscillation is attenuated more considerably using the sub-cavity compared with other methods, and a larger sub-cavity leads to better control performance.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.40 no.8
• /
• pp.647-654
• /
• 2012

Supercooled large droplet issues in aircraft icing have been continually reported due to the important safety considerations. In order to simulate the impingement behavior of large droplets, a two-dimensional and compressible Navier-Stokes code was developed to determine the flow field around the test model. Also, the Eulerian-based droplet impingement model including a semi-empirical approach for the droplet-wall interaction process and droplet break-up was developed. In particular, the droplet-wall interactions were considered as numerical boundary conditions for the droplet impingement simulation in the supercooled large droplet conditions. Finally, the present results were compared with the experimental test data and the LEWICE results. The droplet impingement area and maximum collection efficiency values between present results and wind tunnel data were in good agreements. Otherwise, the inclination of collection efficiency of the present result is over-predicted than the wind tunnel data around a lower surface of the NACA 23012 airfoil.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.1
• /
• pp.16-27
• /
• 2002

A level-set method is used for analyzing the behaviors of gas bubbles in two fluids incompressible viscous flow domain. The governing equations are solved by using a finite volume method. The numerical results are verified by comparing with the experimental and other computational results. Computations for the deformations and motions of one or multi-bubbles in the flow domain with the initial undisturbed free interface are conducted. It can be seen that numerical results for different surface tension and density ratio arise very different behaviors of bubbles. When bubbles rise near the free interface, the free interface gives some great influence on the behaviors of bubbles. The present results computed by a level-set method give useful information about the properties of bubble motions and deformations.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.17 no.3
• /
• pp.158-165
• /
• 2005

In order to analyze incompressible two-phase flow, Level Set method was applied on pseudo-compressibility formulation. Level Set function is defined as a signed distance function from the phase interface, and gives the information of the each phase location and the geometric data to the flow. In this study, Level Set function transport equation was coupled with flow conservation equations, and owing to pseudo-compressibility technique we could solve the resultant vector equation iteratively. Two-phase flow analysis code was developed on general curvilinear coordinate, and numerical tests of bubble dynamics and surging wave problems demonstrate its capability successfully.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 1998.10a
• /
• pp.3-3
• /
• 1998

액체추진제 로켓엔진 연소실에는 고유모드에 대응하는 음향파동이 내재되며 이러한 음향파동은 연소와의 상호작용을 통하여 불안정한 음향에너지를 공급받아 증폭되며 결국에는 연소불안정 상태에까지 이르게 된다. 이와 같은 불안정한 상태에 이르기 위해서는 연소로부터 되먹임되는 불안정 에너지의 양이 충분히 크고 구동 음향파동에 근접한 위상을 가져야 한다. 이와 같은 구동 메커니즘을 구성하는 상세한 물리적 현상들을 규명하고 예측하기 위한 많은 연구들이 보고되었으며, 이들 중 이론적인 시간 지연 모델을 사용하는 음향적인 방법은 매우 경제적인 반면 연소 현상에 대한 상세한 모사가 생략되어 연소 불안정의 구체적인 원인을 규명하는데 어려움이 있고, 파동 방정식에 의하여 연소실 내부의 파동 에너지 증가를 예측하는 방법은 연소기 내에서의 연소 메커니즘에 대한 고려 없이 연소에 의해 발생하는 에너지만을 포함하는 단점과 선형적인 연소 불안정에만 제한된다는 제한이 있다. 음향장과 커플된 기화반응 모델은 분무액적의 기화 과정이 추진제 연소의 지배과정이라는 가정 하에 연소응답을 기화반응으로 대체하는 방법으로, 역시 단시간 내에 결과를 얻을 수 있다는 장점이 있으나 기화반응으로부터 음향파동으로의 에너지 되먹임 과정이 배제되어 있어 정확한 결과를 구하기는 어렵다. 이에 대하여 최근에는 전산 모사적인 방법을 사용하는 대규모의 연소장 해석이 가능하여 짐으로써 음향파동에 의한 외란과 에너지 되먹임과정을 모두 포하마여 수치적인 방법을 사용하여 계산하는 액체추진제 로켓엔진의 고주파 연소불안정 해석방법들이 제시되고 있다.안정성 모드가 있음을 보였다. 밀도 변화가 있는 경우나 밀도 변화가 없는 경우 모두 sinuous 모드의 가장 불안정한 모드가 varicose 모드의 가장 불안정한 모드보다 더 불안정함을 보여주어 후류 유동은 자유 유동에 가까운 위상 속도를 가지는 sinuous 모드에 의해 지배될 것임을 예측할 수 있다. 연소반응이 완전연소에 가까울수록 그리고 압축성 효과가 클수록 유동내부의 온도가 증가하고 점성 또한 증가하여 후류유동은 안정됨을 알 수 있었다 유동변수들의 contour로부터 유동의 특성을 예측한 결과 baroclinic 항이 dilatational 항보다 상대적으로 크며, 중심선 상하에 생기는 vortex를 더욱 성장시킬 것으로 생각된다.냉각 홀의 막임, 연소 입자의 점착 부위 등을 예측하여 보완책을 준비할 수 있도록 하였다.$mm^2$sec였으며, 이는 다른 graphite/epixy 복합재의 확산계수와 유사한 값을 나타내고 있다. 또한 추진제가 충전된 연소관을 절단하여 밀폐한 후 95%RH 습도 조건에 보관함으로써 연소관 내부의 추진제 기계적 특성에 미치는 침투된 습기의 영향도 함께 고찰하였다. 추진제에 따라 차이는 있겠으나 추진제가 충전된 연소관은 순수 복합재 연소관에 비해 습기의 투과 정도가 작으며, 본 연소관에 충전된 RDX/AP계 추진제의 경우 추진제의 습기투과에 의한 추진제 물성 변화는 미미한 것으로 나타났다.의 향상으로, 음성개선에 효과적이라고 사료되었으며, 이 방법이 편측 성대마비 환자의 효과적인 음성개선의 치료방법의 하나로 응용될 수 있으리라 생각된다..

• The Journal of Engineering Geology
• /
• v.15 no.3
• /
• pp.337-348
• /
• 2005

Jurassic granite from Pocheon area were tested to investigate the effect of microcracks on mechanical properties of the granite. Three oriented core specimens were used for uniaxial compressive tests and each core specimen are perpendicular to the axes'R'(rift plane),'c'(grain plane) and'H'(hardway plane), respectively Among vacious elastic constants, the variation of Poisson's ratio as function of the directions was examined. From the related chart between ratio of failure strength and Poisson's ratio, H-specimen shows the highest range in Poisson's ratio and Poisson's ratio decreases in the order of C-specimen and R-specimen. The curve pattern is nearly linear in stage $I\simIII$ but the slope increases abruptly in stage H-3. As shown in the related chart, diverging point of a curve is formed when ratio of failure strength is $0.92\sim0.96$ Stage IV -3 is out of elastic region. The behaviour of rock in the four fracturing stages was analyzed in term of the stress-volumetric strain me. From the stress increment-volumetric strain equations governing the behaviour of rock, characteristic material constants, a, n, Q, m and $\varepsilon_v^{mcf}$, were determined. Among these, inherent microcrack porosity$(a, 10^{-3})$ and compaction exponent(n) in the microcrack closure region(stage I ) show an order of $a^R(3.82)>a^G(3.38)>a^H(2.32)\;and\;n^R(3.69)>n^G(2.79)>n^H(1.99)4, respectively. Especially, critical volumetric microcrack strain($\varepsilon_v^{mcf}$) in the stage W is highest in the H-specimen, normal to the hardway plane. These results indicate a strong correlation between two major sets of microcracks and mechanical properties such as Poisson's ratio and material constants. Correlation of strength anisotropy with microcrack orientation can have important application in rock fracture studies.