• 한국음향학회지
• /
• 제34권1호
• /
• pp.36-45
• /
• 2015

유동광대역소음을 효율적으로 예측하기 위하여 통계적으로 난류를 재생하는 방법에 대한 많은 연구들이 최근에 진행되고 있다. 그 중에서도, FRPM(Fast Random Particle Mesh) 기법은 RANS(Reynolds-Averaged Navier-Stokes) 방정식 해석을 통해 도출된 정상상태 유동장의 난류 운동에너지와 소산 값을 이용하여 특정한 통계적 특성을 가지는 난류를 재생하는 기법으로서 유동광대역소음 문제 등에 성공적인 적용 예에 대해서 보고되고 있다. 하지만 기존의 FRPM 방법은 축류팬과 같이 축 대칭 특성을 갖는 기계의 경우 정상상태의 유동장을 기초로 광대역소음을 예측하는 문제에는 적용할 수 있으나, 원심팬과 같이 볼루트 영역으로 인하여 축 대칭이 성립되지 않는 기계류의 유동광대역소음에는 적용할 수 없다. 본 연구에서는 이러한 FRPM 기법을 확장하여, 원심팬에서 발생하는 광대역소음을 효율적으로 예측하기 위하여 비정상 RANS 방정식의 수치해와 연계하여 광대역소음원으로 고려되는 난류를 특정한 통계적 특성을 가지도록 재생할 수 있는 U-FRPM(Unsteady-FRPM) 기법을 제안하였다. 먼저 전산유체역학을 사용하여 RANS 방정식을 해석함으로써, 원심팬 주위의 비정상상태 유동장 정보를 도출하고, 음향상사법(Acoustic Analogy)을 기초로 도출된 유동소음원을 U-FRPM을 이용하여 모델링하였다. 모델링된 소음원은 경계요소법을 통해 구현되는 선형음향전파모델과 연계하여 수음점에서 광대역소음을 예측하는데 이용되었다. 예측된 결과와 실험결과의 비교를 통해 본 논문에서 제시한 방법의 유효성을 확인하였다.

• 한국항공우주학회지
• /
• 제50권7호
• /
• pp.445-454
• /
• 2022

항공기가 빙점 이하의 습도가 높은 구름대를 지날 때 액적이 항공기와 충돌하면 날개, 동체 등 항공기 구성품에 결빙이 발생한다. 특히 항공기의 날개에 결빙이 증식되면 공력 성능의 저하와 비행 안정성의 감소 등의 치명적인 안전 문제를 초래할 수 있다. 본 연구에서는 항공기 날개에 적용되는 고양력 장치인 다중 익형의 결빙 증식량이 최소가 되도록 형상 최적설계를 수행하였다. 3차원 Reynolds-Averaged Navier-Stokes 지배 방정식을 이용하여 공력해석을 수행하였고, 다물리 전산해석을 통해 결빙의 형상 및 증식량을 예측하였다. 최적설계의 목적함수는 결빙 증식량 최소화로 설정하였고, 설계변수는 Slat과 Flap의 전개 각도와 위치를 정의하는 형상 변수 6개를 선정하였다. 설계 과정에서 목적함수의 평가는 크리깅 근사모델을 사용하여 대체하였고 유전자 알고리즘을 적용하여 최적 형상을 도출하였다. 최적화를 수행한 결과, Slat과 Flap에 최적의 전개 각도와 위치를 적용하였을 때 결빙 증식량이 약 8% 감소하였다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1999년도 추계 학술대회논문집
• /
• pp.127-138
• /
• 1999

본 연구에서는 2 차원 캐비티 유동에서 경계조건 부여 방법의 차이에 따른 유체유동의 이력이 최종적인 정상상태에 미치는 효과를 수치실험을 통하여 관찰하였다. 그 결과 유동의 Reynolds 수가 작은 경우는 유동장이 가지는 이력이 시간의 경과와 함께 소멸하여 하나의 정상유동 상태를 가지는데 반하여, 큰 Reynolds 수에서는 유동장 이력이 소멸하지 않고 유지됨으로 인하여 다수의 정상유동 형태가 나타남이 확인이 되었다. 엄밀한 이론적 해석은 현재 이루워지지 않았으나 유체 운동량의 증가에 따라 Navier-Stokes 방정식의 비선형성이 영향을 미치고 있다고 여겨진다.

• 한국전산응용수학회:학술대회논문집
• /
• 한국전산응용수학회 2003년도 KSCAM 학술발표회 프로그램 및 초록집
• /
• pp.1-1
• /
• 2003

