This is a short review of how CFD contributed to hypersonic flights in the past 50 years. Two unexpected phenomena that occurred in the entry flights of the Apollo and Space Shuttle made us aware of the impact of the high temperature real-gas effects on hypersonic flights: pitching moment anomaly of up to 4 degrees, and radiation overshoot behind a shock wave. The so-called two-temperature nonequilibrium model was introduced to explain these phenomena. CFD techniques were developed to accommodate the two-temperature model. Presently, CFD can predict trim angle of attack to an accuracy of about 1 degree. A concerted effort was made to numerically reproduce the experimentally measured flow-field over a double-cone. As yet, perfect agreement between the experimental data and computation is not achieved. Scramjet technology development is disappointingly slow. The phenomenon of ablation during planetary entries is not yet predicted satisfactorily. In the future, one expects to see more research carried out on planetary entries and space tourism.

This paper presents the state of the art of computational structures technology(CST) and comparison of two computational mechanics - CST and CFD, to the CFD engineers. Classical mechanics is based on the five classical axioms which describe the motion and behaviors of the continuum materials like solid structures and fluids. Computational structures technology uses the finite element method to solve the governing equation, whereas finite volume method is generally used in CFD. A few famous commercial structural analysis programs and DIAMOND/IPSAP will be introduced. DIAMOND/IPSAP is the efficient parallel structural analysis package developed by our research team. DIAMOND/IPSAP shows the better performance than the commercial structural analysis software not only in the parallel computing environments but also in a single computer.

