• 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계 추진제의 경우 추진제의 습기투과에 의한 추진제 물성 변화는 미미한 것으로 나타났다.의 향상으로, 음성개선에 효과적이라고 사료되었으며, 이 방법이 편측 성대마비 환자의 효과적인 음성개선의 치료방법의 하나로 응용될 수 있으리라 생각된다..

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.36 no.2
• /
• pp.161-170
• /
• 2012

There are many factors to consider when attempting to improve the efficiency of fuel cell operation, such as the operation temperature, humidity, stoichiometry, operation pressure, geometric features, etc. In this paper, the effects of the operation pressure were investigated to find the current density and water saturation behavior on a cross section designated by the design geometry. A two-dimensional geometric model was established with a gas channel that can provide $H_2$ to the anode and $O_2$ and water vapor to the cathode gas diffusion layer (GDL). The results from this numerical modeling revealed that higher operation pressures would produce a higher current density than lower ones, and the water saturation behavior was different at operation pressures of 2 atm and 3 atm in the cathode GDL. In particular, the water saturation ratios are higher directly below the collector than in other areas. In addition, this paper presents the dependence of the velocity behavior in the cathode on pressure changes, and the velocity fluctuations through the GDL are higher in the output area than in inlet area. This conclusion will be utilized to design more efficient fuel cell modeling of real fuel cell operation.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.26 no.4
• /
• pp.225-243
• /
• 2014

As the economic matters are involved, applying the WEC, which is used for controlling waves as well as utilizing the wave energy on existing breakwater, is preferred rather than installing exclusive WEC. This study examines the OWC mounted on a pressurized air chamber floating breakwater regarding the functionality of both breakwater and wave-power generation. In order to verify the performance as a WEC, the velocity of air flow from pressurized air chamber to WEC has to be evaluated properly. Therefore, numerical simulation was implemented based on BEM from linear velocity potential theory as well as Boyle's law with the state equation to analyze pressurized air flow. The validity of the obtained values can be determined by comparing the previous results from numerical analysis and empirically obtained values of different shapes. In the actual numerical analysis, properties of wave deformation around OWC system mounted on fixed type and floating type breakwaters, motions of the structure with air flow velocities are investigated. Since, the wind power generating system can be hybridized on the structure, it is expected to be applied on complex power generation system which generates both wind and wave power energy.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.6
• /
• pp.409-422
• /
• 2019

In this study, we carried out the 3D numerical simulation to investigate the hydraulic characteristics of Synchronous Edge wave known as the driving mechanism of beach cusp using the Tool Box called the ihFoam that has its roots on the OpenFoam. As a wave driver, RANS (Reynolds Averaged Navier-Stokes equation) and mass conservation equation are used. In doing so, we materialized short-crested waves known as the prerequisite for the formation of Synchronous Edge waves by generating two obliquely colliding Cnoidal waves. Numerical results show that as can be expected, flow velocity along the cross section where waves are focused are simulated to be much faster than the one along the cross section where waves are diverged. It is also shown that along the cross section where waves are focused, up-rush is moving much faster than its associated back-wash, but a duration period of up-rush is shortened, which complies the typical characteristics of nonlinear waves. On the other hand, due to the water-merging effect triggered by the redirected flow toward wave-diverging area at the pinacle of run-up, along the cross section where waves are diverged, offshore-ward velocity is larger than shore-ward velocity at the vicinity of shore-line, while at the very middle of shoaling process, the asymmetry of flow velocity leaned toward the shore is noticeably weakened. Considering that these flow characteristics can be found without exception in Synchronous Edge waves, the numerical simulation can be regarded to be successfully implemented. In doing so, new insight about how the boundary layer streaming occur are also developed.

• Geotechnical Engineering
• /
• v.13 no.5
• /
• pp.59-74
• /
• 1997

This paper shows the results of a series of model bests on the behavior of single steel pipe pile which is subjected to lateral load in Nak-dong river sand. The purpose of the present paper is to estimate the effect of Non -homogeneous soil, constraint condition of pile head, lateral load velocity, relative density of soil, embedded pile length, and flexural stiffness of pile on the behavior of single pile which is embedded in Nak-dong river strand. These effects can be quantined only by the results of model tests. The nonlinear responses of lateral loadieflection relationships are fitted to 2nd polynomial equations by model tests results. Also, the lateral load of a deflection, yield and ultimate lateral load max. bending moment, and yield bending moment can be expressed as exponential function in terms of relative density and deflection ratio. By comparing Brom's results with model results on the lateral ultimate load, it is found that short and long pile show the contrary results with each other. The contrary results are due to the smaller assumed soil reaction than the soil reaction of the Nakiong river sand at deep point. By comparing lateral behavior on the homogeneous soil with non-homogeneous soil, it is shown that lateral loadieflection relationship is very dependent on the upper relative density. This phenomenon is shown remarkably as the difference between upper and lower relative density increases.

• Journal of Advanced Marine Engineering and Technology
• /
• v.39 no.7
• /
• pp.709-715
• /
• 2015