In this talk, a reduced-order modeling methodology based on centroidal Voronoi tessellations (CVT's)is introduced. CVT's are special Voronoi tessellations for which the generators of the Voronoi diagram are also the centers of mass (means) of the corresponding Voronoi cells. The discrete data sets, CVT's are closely related to the h-means clustering techniques. Even with the use of good mesh generators, discretization schemes, and solution algorithms, the computational simulation of complex, turbulent, or chaotic systems still remains a formidable endeavor. For example, typical finite element codes may require many thousands of degrees of freedom for the accurate simulation of fluid flows. The situation is even worse for optimization problems for which multiple solutions of the complex state system are usually required or in feedback control problems for which real-time solutions of the complex state system are needed. There hava been many studies devoted to the development, testing, and use of reduced-order models for complex systems such as unsteady fluid flows. The types of reduced-ordered models that we study are those attempt to determine accurate approximate solutions of a complex system using very few degrees of freedom. To do so, such models have to use basis functions that are in some way intimately connected to the problem being approximated. Once a very low-dimensional reduced basis has been determined, one can employ it to solve the complex system by applying, e.g., a Galerkin method. In general, reduced bases are globally supported so that the discrete systems are dense; however, if the reduced basis is of very low dimension, one does not care about the lack of sparsity in the discrete system. A discussion of reduced-ordering modeling for complex systems such as fluid flows is given to provide a context for the application of reduced-order bases. Then, detailed descriptions of CVT-based reduced-order bases and how they can be constructed of complex systems are given. Subsequently, some concrete incompressible flow examples are used to illustrate the construction and use of CVT-based reduced-order bases. The CVT-based reduced-order modeling methodology is shown to be effective for these examples and is also shown to be inexpensive to apply compared to other reduced-order methods.

• 한국전산유체공학회지
• /
• 제21권2호
• /
• pp.112-119
• /
• 2016

NEXTSat-1 is the next-generation small-size artificial satellite system planed by the Satellite Technology Research Center(SatTReC) in Korea Advanced Institute of Science and Technology(KAIST). For the control of attitude and transition of the orbit, the system has adopted a RHM(Resisto-jet Head Module), which has a very simple geometry with a reasonable efficiency. An axisymmetric model is devised with two coil-resistance heaters using xenon(Xe) gas, and the minimum required specific impulse is 60 seconds under the thrust more than 30 milli-Newton. To design the module, seven basic parameters should be decided: the nozzle shape, the power distribution of heater, the pressure drop of filter, the diameter of nozzle throat, the slant length and the angle of nozzle, and the size of reservoir, etc. After quasi one-dimensional analysis, a theoretical value of specific impulse is calculated, and the optima of parameters are found out from the baseline with a series of multi-physical numerical simulations based on the compressible Navier-Stokes equations for gas and the heat conduction energy equation for solid. A commercial code, COMSOL Multiphysics is used for the computation with a FEM (finite element method) based numerical scheme. The final values of design parameters indicate 5.8% better performance than those of baseline design after the verification with all the tuned parameters. The present method should be effective to reduce the time cost of trial and error in the development of RHM, the thruster of NEXTSat-1.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2004년도 춘계 학술대회논문집
• /
• pp.146-153
• /
• 2004

The Numerical study has been carried out to investigate the effects of chemical reaction and thermal radiation on the rocket plume flow-field at various altitudes. The theoretical formulation is based on the Navier-Stokes equations for compressible flows along with the infinitely fast chemistry and thermal radiation. The governing equations were solved by a finite volume fully-implicit TVD(Total Variation Diminishing) code which uses Roe's approximate Riemann solver and MUSCL(Monotone Upstream-centered Schemes for Conservation Laws) scheme. LU-SGS (Lower Upper Symmetric Gauss Seidel) method is used for the implicit solution strategy. An equilibrium chemistry module for hydrocarbon mixture with detailed thermo-chemical properties and a thermal radiation module for optically thin media were incorporated with the fluid dynamics code. In this study, kerosene-fueled rocket was assumed operating at O/F ratio of 2.34 with a nozzle expansion ratio of 6.14. Flight conditions considered were Mach number zero at ground level, Mach number 1.16 at altitude 5.06km and Mach number 2.9 at altitude 17.34km. Numerical results gave the understandings on the detailed plume structures at different altitude conditions. The diffusive effect of the thermal radiation on temperature field and the effect of chemical recombination during the expansion process could be also understood. By comparing the results from frozen flow and infinitely fast chemistry assumptions, the excess temperature of the exhaust gas resulting from the chemical recombination seems to be significant and cannot be neglected in the view point of performance, thermal protection and flow physics.

• 한국전산유체공학회지
• /
• 제16권4호
• /
• pp.39-46
• /
• 2011

A numerical study of a turbulent natural convection in an enclosure with the lattice Boltzmann method (LBM) is presented. The primary emphasis of the present study is placed on investigation of accuracy and numerical stability of the LBM for the turbulent natural convection flow. A HYBRID method in which the thermal equation is solved by the conventional Reynolds averaged Navier-Stokes equation method while the conservation of mass and momentum equations are resolved by the LBM is employed in the present study. The elliptic-relaxation model is employed for the turbulence model and the turbulent heat fluxes are treated by the algebraic flux model. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with the deferred correction way to ensure accuracy and stability of solutions. The present LBM is applied to the prediction of a turbulent natural convection in a rectangular cavity and the computed results are compared with the experimental data commonly used for the validation of turbulence models and those by the conventional finite-volume method. It is shown that the LBM with the present HYBRID thermal model predicts the mean velocity components and turbulent quantities which are as good as those by the conventional finite-volume method. It is also found that the accuracy and stability of the solution is significantly affected by the treatment of the convection term, especially near the wall.

• Ocean Systems Engineering
• /
• 제3권3호
• /
• pp.219-236
• /
• 2013

The accurate prediction of motion in waves of a marine vehicle is essential to assess the maximum sea state vs. operational requirements. This is particularly true for small crafts, such as Autonomous Surface Vessels (ASV). Two different numerical methods to predict motions of a SWATH-ASV are considered: an inviscid strip theory initially developed at MIT for catamarans and then adapted for SWATHs and new a hybrid strip theory, based on the numerical solution of the radiation forces by an unsteady viscous, non-linear free surface flow solver. Motion predictions obtained by the viscous flow method are critically discussed against those obtained by potential flow strip theory. Effects of viscosity are analyzed by comparison of sectional added mass and damping calculated at different frequencies and for different sections, RAOs and motions response in irregular waves at zero speed. Some relevant conclusions can be drawn from this study: influence of viscosity is definitely non negligible for SWATH vessels like the one presented: amplitude of the pitch and heave motions predicted at the resonance frequency differ of 20% respectively and 50%; in this respect, the hybrid method with fully non-linear, viscous free surface calculation of the radiation forces turns out to be a very valuable tool to improve the accuracy of traditional strip theories, without the burden of long computational times requested by fully viscous time domain three dimensional simulations.

• 대한조선학회논문집
• /
• 제54권6호
• /
• pp.470-475
• /
• 2017

Since the intake air of gas turbine engine of marine purpose contains water particles, inertial separator for separating the air and water particles are provided. Saw type and wave type separator are now used to separate inflow water particle from the gas. In this paper, the design parameters of saw type separator are studied by numerical simulations. Using the commercial CFD program, Star-CCM+, Lagrangian-Eulerian method was used to perform the analysis of two phase flow of the mist in the air. This method solves Reynolds-Averaged Navier-Stokes equations in Eulerian framework for the continuous phase, while solves equation of motion for individual particles in Lagrangian framework. Lagrangian multiphase method was applied to monitor the particles of different sizes and shapes and to verify collision between particles by chasing particles. Water particles were injected through injectors located at the inlet of the separator and escape mode was used which assumes that the particles attached on the surface of inertial separator were removed from the simulation, effectively escaping the solution domain. Through the numerical computations with the inlet condition of constant water particle size in the wetness fraction of 85%, efficiency of eliminating the water particle and the pressure drop between the inlet and outlet were examined.

• Coupled systems mechanics
• /
• 제6권3호
• /
• pp.287-316
• /
• 2017

This paper studies the particularities of the forced vibration of the hydro-elastic system consisting of a moving elastic plate, compressible viscous fluid and rigid wall. This study is made by employing the discrete-analytical solution method proposed in the paper by the authors (Akbarov and Panakhli (2015)). It is assumed that in the initial state the fluid flow is caused by the axial movement of the plate and the additional lineally-located time-harmonic forces act on the plate and these forces cause additional flow field in the fluid and a stress-strain state in the plate. The stress-strain state in the plate is described by utilizing the exact equations and relations of the linear elastodynamics. However, the additional fluid flow field is described with linearized Navier-Stokes equations for a compressible viscous fluid. Numerical results related to the influence of the problem parameters on the frequency response of the normal stress acting on the plate fluid interface plane and fluid flow velocity on this plane are presented and discussed. In this discussion, attention is focused on the influence of the initial plate axial moving velocity on these responses. At the same, it is established that as a result of the plate moving a resonance type of phenomenon can take place under forced vibration of the system. Moreover, numerical results regarding the influence of the fluid compressibility on these responses are also presented and discussed.