• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2000.11a
• /
• pp.20-20
• /
• 2000

2개의 경사 충격파와 하나의 수직 충격파로 초음속 유동을 압축하는 초음속 공기 흡입구의 수치적 연구를 수행하였다. 지배방정식으로는 Navier-Stokes방정식을 사용하였고 난류모델로는 SST 모델을 사용하였다. 지배방정식의 점성항 계산에는 중심차분법을 사용하였고 대류항 계산에는 풍상차분법인 Roe의 FDS기법을 MUSCL기법과 결합하여 이용하였다. 유한 체적법을 이용하여 차분된 방정식은 LU분할 기법을 이용한 완전 내재적 방법으로 2차 정확도 시간 적분으로 비정상 과정의 연구를 수행하였다. 흡입구 배압을 정해주어야 하는 어려움을 해결하기 위해 흡입구 후면에 노즐을 달고 노즐의 면적을 조절하여 배압이 형성되도록 하였다.(중략)

• Journal of the KSME
• /
• v.29 no.4
• /
• pp.403-413
• /
• 1989

유한요소법의 전산유체 역학분야에 대한 응용현황을 계산방법과 적용례를 중심으로 정리하였다. 유한요소법의 가장 큰 장점은 복잡한 유동영역을 해석하기 위한 불규칙 요소망(unstructured mesh)의 사용이라 볼 수 있으며 적응적 요소망을 이용하여 계산의 정확도를 높일 수 있는 것 또한 강점이라 할 수 있다. 다만 불규칙 요소망 사용으로 인해 수반되는 대수 방정식 계산시간 및 기억용량의 증가는 conjugate gradient 방법 등을 이용하여 반드시 해결되어야만 한다. 지금 까지 유한요소법을 이용한 계산방법을 개발해 오는 과정을 보면 유한차분법에서 오래 전에 개 발된 방법들을 도입한 경우가 많았으며 특히 난류 및 개발된 경우가 많으며 대부분의 경우 이 들을 그대로 도입, 이용하였다. 반대로 최근에 항공기 동체설계 분야를 중심으로 복잡한 형태의 유동영역을 해석이 요구되는 경우 유한차분법, 특히 유한체적법(finite volume method)에 삼각형 유한요소를 이용한 불규칙 요소망을 도입하여 성공적으로 이용하고 있다. 따라서 전산유체 역 학의 발전을 위하여 두 분야의 유기적인 협조가 필요하며 결과적으로 전산유체 역학기법이 완 전히 기계설계의 한 분야로 정립될 수 있도록 많은 노력이 필요하다고 본다.

• Proceedings of the Korea Contents Association Conference
• /
• 2012.05a
• /
• pp.143-144
• /
• 2012

전통적인 VOF 기법을 이용한 내부 조파 방법을 레블셋 기법에 적용하였다. 기하학적으로 유리한 유한요소법을 이용하여, Navier-Stokes 방정식의 공간차분에는 Characteristic Galerkin 기법을, 시간차분에는 Fractional Four-step 기법을 적용하였다. 중심(x=0)에서 전파하는 경우, 외부조파에 의한 영역내 재반사 문제가 해결되어 선형파를 의도한 바대로 잘 조파할 수 있었다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.181-184
• /
• 2012

기존에 제시된 Lin 과 Liu (1999)의 VOF 기법을 이용한 내부 조파 방법을 레블셋 기법에 적용하였다. 기하학적으로 유리한 유한요소법을 이용하여, Navier-Stokes 방정식의 공간차분에는 Characteristic Galerkin 기법을, 시간차분에는 Fractional Four-Step 기법을 적용하였다. 그림에 보인 바와 같이 중심(x=0)에서 전파하는 경우, 외부조파에 의한 영역내 재반사 문제가 해결되어 선형파를 의도한 바대로 잘 조파할 수 있었다.

• 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 주위의 재순환 영역에 대해서 살펴볼 수 있었다.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.11 no.1
• /
• pp.160-169
• /
• 1987

Numerical simulation of the flow of upper convected Maxwell fluid through planar 4:1 contraction has been performed using type dependent difference apprximation of vorticity equation. For creeping flow assumption, the numerical convergence has been achieved up to much higher values of elasticity parameter than those obtained by conventional finite difference method. For non-vanishing Reynolds number flow, it is shown that the corner vortices disappear, which is in good qualitative agreement with extant experimental results. In doing so, spatial distributions of stream function, vorticity and stresses are considered in relation to change of type of vorticity.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.3
• /
• pp.80-86
• /
• 1994