Dynamic Stall is a flow phenomenon which occurs on the retreating side of helicopter rotor blades during forward flight. It also occurs on blades of stall regulated wind turbines under yawing conditions as well as during gust loads. Time scales occurring during this process are comparable on both helicopter and wind turbine blades. Dynamic Stall limits the speed of the helicopter and its manoeuvrability and limits the amount of power production of wind turbines. Extensive numerical as well as experimental investigations have been carried out recently to get detailed insight into the very complex flow structures of the Dynamic Stall process. Numerical codes have to be based on the full equations, i.e. the Navier-Stokes equations to cover the scope of the problems involved: Time dependent flow, unsteady flow separation, vortex development and shedding, compressibility effects, turbulence, transition and 3D-effects, etc. have to be taken into account. In addition to the numerical treatment of the Dynamic Stall problem suitable wind tunnel experiments are inevitable. Comparisons of experimental data with calculated results show us the state of the art and validity of the CFD-codes and the necessity to further improve calculation procedures. In the present paper the phenomenon of Dynamic Stall will be discussed first. This discussion is followed by comparisons of some recently obtained experimental and numerical results for an oscillating helicopter airfoil under Dynamic Stall conditions. From the knowledge base of the Dynamic Stall Problems, the next step can be envisaged: to control Dynamic Stall. The present discussion will address two different Dynamic Stall control methodologies: the Nose-Droop concept and the application of Leading Edge Vortex Generators (LEVoG's) as examples of active and passive control devices. It will be shown that experimental results are available but CFD-data are only of limited comparison. A lot of future work has to be done in CFD-code development to fill this gap. Here mainly 3D-effects as well as improvements of both turbulence and transition modelling are of major concern.

Ice accretion is one of the potential hazards in airplane flight, adversely affecting aircraft aerodynamic. There are two distinct icing analysis that can be simulated. One is predicting the effect of ice on the aerodynamic performance of airfoils when ice geometry is known. The other is simulating ice accretion. This work presents the method of icing accretion analysis. This work presents an Eulerian approach to calculate the droplet collection efficiency on the 2D airfoil. The initial flow solution are obtained the FLUENT and copled with droplet motion in the ambient condition.

The aerodynamic performance of aircraft in icing condition can deteriorate considerably by contamination of aerodynamic and propulsive systems due to icing accretions on aircraft surfaces. A computational analysis based on the Eulerian description was performed on an airfoil to investigate effects of ice accretions on airfoil aerodynamics. A water droplet with liquid water concentration ($0.00075kg/m^3$) and mean volume diameter ($20{\mu}m$) was considered and applied to various angles of attack to investigate the stall angle decrease and the drag increment.

In the present study, viscous flow calculations of helicopter main rotor system in forward flight were made by using an unstructured hybrid mesh solver. Each rotating blade relative to the cartesian frame was simulated independently by adopting unstructured overset mesh technique. For the validation of the present method, calculations for the Caradonna-Tung non-lifting forward flight and the AH-1G main rotor system in forward flight were made. Additional computation was made for the UH-60A rotor in forward flight. Reasonable agreements were obtained between the present results and the experiment.

The aerodynamic characteristics of aircraft winglet with structural deformation was investigated using the static FSI(Fluid-Structure Interaction) system. The system, comprised of CAD, CFD, CSD, VSI, and grid regeneration modules, was constructed. In the process VSI, grid regeneration, and integration modules were developed to combine CSD and CFD modules. As a test model, KC-135A, the double winglet suggested by Whitcomb, was selected and its aerodynamic characteristics for the rigid and deformable models was calculated by applying the static FSI system. As a result, the lift and drag coefficients of test models were reduced to 11% and 1.3%, respectively.

In order to reduce the download around Smart UAV(SUAV) at hovering and transition mode, flow control using synthetic jet has been performed. Many of the complex tilt rotor flow features are captured including wing leading and trailing edge separation, and the large region of separated flow beneath the wing. First, in order to control the trailing edge separation, synthetic jet is located at 30, 95% of flap chord length. The flow control using synthetic jet on flap shows that stall characteristics depending on several mode can be improved through separation vortices resizing. Also, a flap jet and a 0.01c jet which control the separation efficiently are applied at the same time at each test case because controlling the leading edge separation is essential for download reduction. As a result, time averaged download is reduced about 18% comparing with no control case at hovering mode and 48% at transition mode. These research results show that if flow control using leading edge jet and trailing edge jet is used effectively to the SUAV in overall flight mode, flight performance and stability can be improved.

Optimization with metamodel is one of numerical optimization methods. Response surface method is performed for making metamodel. The Hcks-Henne function is used for designing 2D shape of the airfoil and spring analogy is used to change the grid according to the change in shape of the airfoil. Aerodynamic coefficient required for response surface method are obtained by using Navier-Stokes solver with $\kappa-\omega$ shear stress transport turbulence model. For the baseline airfoils, OA 312, OA 309, and OA 407 airfoils select and optimize to improve aerodynamic performance.

An efficient and high-fidelity design approach for wing-body shape optimization is presented. Depending on the size of design space and the number of design of variable, aerodynamic shape optimization process is carried out via different optimization strategies at each design stage. In the first stage, global optimization techniques are applied to planform design with a few geometric design variables. In the second stage, local optimization techniques are used for wing surface design with a lot of design variables to maintain a sufficient design space with a high DOF (Degree of Freedom) geometric change. For global optimization, Kriging method in conjunction with Genetic Algorithm (GA) is used. Asearching algorithm of EI (Expected Improvement) points is introduced to enhance the quality of global optimization for the wing-planform design. For local optimization, a discrete adjoint method is adopted. By the successive combination of global and local optimization techniques, drag minimization is performed for a multi-body aircraft configuration while maintaining the baseline lift and the wing weight at the same time. Through the design process, performances of the test models are remarkably improved in comparison with the single stage design approach. The performance of the proposed design framework including wing planform design variables can be efficiently evaluated by the drag decomposition method, which can examine the improvement of various drag components, such as induced drag, wave drag, viscous drag and profile drag.

New design configurations were developed based on the closed ventilation type igloo nozzle in contrast to the conventional open distribution type. To validate the enhanced performance, the prediction of the exhaust flow patterns of the igloo diffuser was carried out by both experimental and computational simulations. CFD provided detailed parametric studies that significantly reduced the amount of experimentation necessary to alter the design of existing the F-15K igloo exhaust diffuser.

The helicopter development technology is improved by the increasing of computing power and advanced test facilities. The increasing efficiency of fuel consuming by the developing of improved lift-to-drag rotor system is the major issue, the noise reduction for ecology(civil area) and increase of survivability to reduce noise detection(battlefield) also are important. This investigation shows the helicopter external noise flight test methods, conventional military helicopter's(UH-60) experimental results and the numerical modeling method for aero-acoustic of rotor blade and the result of CAA(Computational Aero-Acoustic) for main rotor blade.

The present paper focuses on the analysis of aero-acoustics characteristic by several boundary conditions. In this simulation, a high-order and high-resolution numerical schemes are used for the accurate computation of compressible flow with several boundary conditions including characteristic boundary conditions as well as extrapolation and zonal characteristic boundary condition. These boundary conditions are applied to the computation of two dimensional circular cylinder flows with Mach number of 0.3 and Reynolds number of 400. The computation results are validated with measurement datum and other computation results for the Strouhal frequency of vortex shedding, the mean drag coefficient and root-mean-square lift for the unsteady periodic flow regime. Secondary frequency is predicted by three kinds of boundary conditions characteristic.

In this paper, the aerodynamic shape of air intake was investigated for the efficient cooling of electronic equipments in aircraft pod. As a first step, ESDU method was applied for the basic shape design of air intake considering the operational environments. The second step was to confirm the performance on design point, so the internal flow field of air intake was analyzed using a commercial Navier-Stokes code(FLUENT). And also the aerodynamic characteristics of internal flow at off-design condition was investigated with the variations of airflow rate. The results show that the air intake meets the requirement of target performance under the mission environments.

In the present paper, transition turbulence model is applied to the NACA64(3)618 and detailed flow features are studied. The turbulence model is sensitive to the boundary layer grid quality and y+ of the grid was limited to 1. The prediction of the transition region is dependent on the local flow condition. The pressure coefficient distribution of the transition turbulence model is compared with that of the fully turbulent mode and the drag distribution of the transition turbulence model was compared with that of the wind tunnel test.

In this study, a generic airfoil designed by the inverse method was evaluated with several candidate airfoils as a first step. Each airfoil was compared with respect to aerodynamic performance to meet the requirement of HALE(high altitude long endurance) aircraft. The second step was to optimize the candidate airfoil using the couple of optimization formulations to down select an optimum airfoil. For the analysis of low Reynolds number 2D flow, Drela's MSES was used. After comparing the aerodynamic results, the best airfoil was chosen to construct the baseline 3D wing. The Navier-Stokes code was used to evaluate the overall aerodynamic performance of designed wing with other wings. The results show that the designed wing has the best performance compared with other wings.

A numerical study of the flow characteristics inside a U-type circular tube is carried out in this paper. The numerical simulations carried out by using a Navier-Stokes code which is commercially available. Before detailed numerical simulations, validation of present numerical approach is made by comparing numerical solutions with experimental data. Numerical simulations are performed to study the effect of curvature on the flow characteristics inside a U-type tube. Numerical solutions show that a significant effect on the secondary flow structure in the cross section of the tube, especially in the curved section is shown when the curvature ratio, ratio of curvature to tube diameter, is smaller than about 3.5. As the curvature ratio decreases below 3.5, a counter rotating vortex is found below the primary vortex in the cross section of the tube. Another dramatic change of the flow structure is the formation of streamwise separation zone when the curvature ratio is decreased below 1.25.

The existing code which solves two-dimensional RANS(Reynolds Averaged Navier-Stokes) equations and 2-equation turbulence model equations was modified to enable numerical simulation of various supersonic flows. For this, various boundary conditions have been implemented to the code. Bleed boundary condition was incorporated into the code for calculating wall mean flow quantities. Furthermore, the boundary conditions for the turbulence quantities along rough surfaces as well as porous walls were applied to the code. The code was verified and validated by comparing the computational results against the experimental data for the supersonic flows over bleed region on a flat plate. Using the newly modified code, numerical simulations were performed and compared with other computational results as well as the experimental data for the supersonic flows over an oblique shock with a bleed region.

In this paper, flow simulation algorithms for utilizing unstructured hybrid meshes are introduced. First, various types of meshes are introduced. Advantages and disadvantages of each type of meshes are discussed. Unstructured hybrid mesh approach, that is best suited for high speed viscous flow simulation, is presented. Lastly, various types of flow simulations using unstructured hybrid meshes are introduced.

A two-dimensional unsteady inviscid flow solver has been developed for the simulation of complex geometric configurations on adaptive Cartesian meshes. Embedded condition was used for boundary condition and a predictor-corrector explicit time marching scheme was used for time-accurate numerical simulation. The Cartesian mesh generator, which was previously developed for steady problem, was used grid generation for unsteady flow. The solver was based on ALE formulation for body motion. For diminishing the effects of cut-cells, the cell merging method was used. Using cell merging method, it was eliminated the CFL constraints. The conservation problem, which is caused cell-type variation around region swept by solid boundary, was also solved using cell merging method. The results are presented for 2D circular cylinder and missile launching problem.

A moving mesh system is one of the critical parts in a computational fluid dynamics analysis. In this study, the RBF(Radial Basis Function) which shows better performance than hybrid meshes was developed to obtain the deformed grid. The RBF method can handle large mesh deformations caused by translations, rotations and deformations, both for 2D and 3D meshes. Another advantage of the method is that it can handle both structured and unstructured grids with ease. The method uses a volume spline technique to compute the deformation of block vertices and block edges, and deformed shape.

In this paper automatic 3D wing shape modeling program is introduced. The program is developed in Visual Basic based on Net Framework 3.5 environment by using CATIA COM Library, and it is used together with CATIA system to model 3D wings with or without flaps. With this program users can easily construct wing models by specifying geometry parameters which are usually design variables with the aid of easy-to-use GUI environment, and specifying sectional airfoil data is done either by using analytic shape functions such as NACA series airfoils or by providing input files with point data describing the airfoil shape. When all the input parameters are provided, users can either work further with the model in the CATIA system which would be automatically started by the program or save the resultant model in the format of users choice. Unstructured grid generation program is also briefly described which can make grid generation task for a 3D wing easy and efficient one when used together with the wing modeling program by choosing STL format as the model's output format.

CFD analysis of the fuel injection pattern and the flow field surrounding the liquid propellant injector of a small thruster is performed. A good agreement is shown with PIV test data for the initial configuration. Analysis on various injector shapes is performed to observe the effect of injector shape on the trajectories of liquid droplet. A various shapes of injector is investigated to enhance spray pattern of the small injector.

A preliminary thermal analysis is performed for the optical payload system of a geostationary satellite. The optical payload considered in this paper is GOCI(Geostationary Ocean Color Imager) of COMS of Korea. The radiative and conductive thermal models are employed in order to predict thermal responses of the GOCI on the geostationary orbit. According to the results of this analysis are as follows: 1) the GOCI instrument thermal control is satisfactory to provide the temperatures for the GOCI performances, 2) the thermal control is defined and interfaces are validated, and 3) the entrance baffle temperature is found slightly out its specification, therefore further detailed analyses should be continued on this element.