Underwater gliders do not typically have separate propellers for forward motion. They generate propulsive forces based on the difference between their buoyancy and gravity. They can control the volume from the buoyancy engine to adjust the propulsive force. In addition, the attitude of the underwater glider is controlled by a rubberless motion controller. The motion controller can change the mass center and moment of inertia of the inner moving mass. Owing to the change in these parameters, the attitude of the underwater glider is changed. In this study, we derive nonlinear, six degree of freedom (DOF) mathematical models for the motion controller and buoyancy engine. Using these equations, we perform dynamic simulations of the proposed underwater glider, and verify the suitability of the design and dynamic performances of the proposed underwater glider. We then perform the motion control simulation for the pitch and roll angle, and analyze the dynamic performance according to the pitch and roll angles.

• Journal of Navigation and Port Research
• /
• v.39 no.1
• /
• pp.69-75
• /
• 2015

The concrete purpose of this study is to suggest actually a debt ratio to optimize the capital structure providing a kind of approach to estimate the proper debt ratio with an analytical model and empirical data in Korean shipping industry. The mathematical and analytical model is started from the first equation about ROE, return of net operating income on equity, with an independent variable, debt ratio. It is constructed with several parameters, ROS(return of operating income on sales), TAT(total assets turnover), and NFCL(net finance cost to liabilities). There could not be a certain relationship between debt ratio and ROS or TAT, while some correlation or causality between debt ratio and NFCL. In other words, most of firms with high debt ratio is likely to burden higher finance cost than others with low one. In this case, there is a linearity relationship between debt ratio and NFCL, so then the second equation considering this relation could be included within the analytical approach of this paper. To be short, if the criteria of adequate debt ratio has to be defined as some level of debt ratio to optimize ROE, the ROE could be illustrated as a quadratic equation to debt ratio from two equations. Next, this research estimated those parameters' numbers through the single regression method with data over 12 years of Korean shipping industry, and identified empirically the fact that optimal debt ratio would be approximately 400%. To conclude, if that industry's sales and operating incomes are stable, the debt ratio could be accepted until twice of 200% had forced in order to guarantee its financial safety in past time.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.3
• /
• pp.25-34
• /
• 1993

The long-term behavior, such as in excavation problems of weak medium, can be dealt with by the elasto-viscoplasticity models. In this paper, a combined formulation of elasto-viscoplasticity using boundary elements and finite elements without using internal cells is presented. The domain integral introduced due to the viscoplastic stresses is transformed into a boundary integral applying direct integration in cylindrical coordinates. The results of the developed boundary element analysis are compared with those from the explicit solution and from the finite element analysis. It is observed that the boundary element analysis without internal cells results in some error because of its deficiency in handling the nonlinearity in local stress concentration. Therefore, a coupled analysis of boundary elements and finite elements, in which finite elements are used in the area of stress concentration, is developed. The coupled method is applied to a time dependent inelastic problem with semi-infinite boundaries. It results in reasonable solution compared with other methods where relatively higher degree of freedoms are employed. Thus, it is concluded that the combined analysis may be used for such problems in the effective manner.

• Composites Research
• /
• v.26 no.6
• /
• pp.363-372
• /
• 2013

Numerical simulation for Resin Transfer Molding manufacturing process is attempted by using the eXtended Finite Element Method (XFEM) combined with the level set method. XFEM allows to obtaining a good numerical precision of the pressure near the resin flow front, where its gradient is discontinuous. The enriched shape functions of XFEM are derived by using the level set values so as to correctly describe the interpolation with the resin flow front. In addition, the level set method is used to transport the resin flow front at each time step during the mold filling. The level set values are calculated by an implicit characteristic Galerkin FEM. The multi-frontal solver of IPSAP is adopted to solve the system. This work is validated by comparing the obtained results with analytic solutions. Moreover, a localization method of XFEM and level set method is proposed to increase the computing efficiency. The computation domain is reduced to the small region near the resin flow front. Therefore, the total computing time is strongly reduced by it. The efficiency test is made with a simple channel flow model. Several application examples are analyzed to demonstrate ability of this method.

• Journal of Korea Water Resources Association
• /
• v.45 no.6
• /
• pp.557-568
• /
• 2012

In order to understand the hydrological processes on the mountainous forest, the configuration of soil evaporation (E) out of evapotranspiration (ET) is a challenging and important topic. In this study, we attempted to understand the soil evaporation process for a humid forest hillslope via measuring and analyzing soil moistures with a sampling interval in 2 hours at three locations for 10 days between May 22th and 31th 2009. Two methods were used to estimate soil evaporation in every 2hr; one is a method using soil moisture measurement ($E_{SM}$), the others methods are based on Penman equation (Penman (1948), Staple (1974), Konukcu (2007), Equilibrium Penman ($E_{equili}$)). As a critical parameter in determining $E_{SM}$, the dry surface layer (DSL), was estimated using energy balance equation. The accumulated soil evaporation ($E_{SM}$) of A, B, C points were estimated as 2.09, 1.08 and 2.88 mm, respectively. The estimated evaporation of Penman (1948), Staple (1974), Konukcu (2007), $E_{equili}$ were 4.91, 8.80, 8.63 and 3.28 mm. The proposed method with soil moisture measurement showed lower soil evaporations than the other conventional methods. The increasing soil temperature and interaction between soil and atmosphere due to existence of litter and DSL are considered as dominant factors for soil evaporation. The $E_{SM}$ has the apparent lag time between 2 and 4 hr compared with $E_{equili}$ and net radiation. The DSL and surface resistance ($r_s$) were increased as soil moisture was decreased for in this study. The estimated DSL through the temporal distribution analysis of soil moisture and tension measurements was also similar to that of the energy balance relationship.