A numerical method for simulations of inviscid incompressible flow fields around a ship advancing on the free surface is developed. A body fitted coordinate system, generated by numerically solving elliptic type partial differential equations is used to conform the ship and free surface configurations. Three dimensional Euler equations transformed to the non-staggered body fitted coordinate system are discretised by finite difference method. Time and spatial derivatives are discretised by forward and centered differencings, respectively, and artificial dissipations are added to discretised convection terms for improvements of numerical stability. At each time steps, free surface elevations are recomputed to satisfy nonlinear free surface conditions. Poisson equations for pressure field are solved iteratively and the velocity field for next time step is extrapolated. To verify the developed numerical method, flow fields around a Wigley model are simulated(Fn=0.250-0.408) and compared with experimental data to show good agreements.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2009.05a
• /
• pp.815-821
• /
• 2009

하천의 2차원 흐름 및 하상변동, 오염확산 해석을 위한 유체의 수치해석법에는 유한요소법, 유한차분법, 유한차분법의 변형인 유한체적법, 경계적분법 등이 있으며, 국내의 경우 비구조적 요소망(unstructured mesh)을 이용하여 복잡한 형상을 표현하기가 상대적으로 용이한 유한요소법이 널리 사용되고 있다. 하천을 유한 요소화 하는 전처리 과정은 전체 해석 과정을 자동화 하는데 있어 필수적인 요소이며, 주로 삼각 요소망 또는 사각 요소망을 이용하여 해석을 수행하게 된다. 삼각 요소망의 경우 상대적으로 자동화하기 쉬운 반면 사각 요소망의 생성은 절점 생성 자체가 삼각 요소망 보다 더 많은 기하학적 제한 요소를 가지고 있기 때문에 상대적으로 완성도 높은 알고리즘을 구현하기가 어렵다 할 수 있다. 이에 따라 본 연구에서는 2차원 상에서 사각 요소망(quadrilateral elements)을 생성할 수 있는 Paving method를 중심으로 한 요소망 생성 알고리즘에 대해 고찰하고, 국내 최초의 범용 수치해석 모형인 RAMS(River Analysis and Modeling System)에 적용하였다. Paving method는 1990년에 Blacker and Stephenson에 의해 제안되었으며, Sandia National Laboratories에 의해 완성되었다. Paving Method는 advancing front style의 요소망을 생성하게 되고, 바깥쪽에서 안쪽으로 element layer를 생성하면서 채워나간다. 본 연구에서는 기존의 요소망 생성 프로세스에서 element 삽입 전의 검증 기능을 강화한 새로운 버전의 paving method를 적용하엿다.

• Publications of The Korean Astronomical Society
• /
• v.5 no.1
• /
• pp.40-64
• /
• 1990

복사전달식, 통계평형식 및 전하 입자 보존식을 동시에 만족시키는 비열평형상태의 대기에 대하여 그 방출원자스펙트럼을 수치계산하였다. 등온, 중압의 비열평형 태양홍염을 대상모델로 선정하였는데 여기에 적응한 제한조건중 복사전달에 관련된 편미분방정식은 3점 근사 차분법에 의해 정리하였고 홍염중심에 대칭성을 가정하여 경제조건을 부여하였으며 대기의 물리상태에 관련된 비선형 연립방정식의 해는 완전선형화 기법을 통한 퓨트리에 소거법을 이용해 구하였다.

• Journal of the Korea Convergence Society
• /
• v.12 no.2
• /
• pp.29-34
• /
• 2021

In this paper, a cylindrical photonic crystal waveguide with a low-index core is first proposed. The core can be filled with air, liquid, or arbitrary dielectric materials. Exact analyses for the electromagnetic field characteristics of guided modes, by using appropriate Bessel functions and applying the boundary conditions, are performed to find out the guiding characteristics of the proposed waveguide. For verification and usage in design and manufacturing process, the computer-calculation of the waveguide transmission characteristics is also performed by applying the rigorous full-vectorial finite difference method. Providing variations of the effective area for the fundamental mode of the designed waveguide with different numbers of cladding layers, ranging from 2.6056 ㎛2 to 5.9673 ㎛2 over the operation wavelength, generally as the core refractive index n1 is higher, the mode area becomes smaller and the result leads to more optimistic effect for nonlinear device applications.