• Journal of the Korean Society of Propulsion Engineers
• /
• v.22 no.3
• /
• pp.46-52
• /
• 2018

A simple Crocco's $n-{\tau}$ time delay model and linear analysis of fluid flow coupled with acoustics are combined to investigate the high frequency combustion instability in the combustion chamber of LOX/hydrocarbon engines. The partial differential equation of the velocity potential is separated into ordinary differential equations, and eigenvalues that correspond to tangential resonance modes in the cylindrical chamber are determined. A general solution is obtained by solving the differential equation in the axial direction, and boundary conditions at the injector face and nozzle entrance are applied in order to calculate the chamber admittance. Frequency analysis of the transfer function is used to evaluate the stability of system. Stability margin is determined from the system gain and phase angle for the desired frequency range of 1T mode. The chamber model with variable baffle length and configurations are also considered in order to enhance the 1T mode stability of the combustion chamber.

• Journal of the Korean Institute of Gas
• /
• v.7 no.4 s.21
• /
• pp.44-52
• /
• 2003

In this work dynamic heat transfer in a CPFS (cable penetration fire stop) system built in the firewall of nuclear power plants is three-dimensionally investigated to develop a test-simulator that can be used to verify effectiveness of the sealant. Dynamic heat transfer in the fire stop system is formulated in a parabolic PDE (partial differential equation) subjected to a set of initial and boundary conditions. First, the PDE model is divided into two parts; one corresponding to heat transfer in the axial direction and the other corresponding to heat transfer on the vertical planes. The first PDE is converted to a series of ODEs (ordinary differential equations) at finite discrete axial points for applying the numerical method of SOR (successive over-relaxation) to the problem. The ODEs are solved by using an ODE solver In such manner, the axial heat flux can be calculated at least at the finite discrete points. After that, all the planes are separated into finite elements, where the time and spatial functions are assumed to be of orthogonal collocation state at each element. The initial condition of each finite element can be obtained from the above solution. The heat fluxes on the vertical planes are calculated by the Galerkin FEM (finite element method). The CPFS system was modeled, simulated, and analyzed here. The simulation results were illustrated in three-dimensional graphics. Through simulation, it was shown clearly that the temperature distribution was influenced very much by the number, position, and temperature of the cable stream, and that dynamic heat transfer through the cable stream was one of the most dominant factors, and that the feature of heat conduction could be understood as an unsteady-state process.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.3B
• /
• pp.231-237
• /
• 2009

This study presents a numerical model to simulate density currents developing two dimensionally. The ${\kappa}-{\varepsilon}$ model is used for the turbulence closure. Elliptic flow equations are solved by the finite volume method. In order to investigate the applicability of the numerical model, discontinuous density currents are simulated numerically. The vortices due to the instability at the interface are simulated, showing a good agreement with the experimental visualizations in the literature. It is also investigated that the transition from slumping phase to inertial phase occurs when a bore generated at the end wall overtakes the front. However, the propagation of the density current is retarded compared with the experimental results. Two-dimensional modeling seems to have an effect on underestimating the front velocity of the density current.

• Journal of Environmental Impact Assessment
• /
• v.8 no.4
• /
• pp.47-57
• /
• 1999

This study is aimed at the development of comprehensive environmental management system, which can be operated on the basis of world wide web, as a topic of G7 project. Even though there should be lots of works remaining to achieve this goal, preliminary products can be summarized as follows : 1) integrated environmental information management system, 2) web based control engine, 3) surface water environment management system, 4) subsurface water environment management system, 5) sewer and waterworks management system. The core methodology of the engine is the generalized multidimensional finite element matrices to depict the terms in the analysis of various partial differential equations. Spatial information management system (ArcView) and Visual Basic were extensively employed to construct GUI oriented web based engine. The developed systems were composed of very intense computer codes due to the necessity of combinatory management of environmental problems. The web based engine could be served as a decision tool for the integrated management of environmental projects in Korea.

• Journal of Environmental Impact Assessment
• /
• v.10 no.4
• /
• pp.319-326
• /
• 2001

