• Proceeding of EDISON Challenge
• /
• 2013.04a
• /
• pp.354-359
• /
• 2013

일반적으로 건축물의 설계시 풍동 실험을 통한 풍환경의 평가를 수행하고 있으며, 이는 환경 영향 평가법에서 정한 건축 사업 시행 시 수반되어야 할 자연환경, 생활환경 그리고 사회경제환경의 영향 평가의 일환으로 실시되고 있다. 그러나, 풍동 실험의 경우 여러 가지 현실적 제약조건으로 설계와 실험의 피드백 (Feedback)이 원활하지 못하며, 특히 대상 건축물이 공장과 같이 대기 오염원이 되는 경우 실험은 더욱 어려운 형편이다. 이에 대한 보완책으로 전산 유체 역학을 이용한 건축물의 풍압 해석에 의한 풍하중 추정이나 인접 지형-지물의 영향을 고려한 건축물 주위의 풍환경 평가가 있다. 전산 모사에 의해 풍동 실험의 미비점을 보완하고, 보다 상세한 정보를 확보함으로써 건축물의 구조적 안전성의 증대와 환경 피해 감소를 기할 수 있다. 그러나 복잡한 지형-지물이나 건축물 주위의 풍환경에 대한 전산 모사는 주로 두 가지의 기술적 어려움을 수반하게 되다. 그 중 하나는 고정 경계면을 이루는 형상의 복잡성으로 인해 기존에 많이 이용하고 있는 Body-fitted 격자계를 이용하는 경우, 격자 생성 과정이 매우 복잡하고 어려울 뿐 만 아니라 생성된 격자가 주로 비정렬 (unstructured) 특성을 갖게 되어 수치해석 과정의 효율을 저하시키는 요인이 되며, 격자의 형상도 수치해석의 수렴성을 저하시키는 예가 많다. 다른 어려움으로 풍환경은 전형적인 난류 유동장으로서 난류의 전산 해석은 아직도 해결하지 못한 부분이 많다는 점이다. 이에 본 논문에서는 복잡한 지형-지물이나 건축물의 풍하중과 풍환경의 전산 모사 기술 확보를 위하여 수행중인 연구의 일환으로 물체 형상의 기하학적 복잡성의 극복을 위한 가상경계법 (Immersed Boundary Method)과 난류 유동장의 물리적 엄밀성을 높이기 위한 다와동 모사 (Large Eddy Simulation)을 이용한 물체 형상과 무관한 유동장 해석 기술 개발에 대하여 다루고자 한다. 먼저 최근에 유동 해석에 이용되는 방법인 가상경계법(IBM)은 물체를 포함한 전체 전산 영역을 직교 좌표계에 의해 이산화하고, 유동장내 존재하는 물체의 표면에서의 점착 조건을 만족시키기 위하여 지배 방정식에 적절한 외력을 추가로 고려하는 방법이다. 본 연구에서는 가상경계법을 이용하여 경계층에 위치한 건물 형상의 각진 물체 주위 사이에 형성되는 공동 내부의 비정상 유속 및 압력에 대한 전산 해석을 수행하고, 풍상측 전면에 형성되는 경계층에 의한 영향을 분석하였다.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• v.y2005m4
• /
• pp.423-426
• /
• 2005

Numerical analysis of a gaseous jet injected into a closed nozzle was carried out considering the valve operation to simulate the interior pressure development of closed nozzle. Pressure data was obtained through valve operation and we compared numerical analysis with experiment data. Commercial code was used for the present calculation. Numerical results of a pressure development since its start showed good agreement for experiment data of the nozzle operation.

• The Korean Journal of Rheology
• /
• v.5 no.1
• /
• pp.34-48
• /
• 1993