COMS (Communication, Ocean and Meteorological Satellite) is a geostationary satellite and has been developing by KARI for communication, ocean and meteorological observations. It will be tested under vacuum condition and very low temperature in order to verify thermal design of COMS. The test will be performed by using KARI large thermal vacuum chamber, which was developed by KARI, and the COMS will be the first flight satellite tested in this chamber. The purposes of thermal balance test are to correlate analytical model used for design evaluation and predicting temperatures, and to verify and adjust thermal control concept. KARI has plan to use heating plates to simulate space hot condition especially for radiator panels such as north and south panels. They will be controlled from 90K to 273K by circulating GN2 and LN2 alternatively according to the test phases, while the shroud of the vacuum chamber will be under constant temperature, 90K, during all thermal balance test. This paper presents thermal modelling including test chamber, heating plates and the satellite without solar array wing and Ka-band reflectors and discusses temperature prediction during thermal balance test.

The cavitating flow simulation is of practical importance for many engineering systems, such as marine propellers, pump impellers, nozzles, injectors, torpedoes, etc. The present work has focused on the simulation of cavitating flow past cylinders with strong side flow. The governing equation is the Navier-Stokes equation based on homogeneous mixture model. The momentum and energy equation is in the mixture phase while the continuity equation is solved liquid and vapor phase, separately. An implicit dual time and preconditioning method are employed for computational analysis. The results from the present solver have been in a fairly good agreement with the experimental data and other numerical results. After the code validation the strong side flow was applied to include the wake flow effect of the submarine.

