• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 1997.11a
• /
• pp.17-17
• /
• 1997

극초음속 여객기와 군사용 항공기에 대한 수요가 증가함에 따라서 새로운 개념의 다양한 추진기관이 연구가 진행되고 개발되어 왔다. 초음속 항공기의 속도 영역은 마하 10-20 정도가 되는데 이 속도 한계를 극복하기 위하여 초음속 연소 램제트 엔진(SCRamjet; Supersonic Combustion Ramjet)이 제안되었다. 스크램 제트를 개발하기 위해서는 연료와 산화제의 혼합 효율 문제, 화염의 안정화 문제, 벽면의 냉각에 관한 문제 등 몇 가지 기본적인 문제들을 해결해야 한다. Univ of Michigan에서 실험한 연소기를 모델로 본 연구에서는 연료와 공기의 혼합에 관한 수치 연구를 수행하였다. 다원 혼합기체에 관한 축대칭 Navier-Stokes 방정식을 지배 방정식을 이용하였고 비평형 화학반응식을 고려하였다. 공간 차분에는 유한 체적법을 이용하였다. 대류 플럭스 항은 Roe의 Upwind FDS 기법을 사용하여 차분하였고 점성항에는 중심 차분법을 이용하였다. 시간 적분법으로는 근사 자코비안과 LU분할 기법을 이용한 완전 내재적 방법이 쓰였다. 난류 모델로는 Mentor에 의해 제안된 2 방정식 k-$\varepsilon$/k-$\omega$ 혼합모델을 사용하였다. 유동장이 실험에서의 찍은 사진과 유사한 모습의 충격파 간섭을 수치 모사하였고 수소가 확산되는 모습과 함께 노즐 lip 주위의 재순환 영역에 대해서 살펴볼 수 있었다.

• Proceedings of the KIEE Conference
• /
• 2003.07d
• /
• pp.2221-2223
• /
• 2003

다양한 공정에 범용 적으로 적용할 수 있는 PID 제어기 동조 방법이 많이 연구되고 있다. 한 가지 방법으로서 실제공정에 대한 축소모델을 구하고 축소모델로부터 파라메터를 동조하는 방법이 사용되고 있다. 본 논문에서는 시간응답으로부터 최소 자승법을 이용하여 고차 시스템을 2차계로 축소하는 방법을 제안하였다. 2차계 시스템을 차분 방정식으로 표현하고, 시간응답의 입출력 데이터로부터 최소 자승법을 이용하여 차분 방정식의 계수를 구하였다. 차분 방정식을 이산시간상태방정식으로 변환, 이산시간 상태방정식을 연속시간 상태방정식으로 변환하는 과정을 유도하여, 차분 방정식의 계수로부터 연속시간 상태방정식을 구하는 방법을 제시하였다. 시뮬레이션을 통하여 축소모델의 정확성 및 실제공정과의 특성을 비교하였다.

• Journal of the Korean Society of Safety
• /
• v.10 no.2
• /
• pp.129-133
• /
• 1995

A new approach to the analysis of thin. rectangular window giass glass supported on flexible gaskets. and subjected to uniform lateral pressures was evolved. Based on the Von Karman theory of plates and using the finite difference method. a computer program which determines the deflections and stresses in simply supported thin glass plates was developed.

• Proceedings of the Korean Nuclear Society Conference
• /
• 2002.10a
• /
• pp.51.1-51
• /
• 2002

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.6
• /
• pp.583-590
• /
• 2016

This paper presents a dynamic crack propagation algorithm with Rayleigh damping effect based on the MLS(Moving Least Squares) Difference Method. Dynamic equilibrium equation and constitutive equation are derived by considering Rayliegh damping and governing equations are discretized by the MLS derivative approximation; the proportional damping, which has not been properly treated in the conventional strong formulations, was implemented in both the equilibrium equation and constitutive equation. Dynamic equilibrium equation including time relevant terms is integrated by the Central Difference Method and the discrete equations are simplified by lagging the velocity one step behind. A geometrical feature of crack is modeled by imposing the traction-free condition onto the nodes placed at crack surfaces and the effect of movement and addition of the nodes at every time step due to crack growth is appropriately reflected on the construction of total system. The robustness of the proposed numerical algorithm was proved by simulating single and multiple crack growth problems and the effect of proportional damping on the dynamic crack propagation analysis was effectively demonstrated.

• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.8-18
• /
• 1999

