• 한국전산유체공학회지
• /
• 제15권2호
• /
• pp.21-27
• /
• 2010

Immersed boundary lattice Boltzmann method (IBLBM) has been applied to simulate a turbulent flow over circular cylinder in a flow field effectively. Although IBLBM is very effective method to simulate the flow over a complex shape of obstacle in the flow field regardless of the constructed grids in the calculation domain, the results, however, become numerically unstable in high reynolds number flow. The most effective suggestion to archive the numerical stability in high Reynolds number flow is applying the multiple relaxation time (MRT) model instead of single relaxation time(SRT) model in the collision term of lattice Boltzmann equation. In the research MRT model for IBLBM was introduced and comparing the numerical results obtained by applying SRT and MRT. The hydraulic characteristic of cylinder in a flow field between two parallel plate at the range of $Re{\leqq}2000$represented and it is also compared the drag and lifting coefficients of the cylinder calculated by IBLBM with SRT and MRT model.

• Advances in aircraft and spacecraft science
• /
• 제3권4호
• /
• pp.367-378
• /
• 2016

Thrust Vectoring is a dynamic feature that offers many benefits in terms of maneuverability and control effectiveness. Thrust vectoring capabilities make the satisfaction of take-off and landing requirements easier. Moreover, it can be a valuable control effector at low dynamic pressures, where traditional aerodynamic controls are less effective. A numerical investigation of Fluidic Thrust Vectoring (FTV) is completed to evaluate the use of fluidic injection to manipulate flow separation and cause thrust vectoring of the primary jet thrust. The methodology presented is general and can be used to study different techniques of fluidic thrust vectoring like shock-vector control, sonic-plane skewing and counterflow methods. For validation purposes the method will focus on the dual-throat nozzle concept. Internal nozzle performances and thrust vector angles were computed for several range of nozzle pressure ratios and fluidic injection flow rate. The numerical results obtained are compared with the analogues experimental data reported in the scientific literature. The model is integrated using a finite volume discretization of the compressible URANS equations coupled with a Spalart-Allmaras turbulence model. Second order accuracy in space and time is achieved using an ENO scheme.

• 대한기계학회논문집A
• /
• 제28권6호
• /
• pp.774-782
• /
• 2004

This article presents a numerical and experimental investigation for the single-phase forced laminar convective heat transfer through arrays of micro-channels in micro heat exchangers to be used for cooling power-intensive semiconductor packages, especially the stacked multi-chip modules. In the numerical analysis, a parametric study was carried out for the parameters affecting the efficiency of heat transfer in the flow of coolants through parallel rectangular micro-channels. In the experimental study, the cooling performance of the micro heat exchanger was tested on prototypes of stacked multi-chip modules with difference channel dimensions. The simulation results and the experiment data were acceptably accordant within a wide range of design variations, suggesting the numerical procedure as a useful method for designing the cooling mechanism in stacked multi-chip packages and similar electronic applications.

• 대한기계학회논문집B
• /
• 제21권3호
• /
• pp.414-426
• /
• 1997

Numerical computations were conducted to characterize the three-dimensional laminar flow through an injector orifice having an inclined angle of 30 .deg.. For this study, the incompressible Navier-Stokes equations in generalized curvilinear coordinates, using a pseudocompressibility approach for continuity equation, were solved. The computations were performed using the finite difference implicit, approximately factored scheme of Beam and Warming and multi-block grids of complete continuity at block interfaces. The multi-block computations were validated for the steady state using direct comparison of multi-block solutions with equivalent single-block ones, including 2-D 180.deg. TAD and 3-D 90.deg. pipe bend. The comparisons between the numerical solutions and the flow field measurements for a tube with sudden contraction were presented in this work for solution validation. Computational results showed the nature of complex flow fields within the inclined injector orifice, including strong pressure-driven secondary flows in the cross stream induced by the effect of streamline curvature. In addition, asymmetric secondary flows were induced in the Reynolds number range above assumed laminar flow regime considered. However, turbulence calculations and grid dependency studies are needed for more accurate computations.

• 대한기계학회논문집
• /
• 제18권3호
• /
• pp.529-539
• /
• 1994

In laser heat treatment process of steels, the thin layer of substrate is rapidly heated to the austenitizing temperature and subsequently cooled at a very fast rate due to the self-quenching effect. Consequently, it is transformed to martensitic structure which has low magnetic permeability. This observation facilitates the use of a sensor measuring the change of electromagnetic field induced by the hardening layer. In this paper, the eddy-current electromagnetic field is analyzed by a finite element method. The purpose of this analysis is to investigate how the electrical impedance of the sensor's sensing coil varies with the change in permeability. To achieve this, a numerical model is formulated, taking into consideration the hardening depth, distance of the sensor from the hardened surface and the frequency driving the sensor. The results obtained by numerical simulation show that the eddy-current measurement method can feasibly be used to measure the changing hardening depth within the frequency range from 10 kHz to 50 kHz.