We developed a numerical method to analyze the contact-line problems, incorporating a dynamic contact angle model. We used level set method to capture free surface. The method is applied to the analysis of dynamic behavior of a droplet in DC electrowetting. The result is compared with an experimental data and result of perturbation method.

The water droplet motion and the interaction between the droplets in a PEMFC air flow channel with multiple pores, through which water emerges, is studied numerically by solving the equations governing the conservation of mass and momentum. The liquid-gas interface is tracked by a level set method which is based on a sharp-interface representation for accurately imposing the matching conditions at the interface. The method is modified to implement the contact angle conditions on the walls and pores. The dynamic interaction between the droplets growing on multiple pores while keeping the total water flow rate through pores constant is investigated by conducting the computations until the droplet motion exhibits a periodic pattern. The numerical results show that the droplet merging caused by increasing the number of pores is not effective for water removal and that the contact angle of channel wall strongly affects water management in the PEMFC air flow channel.

The effect of the Dirichlet boundary condition for the redistance equation of level set method on the solutionof sloshing problem is investigated by adopting four Dirichlet boundary conditions. For the solution of the incompressible Navier-Stokes equations, P1P1 four-step fractional finite element method is employed and a least-square finite element method is used for the solutions of the two hyperbolic type equations of level set method; advection and redistance equation. ALE (Arbitrary Lagrangian Eulerian) method is used to deal with a moving computational domain. It has been shown that the free surface motion in a sloshing tank is strongly dependent on the type of the Dirichlet boundary condition and the results of broken dam and sloshing problems using various Dirichlet boundary conditions are discussed and compared with the existing experimental results.

Shock focusing is related with explosive release of shock wave energy on a narrow spot in a short duration of time triggering a spontaneous high pressure near the focal point. It is well known that reflection of planar incident shock wave from the metallic concave mirror such as ellipsoidal, paraboloidal or hemispherical cavities will focus on a focal point. We intend to improve the computational results using a wave propagation algorithm and to resolve the mushroom-like structure. For computation of the concave cavity flow, it is not easy to use a single-block mesh because of the many singular points in geometry and coordinates. We have employed a uniform Cartesian-grid method for the wave propagation algorithm.

In the present study, a numerical simulation for the diffusion of hydrogen jet in a enclosure was performed to aid the leakage test of the hydrogen for the safety of the hydrogen vehicle. The temporal and spatial distributions of the hydrogen concentration in the test chamber are predicted from the present numerical analyses. Flammable region of 4-74% and explosive region of 18-59% hydrogen by volume was identified from the present results. Factors influencing the diffusion of the hydrogen jet were examined to evaluate the effectiveness of forced ventilation for relieving the accumulation of the leaked hydrogen gas in the chamber, which include location of open windows, size of leakage nozzle, and leakage rate among others. The distribution of the concentration of the leaked hydrogen for various cases can be used as a database in various applications for the hydrogen safety.