This study is aimed at the development of a comprehensive web-based partial differential equation solver (WPDES) using multidimensional finite element method, which can be operated on the basis of world wide web. Overall issues of engineering and environmental information management and facility control could be implemented using this solver. This paper describes the development technique of the model, which is first part on development of partial differential equation solver. Conventional commercial general solver of computational fluid dynamics problems were investigated. All the relevant environmental models were analyzed to develop integrated environmental management system using WPDES. The governing equations and the parameters of investigated models were analyzed and integrated. Several numerical modules were invented for each partial differential term in partial differential equation of many related modeling problems. Each module was coded in the fashion of object oriented method, and was combined independently for the overall governing equation. WPDES has unique characteristic, which can analyze the problem through the suitable combination of modules without development of additional models for each environment problem with different governing equation, main variables, and parameters.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.3
• /
• pp.453-462
• /
• 1994

Variation of seismic wave field in a multi-layered attenuating elastic half space is studied by the propagator matrix method and point source models of which fault-slip functions are defined as ramp functions. In this paper, the earth is modeled as being composed of horizontally stratified layers, with uniform material properties for each layer. The partial differential equations for the seismic motion in each layer are solved using a Fourier Hankel transform approach. Time histories and frequency contents of accelerations and displacements due to a vertical dip-slip and strike-slip point source located in the underlain half space are calculated at the layer interfaces using the developed programs and their characteristics are represented.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.1 no.3
• /
• pp.156-163
• /
• 1977

길이방향의 열전도효과가 비교적 큰유체- 고체 열교환계내에 일정한 조ㅍ은 온도를 가진 유체가 유입될 때 이에 따른 비정상 온도분포가 해석되었다. 이러한 계를 기술하는 연립편미분 방정식의 수치해가 유한차분법에 의해 구해졌으며, 열전도가 무시된 경우에 대한 해석적 해를 구하여 위의 결과를 검증하였다. 한편 이상화된 모델에 대한 실험을 행하고 이 결과를 앞서 구한 수치해와 비교 검토하였다. 여기서 구해진 온도분포를 이용하여 축열효율${\eta}_A$및 ${\eta}_B$를 ${\eta}_A$는 시간 t에서의 축열량과 정상상태에 도달했을 때의 축열량의 비로 , 그리고 ${\eta}_B$는 시간 t 동안의 축열된 양과 공급된 열량의 비고 정의하여 온도분포와 효율들에 애하여 시간변수 .tau., 위치변수 .xi.와 무차원 열확산율 .betha.가 미치는 영향을 조사하였다. 이에 따르면 온도는 .betha. 가 적어짐에 따라 점점 가파르게 분포되고 한편 ${\eta}_A$는 .tau.가 점점 커지거나, .xi.가 작아질수록 좋아지는 것을 알 수 있다. 또한 ${\eta}_B$는${\eta}_A$와는 정반대의 경향을 보여주어 축열계의 최적상태가 존재함을 알 수 있다.

• Journal of the Institute of Electronics Engineers of Korea CI
• /
• v.42 no.1
• /
• pp.79-87
• /
• 2005

This paper presents an active contour model for extracting accurate facial regions in complex environments. In the model, a contour is represented by a zero level set of level function φ, and evolved via level set partial differential equations. Then, unlike general active contours, skin color information that is represented by 2D Gaussian model is used for evolving and slopping a curve, which allows the proposed method to be robust to noise and varying pose. To assess the effectiveness of the proposed method it was tested with several natural scenes, and the results were compared with those of geodesic active contours. Experimental results demonstrate the superior performance of the proposed method.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2006.06a
• /
• pp.131-135
• /
• 2006

A seismic waveform inversion for prestack seismic data based on the Gauss-Newton method is presented. The Gauss-Newton method for seismic waveform inversion was proposed in the 80s but has rarely been studied. Extensive computational and memory requirements have been principal difficulties. To overcome this, we used different sizes of grids in the inversion stage from those of grids in the wave propagation simulation, temporal windowing of the simulation and approximation of virtual sources for calculating partial derivatives, and implemented this algorithm on parallel supercomputers. We show that the Gauss-Newton method has high resolving power and convergence rate, and demonstrate potential applications to real seismic data.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.2
• /
• pp.103-114
• /
• 1992

A study result of buckling reliability is presented for the axially compressed imperfect elastic cylinder. Multi-mode analysis program is developed from Karman-Donnell Equation for the calculation of the buckling load of the cylindrical she1l. Geometric intial imperfection is approximated by double Fourier series of which coefficients are assumed random variables with jointly normal distribution characteristics. Crude Monte Carlo simulation technique is used to calculate the probabilistic failure properties of several cases with various imperfection Conditions.