최근 CAD/CAM의 발전과 더불어 사출성형공정은 여러분야에 폭넓게 응용되고 있 다. 사출성형공정은 크게 충전과정(filling stage), 냉각과정(cooling stage), 보압과정(packing stage)로 나누어 지는데 이중 충전과정은냉각과정과 보압과정에서 나타날 물리적인 현상과 최종 성형품의 기계적 성질에 중요한 영향을 끼치게 된다. 충전과정의 수치 해석 방법은 대 표적으로 control volume method, branching flow method, transient moving boun-dary method로 구분된다. 본 연구에서는 격자의 형태를 양호하게 형성시키고 유동선단의 형태를 개선하기위한 기법인 Spline 곡선을 이용한 선단격자 재구성(frontal remeshing using spline curve)과 수치해석에 소요되는 시간을 줄이기 위하여 벽면경계조건 처리를 위한 선단 유동 장생성(frontal flow field construction for wall boun-dary condition treatment)기법을 개발 하고 transient moving voundary method에 적용시켜 원형 평판과 인장 및 굽힘시편 그리고 두께가 변하는 사각 형상을 가진 캐비터에서의 충전과정을 수치해석하였다. 그결과 압력 분 포, 온도분포, 속도장, 유동선단의 진전형태 등이 기존에 제출된 해석결과와 비교하여 볼 때 만족스러운 수치해석결과를 보였다.

• The KSFM Journal of Fluid Machinery
• /
• v.2 no.1 s.2
• /
• pp.64-72
• /
• 1999

CSCM upwind flux difference splitting compressible Navier-Stokes method has been used to predict the transonic flows in a centrifugal compressor diffuser. The modified cyclic. TDMA and the mass flux boundary conditions were used as boundary conditions of the diffuser analysis. Broad flow separation on the suction surface near the hub and shroud was observed from the results of the mass flow rates 5.8, 6.0 and 6.2kg/s at 27000 rpm. The three-dimensional flow analysis predicted successfully that the static pressure increased and the total pressure decreased through the flow passage of the channel diffuser when compared to two-dimensional analysis due to the strong effect of the three-dimensional flow. The mass averaged loss coefficients and pressure coefficients were also studied.

• Journal of Ocean Engineering and Technology
• /
• v.26 no.5
• /
• pp.55-62
• /
• 2012

In this paper, the rigid lid boundary condition is applied to simulate the influence of floating structures such as ships or pontoons, and the pressure term in both the momentum equations and continuity equation are modified. The pressure of a floating structure under the free surface is dependent on the draft of the structure, generally called a ship. If the free surface is covered by a floating structure, the free surface cannot move freely. The water level should be fixed, using a rigid lid boundary condition. This boundary condition is implemented by reducing the storage area of the grid cell with a factor between zero and one. The numerical model developed by Hong (2009) is verified through a comparison with experimental results, and the influence of the reduction factor is investigated using the verified numerical model.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.15 no.5
• /
• pp.66-71
• /
• 2011

The spray distribution at high pressure condition was measured by the Tomography method. The constructed spray distribution was compared with the images by Indirect Photography method so that the spray size confirmed and took the boundary of the distribution. It confirmed that the Tomography applies to construct the distribution at high pressure.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.112-122
• /
• 1993

A potential-based panel method is presented for the analysis of a partially or supercavitating two-dimensional hydrofoil at a finite submergence beneath a free surface, treating without approximation the effects of the finite Froude number and the hydrostatic pressure. Free surface sources and normal dipoles are distributed on the foil and cavity surfaces, their strength being determined by satisfying the kinematic and dynamic boundary conditions on the foil-cavity boundary. The cavity surface is determined iteratively as a part of the solution. Numerical results show that the wave profile is altered significantly due to the presence of the cavity. The buoyancy effect due to the hydrostatic pressure, which has usually been neglected in most of the cavitating flow analysis, is found playing an important role, especially for the supercavitating hydrofoil; the gravity field increases the cavity size in shallow submergence, but decreases it when deeply submerged, while the lift reduces at all submergence depth.

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