The objective of this study is to get simulation data about pulsatile flow around an interior solid body inside a bifurcated tube. All the processes were based on CFD method, with a commercial FVM code, SC/Tetra ver. 6.0 for solving, and with CATIA R16 for generating geometries. The bifurcated tube models were drawn with the bifurcated angle of 45 degrees, considering Murray's law about the diameter ratio. With various locations of the object, the effects of flow on the drag were considered. For the pulsating flow condition, the velocity wave profile was given as the inlet boundary condition. To validate all the result, the simulation was compared with the existing data of the other papers first. Overall flow field of both data were similar, but there was some difference at a zero velocity. Therefore the next simulation was continued with the sine wave profiles where there is no negative flow, and then the data was compared with one of the pulmonary artery velocity where there is negative flow. The final process was to calculate flow variables such as the wall shear stress (WSS) and to compute the drag of the solid object.

A parametric study for for the water entry of a two dimensional symmetric wedge with deadrise angle of 10 degrees was carried out to find out the most dominant parameter. Water entry problem with constant velocity is simplified as the stationary wedge in the way of the upcoming water surface. The calculated impact loads showed that the effect of the viscosity was not so important in this problem. For a given grid system a suitable time step size can be found. The most sensitive parameter was found to be the grid size.

Global warming issues due to the $CO_2$(Carbon Dioxide) become increasing since the Industrial Revolution. After the Kyoto protocol at 1997, nations which have the prearranged quota drives their national project for the reduction of $CO_2$. Korean Government start to the related big projects in the view of three concepts which have consist of the $CO_2$ exhaust reduction on land, $CO_2$ capture and $CO_2$ storage. Furthermore, the storage method putting into depleted region underground is accepted by the London Convention while the ocean diluted method discharging the liquid $CO_2$ into the deep ocean using the long pipe which is towed by the surface vessel is underway for the research steps which means that there are many potentials for the R&Ds that need for the breakthrough. In this paper, the role and example of the Computational Fluid Dynamics for the feasibility study of the $CO_2$ ocean sequestration is mentioned.

본 연구에서는 대심도 지하역사에서 화원 위치에 따른 연기거동의 특징을 분석하였다. 전산수치해석에 FDS code가 사용되었다. 화재에 의한 유체거동을 모사하기 위하여 난류모델은 LES를 사용하였으며, 빠른 계산을 수행하기 위하여 병렬수치해석기법을 사용하였다. 본 연구를 통하여 대심도 역사에서 화재시 피난을 저해하는 화원의 위치에 대하여 검토하였다.

지하 복수층 주차장은 도시의 밀집화 및 자동차의 급격한 증가로 인한 토지부족 및 주차난을 해소할 수 있는 좋은 방법이다. 하지만 인체유해 배출가스의 집중으로 환기의 중요성이 대두되는 추세에도 불구하고 실제로는 시공이 완료 된 후에나 오염도 측정이 가능하다. 이러한 이유로 CFD 시뮬레이션으로 시공이전의 지하 복수층 주차장의 공기오염도를 예측했다. CFD 시뮬레이션은 예산과 시간을 절감한 공기오염도 예측이 가능하다. 실제의 당산근린공원의 지하주차장을 대상으로 오염도를 예측했다. 급기에만 강제유동을 일으키는 2종과 배기에만 강제유동을 일으키는 3종 기계환기법, 그리고 급배기에 강제유동을 일으키는 1종기계환기법을 각각 적용하여 공기오염도를 계산했다. 지하1층과 지하2층의 입구 및 출구유속을 10, 15m/s로 하여 실내평균연령과 실내평균환기효율을 살펴보았다. 계산결과는 급기와 배기에 강제유동을 동시에 일으키는 1종 기계환기법의 값과 비교하였다.

In the present study, a numerical analysis on the sloshing in a tank with the harmonic motion was investigated. A VOF method was used for two-phase flows inside the sloshing tank and a source term of the momentum equation was applied for the harmonic motion. This numerical method was verified by comparing its results with the available experimental data. The sloshing in a tank causes the instability of the fluid flows and the fluctuation of the impact pressure on the tank. By these phenomena of the tank sloshing, the sloshing problems such as the failure and the noise of system can be generated. For the reduction of these sloshing problems, the various baffles such as the horizontal/vertical plate baffles and the porous baffles inside the tank are installed. With the installations of these baffles, the characteristics of the liquid behavior in the sloshing tank, the impact pressure on the wall, the amplitude of the free surface near the wall and the sloshing noise were numerically analyzed.