• 한국가시화정보학회지
• /
• 제14권3호
• /
• pp.32-37
• /
• 2016

In this study, the mixing process of two distinct flow is numerically investigated. Two flow with different physical properties (resin and hardener) are mixed through the opposing mixing jets. At a high pressure mixing process, the high speed flow is provided by two in-line nozzles. In the case of numerical modeling, Reynolds-Averaged Navier-Stokes Equations (RANS) is conducted to model the flow pattern inside the chamber. Additionally, SST k-omega turbulence model is selected to predict the kinetic energy of flow in impingement zone. The results show that mixing of two distinct flows would be efficient if the velocity of jet is high enough and nozzle diameter is a predominant parameter. Also, this velocity would create higher shear stress between two distinct flows which increases the mixing quality as well as strength of formed vortices. Eventually, the histogram of concentration fraction of resin is examined in order to show the quality of mixing and the range of concentration fractions in the output of chamber.

• 비파괴검사학회지
• /
• 제26권6호
• /
• pp.403-410
• /
• 2006

Eddy current testing (ECT) signals produced by a differential bobbin coil probe vary according to probe design parameters such as the number of turns, geometry and coil gap size. In the present study, the characteristics of a differential bobbin coil probe signals are investigated by numerical simulation in order to determine the optimum coil gap. For verification of numerical simulation accuracy, a specially designed bobbin probe of which the coil gap can be adjusted is fabricated and a series of experiments to acquire signals from two kinds of standard tubes with the variation in coil gap is performed. Then, the experimental signals are compared to the simulation results. Based on this investigation, a decision on the optimum range of coil gap is made. The theoretically predicted signals agree very well to the experimental signals. In fact, this excellent agreement demonstrates a high potential of the simulation as a design optimization tool for ECT bobbin probes.

• 한국전산구조공학회논문집
• /
• 제17권1호
• /
• pp.59-65
• /
• 2004

응력확대계수는 균열진전경로의 수치해석적 연구에서 널리 사용되고 있다. 그러나 많은 경우에서 균열선단주위 응력의 급수전개식 중 이어지는 항은 정량적으로 중요하다. 따라서 본 연구에서는 이의 항을 계산하기 위하여 등매개 2차특이요소를 이용한 유한요소해석을 수행하였다. 일례로 단축하중을 받는 무한 직방성체 내 경사균열에 대하여 균열요소크기와 균열경사각을 달리 주어 가며 해석을 수행하였으며, 수치해석결과는 이론해와 비교하여 잘 일치하고 있다.

• 조명전기설비학회논문지
• /
• 제28권2호
• /
• pp.76-83
• /
• 2014

This paper deals with the numerical method of calculating the frequency-dependent impedances of grounding electrodes. The proposed electromagnetic field approach is based on the solutions to Maxwell's equations obtained from the method of moment in the frequency domain. In order to evaluate the quality of the proposed simulation method, the frequency-dependent impedances of horizontally-buried ground electrodes were presented. The program for calculating the current distributions and impedances of grounding electrodes was implemented in MATLAB. The grounding impedances of two 10m and 50m long horizontal ground electrodes were measured and simulated in the frequency range from 100Hz to 10MHz for easy analysis and comparison. Also the simulated results were compared with those calculated from a sophisticated computer program CDEGS (HIFREQ module). As a result, the resultant results of frequency-dependent impedances obtained by using the numerical simulation method proposed in this work are in good agreement with experimental data. The validity of the approach techniques was confirmed.

• 한국산학기술학회논문지
• /
• 제4권3호
• /
• pp.274-278
• /
• 2003

터널의 3차원 해석을 통한 발파진동 평가시에 최대진동속도를 산정하기 위해 구조물 하부에서의 진동응답이 감소하는 터널직경의 3D위치까지의 경계를 설정하여 수치해석을 수행하였다. 구조물하부에서 터널발파시 진동성분이 장약으로부터 구형형상으로 전파되어나갈 때 구조물에서 발파막장의 거리가 터널직경의 약 1-2배일 때 최대진동속도를 나타내었으며, 수치해석시에는 속도성분영역을 확인하여 정확한 구형의 발파진동속도영역의 산정이 중요한 것으로 나타났다. 또한, 속도성분과 기존의 발파진동전파식을 비교한 결과 터널심도 15m, 25m일 경우에는 유사하게 나타났으나, 터널심도 35m일 경우에는 기존 발파진동전파식과는 상당한 차이를 보이므로 터널심도가 터널직경의 5배 이상일 경우에는 기존발파진동전파식에 의한 속도성분이 수치해석에 비해 과다하게 산정됨을 알 수 있었다.