In devising a numerical approximation for the convective spatial transport of a fluid mechanical quantity, it is noted that the convective motion of a scalar quantity occurs in one-way, or from upstream to downstream. This consideration leads to a new scheme termed a pure upwind difference scheme (PUDS) in which an estimated value for a fluid mechanical quantity at a control surface is not influenced from downstream values. The formal accuracy of the proposed scheme is third order accurate. Two typical benchmark problems of a wall-driven fluid flow in a square cavity and a buoyancy-driven natural convection in a tall cavity are computed to evaluate performance of the proposed method. for comparison, the widely used simple upwind scheme, power-law scheme, and QUICK methods are also considered. Computation results are encouraging: the proposed PUDS sensitized to the convection direction produces the least numerical diffusion among tested convection schemes, and, notable improvements in representing recirculation of fluid stream and spatial change of a scalar. Although the formal accuracy of PUDS and QUICK are the same, the accuracy difference of approximately a single order is observed from the revealed results.

• Tunnel and Underground Space
• /
• v.30 no.5
• /
• pp.484-495
• /
• 2020

To numerically simulate the advance of EPB TBM, various type of numerical analysis methods have been adopted including discrete element method (DEM), finite element method (FEM), and finite difference method (FDM). In this paper, an EPB TBM driving model was proposed by using coupled DEM-FDM. In the numerical model, DEM was applied in the TBM excavation area, and contact properties of particles were calibrated by a series of triaxial tests. Since the ground around the excavation area was coupled with FDM, the horizontal stress considering the coefficient of earth pressure at rest could be applied. Also, the number of required particles was reduced and the efficiency of the analysis was increased. The proposed model can control the advance rate and rotational speed of the cutter head and screw conveyor, and derive the torque, thrust force, chamber pressure, and discharging during TBM tunnelling.

• Journal of the Korea institute for structural maintenance and inspection
• /
• v.9 no.4
• /
• pp.159-168
• /
• 2005

This study carries out a finite difference nonlinear analysis of anisotropic advanced composite plate structures with various layer sequences. In the numerical analysis of various mechanical problems involving complex partial differential equations, the finite difference method (FDM) developed in this study has an advantage over the finite element method in its ability to avoid mesh generation and numerical integration. Many studies in FDM have been made on clamped or simple boundary conditions using merely an energy approach. These approaches cannot be satisfied, however, with pivotal points along the free boundary. Therefore, this study addresses the nonlinear problem of anisotropic plates by adopting a finite difference modeling elimination of pivotal difference points in the case of a free boundary condition. Complex nonlinear behaviors of composite plate structures for various parameters, especially for layer sequences, are analyzed using the proposed approach.

• Journal of Energy Engineering
• /
• v.17 no.4
• /
• pp.233-240
• /
• 2008

Conventionally, the second-order self-adjoint neutron transport equations have been studied using the even parity and the odd parity equations. Recently, however, the SAAF(self-adjoint angular flux) form of neutron transport equation has been introduced as a new option for the second-order self-adjoint equations. In this paper we validated the SAAF equation mathematically and clarified how it relates with the existing even and odd parity equations. We also developed a second-order SAAF differencing formula including DSA(diffusion synthetic acceleration) from the first-order difference equations. Numerical result is attached to show that the proposed methods increases accuracy with effective computational effort.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.178-178
• /
• 2021

최근 들어 기후변화로 인한 극심한 가뭄 피해가 한반도에 발생하고 있다. 가뭄 상황에 대비하여 댐을 안정적으로 운영하기 위해서는 갈수빈도 유입량에 대한 분석이 필수적이다. 갈수빈도해석의 경우, 홍수빈도해석과 유사하게 확률밀도함수의 극값에 대한 확률값을 산정하며, 확률 분포형의 역함수에 비초과확률을 대입하여 산정한다. 그러나 홍수와 달리 가뭄은 지속기간이 긴 특성 탓에 자기상관을 고려해야하며, 댐 및 저수지 등 대규모 시설물의 경우 일반적인 하천과 달리 저류효과로 인해 누적 유량에 대한 고려가 필요하다. 이에 K-water는 자체 제작한 누가차분법 및 Disaggregation 두 가지 방법을 채택하여 실무에서 사용해왔다. 그러나 누가차분법을 사용할 경우, 빈도유입량이 지나치게 크게 산정되는 문제가 있으며, Disaggregation 방법을 사용하는 경우, 특정 빈도 이상의 극한가뭄에서 유입량의 차이가 유의미하지 않아 산정된 빈도유입량과 최근 발생한 극심한 가뭄의 실측유입량간 큰 차이가 발생하고 있다. 따라서 본 연구에서는 자기상관을 고려한 선형회귀모형에 근거하여 빈도유입량을 배분하는 방법을 제안한다. 또한, 앞서 서술한 네 가지 빈도유입량 방법(월빈도분석, 누가차분법, K-water Disaggregation, 자기상관 선형회귀모형)에 대한 수식적 비교를 수행하며, 국내 댐 유역에 적용 및 평가를 통해 자료 특성에 따른 적절한 빈도유입량 산정방식에 대한 기준을 제안한다. 본 연구를 통해 가뭄특성을 고려한 합리적인 댐 유입량을 산정함으로써 보다 유연한 수자원시설물의 가뭄대응이 이루어질 것으로 기대된다.