This paper presents three-dimensional flow analysis for a mixed-flow pump which consists of a rotor and a stator. Reynolds-averaged Navier-Stokes equations with shear stress transport turbulence model are discretized by finite volume approximations and solved by the commercial CFD code CFX 11.0. Structured grid system is constructed in the computational domain, which has O-type grids near the blade surfaces and H-type grids in other regions. Validation of the numerical results was performed with experimental data for head coefficients and hydraulic efficiencies at different flow coefficients. This paper shows that the pump characteristics can be predicted effectively by numerical analysis.

Traditional vacuum cleaner adoptsa highly rotating centrifugal impeller for generating suction region at lower pressure. The rotational speed is very high at 35,000 to 40,000 rpm and accessory structure such as a circular brush at the base plate of the cleaning devices is equipped for gathering dusts. Highly rotating impeller is effective for low pressure generation but causes noise problems. Recently, the unwanted noise is removed by installation of air-pump, instead of a centrifugal impeller, and the internal flows of the modified system are estimated in numerical and experimental approach, respectively.

The objective of this study is to investigate the suitable design for a domestic CWP pump, which is used in cooling-water intakes for the unit 3 and 4 of Yeonggwang nuclear power plant. All the simulations are performed, using CFD method with a commercial code STAR-CCM+ version 3.02. After modeling a present design of the pump, the flow around the rotating blade was calculated by using quasi-static method and sliding mesh method with the almost same condition as an actual state. Based on fundamental simulations with various depth of sea water, the reference pressure for the boundary condition of the present study was decided. To verify the reliability of the calculation results, the suction flow rate of the data was compared with that of the experimental data. As a result of this comparison, it is confirmed that two results are fairly consistent. For the improvement of the suction flow rate, computational analysis was done by changing a flow channel and blade shapes. It is shown that the suction flow rate of the new pump was improved.

The gridless (or meshfree) methods, such as MPS, SPH, FPM an so forth, are feasible and robust for the problems with moving boundary and/or complicated boundary shapes, because these methods do not need to generate a grid system. In this study, a gridless solver, which is based on the combination of moving least square interpolations on a cloud of points with point collocation for evaluating the derivatives of governing equations, is presented for two-dimensional unsteady incompressible Navier-Stokes problem in the low Reynolds number. A MAC-type algorithm was adopted and the Poission equation for the pressure was solved by successively in the moving least square sense. Some weighing functions were tested in order to investigate the up-winding effect for the convection term. Some typical problems were solved by the presented solver for the validation and the results obtained were compared with analytic solutions and the numerical results by conventional CFD methods, such as FVM.

Numerical algorithms for solving the incompressible Navier-Stokes equations using P2P1 finite element are compared regarding the eigenvalues of matrices. P2P1 element allocates pressure at vertex nodes and velocity at both vertex and mid nodes. Therefore, compared to the P1P1 element, the number of pressure variables in the P2P1 element decreases to 1/4 in the case of two-dimensional problems and to 1/8 in the three-dimensional problems. Fully-implicit-integrated, semi-implicit- integrated and semi-segregated finite element formulations using P2P1 element are compared in terms of elapsed time, accuracy and eigenvlue distribution (condition number). For the comparison,they have been applied to the well-known benchmark problems. That is, the two-dimensional unsteady flows around a fixed circular cylinder and decaying vortex flow are adopted to check spatial accuracy.

In the present study, a grid deformation technique has been incorporated into the unsteady compressible and incompressible viscous flow solvers on unstructured hybrid meshes. An algebraic method based on the basis decomposition of normal edge vector was used for the deformation of viscous elements, and a ball-vertex spring analogy was adopted for inviscid elements among several spring analogy methods due to its robustness. The present method was validated by comparing the results obtained from the grid deformation and the rigid motion of entire grids. Fish swimming motion of an NACA0012 airfoil and flapping wing motion of a generic fighter were simulated to demonstrate the robustness of the present grid deformation technique.

The local preconditioning method has both robust convergence and accurate solutions by using local flow properties for parameters in the preconditioning matrix. Preconditioning methods have been very effective to low speed inviscid flows. In the viscous and turbulent flows, deterioration of convergence should be overcame on the high aspect ratio grids to get better convergence and accuracy. In the present study, the local time stepping and min-CFL/max-VNN definitions are applied to compare the results and we propose the method that switches between two methods. The min-CFL definition is applied for inviscid flow problems and the min-CFL/max-VNN definition is implemented to viscous and turbulent flow problems.