고체 추진제를 사용하는 추진 시스템을 개발하는데 가장 커다란 문제로 인식되고 있는 것은 추진제의 연소 특성을 이해하는 일이다. 그 중에서도 연소실의 압력 진동과 추진제 벽면으로 흡수되는 복사 열전달에 의한 연소율(burning rate)의 변화로 인하여 발생하는 연소 불안정에 대한 이해는 아직도 완전히 규명되지 않고 있다. 고체 추진제의 연소 불안정에 대한 이론적 해석은 준-정상 1차원 해석(Quasi-Steady Homogeneous One-Dimension) 방법에 의하여 단순화된 지배방정식을 해석하는 것이 일반적으로 잘 알려져 있는 방법이다. 이 가정은 고체 추진제가 연수되는 영역을 두께가 매우 얇은 영역의 표면반응영역(surface reaction layer)과 화학반응이 없는 응축상태영역(condensed phase zone) 그리고 기체상태의 연료와 화염이 존재하는 기체상태영역(gas phase zone) 등의 3영역으로 구분하며, 기체상태영역에서 발생하는 교란에 대한 응축상태영역의 반응시간 크기(response time scale)가 매우 크기 때문에 응축상태영역의 반응은 준 정상적으로 일어난다고 가정하는 것이다.그러나, 연소실의 온도가 $3000^{\circ}K$ 정도의 높은 온도이어서 복사 열전달에 의한 고체 추진제의 가열이 중요한 열전달 방법으로 작용하게 되므로 이를 무시한 이론적 해석은 물리적인 중요성이 약하여질 수밖에 없다. 본 연구에서는 기체영역으로부터 전달되는 복사 열전달은 투명(transparent)한 표면반응영역을 통과하여 응축상태영역에서 모두 흡수되며 추진제 표면에서의 복사열방출(emission)을 고려하였다. 또한 연소불안정 현상을 해석하기 위하여 표면반응영역에서의 경계조건은 선형교란량으로 대치하는 Zn(Zeldovich-Novozhilov) 방법을 사용하였다. 이 방법은 기체상태영역에 대한 구체적인 해석없이도 연소불안정 현상을 해석할 수 있는 장점이 잇다. 즉 응축상태영역에서의 연소율과 표면온도는 각각 기체영역으로부터 전달되는 온도구배와 연소압력, 그리고 복사 열전달의 함수관계이므로 선형교란에 의한 추진제표면에서의 교란경계조건을 얻을 수 잇으며, 응축영역의 교란지배방정식과 함께 사용하여 압력교란과 복사 열전달의 교란에 대한 연소율의 교란 증감 여부를 판단하여 연소 불안정 현상을 해석할 수 있다.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2017.05a
• /
• pp.881-888
• /
• 2017

In general, the choking phenomenon occurs by flow acceleration for a turbine at high pressure ratio condition. In choking condition, total pressure ratio increases without mass flow rate variation. It is hard to predict choking characteristics by using conventional meanline analysis which used mass flow inlet boundary condition. In the present study, the algorithm for predicting choking point is developed to solve the problem. Moreover, performance estimation algorithm after choking is presented by reflecting the flow behaviour of flow expansion at choked nozzle or rotor. The analysis results are compared with 3D CFD analysis and experimental data to validate present method.

• Journal of the Society of Naval Architects of Korea
• /
• v.32 no.3
• /
• pp.62-71
• /
• 1995

In the numerical analysis of incompressible unsteady Navier-stokes equation, large time is required for solving the pressure Poisson equation of the elliptic type at each time step. In this paper, a numerical analysis by the direct method is carried out to solve the pressure Poisson equation and the computing time is analyzed as mesh size increases. The pressure Poisson equation can be transformed to the boundary value problem by the Green theorem. The computing time for the convolution type of the domain integral can be reduced by using F.F.T. and the computing time in the direct method depends entirely on obtaining the solution of the boundary value problem. The numerical analysis on the known solutions is carried out and compared for the verification of the direct method. And the numerical analysis on the body boundary and domain decomposition problem are carried out with the computing time less than O($n^{3}$) in the (n.n) mesh.