A temperature equation which is derived from an enthalpy transport equation by using an assumption of a constant specific heat is very attractive for analyses of heat and fluid flows. It can be used for an analysis of a solid-fluid conjugate heat transfer, and it does not need a numerical method to find temperature from a temperature-enthalpy relation. But its application is limited because of the assumption. A new method is derived in this study, which is a temperature-explicit formulation of the energy equation. The enthalpy form of the energy equation is used in the method. But the final discrete form of the equation is expressed with temperature. It can be used for a solid-fluid conjugate heat transfer and multiphase flows. It is found by numerical tests that it is very efficient and as accurate as the standard enthalpy formulation.

In level set method, material properties are made to change smoothly across an interface of two materials with different properties by introducing an interpolation or smoothing scheme. So far, the weighted arithmetic mean (WAM) method has been exclusively adopted in level set method, without complete assessment for its validity. We showed here that the weighted harmonic mean (WHM) method for rate constants of various rate processes, including viscosity, thermal conductivity, electrical conductivity, and permittivity, gives much more accurate results than the WAM method. The selection of interpolation scheme is particularly important in multi-phase electrohydrodynamic problems in which driving force for fluid flow is electrical force exerted on the phase interface. Our analysis also showed that WHM method for both electrical conductivity and permittivity gives not only more accurate, but also more physically realistic distribution of electrical force at the interface. Our arguments are confirmed by numerical simulations of drop deformation under DC electric field.

가압 경수로의 주요 기기에서 발생할 수 있는 과도 2상 유동(Two-phase flow) 현상에 대한 해석을 수행하기 위해 원자로 기기 열수력 해석 코드를 개발 중에 있다. 개발 중인 기기 열수력 해석 코드는 지배 방정식으로 Two-phase, three-field model을 사용하고 있으며, 복잡한 기하학적 형상의 원자로 기기를 모사하기 위해 비정렬 격자계(Unstructured grid)를 활용하고 있다. 수치해석 기법으로는, 원자로 계통 해석코드 RELAP5가 사용 중이며 대부분의 원자로 내 2상 유동 조건에서 안정적이며 정확하다고 알려진 Semi-implicit 방법을 적용하였다. 그러나 기존의 Semi-implicit 방법은 1차원, 엇갈림격자(Staggered grid)에 대해 개발되었기 때문에, 이를 다차원, 비정렬, 비엇갈림 격자(Non-staggered grid)에 적용하기 위해 기존의 Semi-implicit 방법을 수정하였다. 본 논문에서는 Semi-implicit 방법의 대류항을 이차정확도를 갖도록 확장하였으며, 이차정확도에 의한 수치확산의 감소를 평가하기 위해 수행된 수치시험의 결과를 기술하였다. 이차정확도 및 일차정확도로 계산된 값을 해석해 또는 격자 수렴성 시험을 통해 평가해 본 결과, 이차정확도 계산시 수치 확산의 감소 확인하였다.

Numerical analysis of two dimensional incompressible laminar free jet flow was carried out by using finite difference SMAC scheme. Flow characteristics of free jet flow such as jet width, similarity of jet velocity and hypothetical origin were investigated for different Reynolds numbers of Re=30 and 100. The reliability of predictions were confirmed by comparison with exact solution. Non-dimensional velocity distribution showed similarity of jet flow velocity after initial region. In the region of laminar flow, the location of hypothetical origin becomes more distant with Reynolds number.

In this study, the heat transfer and flow field of condenser used on Kim-chi refrigerator is analysed with numerical method. Main objective is to present the base data for designing new condenser model with improvement of heat transfer performance. For CFD analysis, a commercial code, STAR-CCM+ was used. The water was used for the inner working fluid and the air was used for the outer fluid. The condenser type used in this study is a flat plate fin-and-tube heat exchanger. As factors for performance analysis, the effect of condenser geometry and air velocity was investigated. As a result, it has been observed that there is a suitable fin pitch with which heat transfer performance of condenser is maximized.

A numerical code for the interior ballistics has been investigated. The method of the ghost cell extrapolation has been used for the moving boundary with the projectile movement. The porosity effect and the Ergun's Equation have been used in the numerical calculation for the grain combustion. The calculation results of the numerical code have been compared and verified through those of the lumped parameter method. Computerization techniques of the numerical analysis for the interior ballistics have been developed.

A computational investigation is conducted on free convection from a thin plate having a surface heat source. The thermal configuration simulates the recently-proposed transparent film heater made of a single-walled carbon nanotube film on a glass substrate. The Navier-Stokes computations are carried out to study laminar free convection from the heater. Parallel numerical experiments are performed by using a simplified design analysis model which solve the conduction equation with the boundary conditions utilizing several existing correlations for convective heat transfer coefficient. Comparison leads to the most suitable boundary condition for the thermal model to evaluate the performance evaluation of a transparent thin-film heater.

Free surface impinging jet on a moving plate, which is applicable to cooling of hot metals in a steel-making process, is investigated numerically by solving the Navier-Stokes equations in the liquid and gas phases. The free surface of liquid-gas interface is tracked by a level-set method which is improved by incorporating the ghost fluid approach based on a sharp-interface representation. The method is further improved by employing a nonequilibrium $\kappa-\varepsilon$ turbulence model including the effect of low Reynolds number. The computations are made to investigate the effects of the nozzle pitch, moving velocity of plate and jet velocity on the interfacial motion and the associated flow and temperature fields.

An oxy-fuel boiler has been developed to capture $CO_2$ from the exhaust gas. FGR (flue gas recirculation) is adopted to be compliant with the retrofit scenario. Numerical simulations have been performed to study the detailed physics inside the combustion chamber of the boiler. The temperature field obtained from the simulation agrees with the flame image from the experiment. The FGR combustion yields similar heat transfer characteristics with the conventional air combustion while the flame is formed further downstream in case of the FGR combustion.

Operating medium or large scale industrial boilers in partial load condition, the burner should undergo the off-design points resulting in poor exhaust gas characteristics. To obtain the stable turn down performance, two or more burners can be used for the industrial boiler. In case multiple burners are adopted, the heat transfer can be enhanced by arranging the burners properly. In the present study, numerical simulations have been conducted for the combustion chamber of a 2 t/h class industrial boiler in order to clarify the heat transfer characteristics at the combustion chamber.

The purpose of the present study is to investigate the flow resistance and the heat transfer characteristics of oval fin-tube heat exchanger. Six kinds of oval fin having the same fin area and the different diameter ratio tested numerically. Test data for the heat transfer, pressure drop and fin temperature were shown and discussed. The pressure drop and heat transfer increased for increasing the oval fin diameter ratio(diameter of span-wise direction to diameter of longitudinal diameter) up to 50% and 45% respectively.

Plate Heat Exchnager(PHE) has recently become widely adopted for cogeneration systems owing to its small installation space and high thermal efficiency. The gap between plates can be changed depending on its assemble type, i.e. gasket or blazing. The gap is known to affect thermal efficiency and working pressure drop in PHE with complicated geometrical features. Numerical simulation techniques have been developed to deal with PHE with complex configuration of chevron plates. The present study is aiming at identifying the gap effect on pressure drop and thermal efficiency of the PHE. The numerical simulation results show that the gap has relatively large effects on working pressure drop than thermal efficiency in performance of PHE.

Most material of engineering interest undergoes solidification process from liquid to solid state. Identifying the underlying mechanism during solidification process is essential to determine the microstructure of material which governs the physical properties of final product. In this paper, we expand our previous two-dimensional numerical technique to three-dimensional simulation for computing dendritic solidification process with fluid convection. We used Level Contour Reconstruction Method to track the moving liquid-solid interface and Sharp Interface Technique to correctly implement phase changing boundary condition. Three-dimensional results showed clear difference compared to two-dimensional simulation on tip growth rate and velocity.

Recently, the manufacturing of micro-cavity by means of laser surface texturing (LST) technique and low friction study by the LST have been in great progress. Most of current works have been dealing with the effect of cavity on friction and wear. The main objective of the present study was to investigate numerically two-dimensional lubrication characteristics of micro-dimple shapes fabricated on solid surfaces, and this study utilized the commercial CFD code (Fluent V.6.3). For the evaluation, preliminary simulation was conducted and numerical predictions were compared with the analytic solution obtained from the Reynolds's equation. Mainly, the present study investigated the influence of dimple depth, pattern shapes, and film thickness on lubrication characteristics related to the reduction of friction. It is found that the existence of micro-dimpled surface makes it possible to substantially reduce the friction forces exerted on the surfaces. In particular, substantial decrease in shear stresses was observed as the lubricant film thickness decreases. For instance, in the case of the film thickness of 0.01 mm, the estimated shear stress decreases up to about 40%. It indicates that the film thickness would be important factor in designing the micro-dimpled surfaces. Furthermore, it was observed that such a optimum dimple depth would be present because the dimple depth larger than the optimum value did no longer affect the reduction in shear stresses.

Studies on the effect of the wall-ion, wall-water, water-ion and ion-ion interaction on properties of water and ions in nano-channels have been performed through the use of different kinds of ions or different models of potential energy between wall-ion or wall-water. On this paper, we address the effect of water-wall interaction potential on the properties of confined aqueous solution by using the molecular dynamics (MD) simulations. As the interaction potential energies between water and wall we employed the models of the Weeks-Chandler-Andersen (WCA) and Lennard-Jones (LJ). On the MD simulations, 680 water molecules and 20 ions are included between uniformly charged plates that are separated by 2.6 nm. The water molecules are modeled by using the rigid SPC/E model (simple point charge/Extended) and the ions by the charged Lennard-Jones particle model. We compared the results obtained by using WCA potential with those by LJ potential. We also compared the results (e.g. ion density and electro-static potential distributions) in each of the above cases with those provided by solving the Poisson-Boltzmann equation.