Fluid flows within microchannels are characterized by low Reynolds numbers. Therefore the effect of mixing is a crucial factor in design of the channels. Since the action of the electro-osmotic or magnetic forces used in the mixing enhancement is usually periodic in the three-dimensional channel configuration, use of the various concepts of chaotic advection is reasonable in the quantification of the stirring effect. In this paper, the details of the method of material-stretching mapping is explained. The actual application of the method to the screw extruder is also presented.

The paper introduces the state-of-the art of WIG ship and application field of the CFD to WIG ship development. WIG ships are highly efficient and fast transport vehicles which take advantage of benefits of ground effect. Chapter 1 describes the characteristics of WIG ship. Chapter 2 reviews the research works of Russia, German, Chinese, Korea and etc. Chapter 3 explains the kind of WIG ship categorized by the main operational mode and take-off system. Chapter 4 describes about the application field of CFD to WIG ship development procedure.

Numerical analysis were made for the unsteady flow fields of the rotor system of a Tilt-Rotor aircraft in cruise mode. The Reynolds-averaged thin-layer Wavier-Stokes equations were discretized by Roe's upwind differencing scheme and integrated in time by the LU-SGS algorithm. The computational domain of the rotor system was constructed by seven multi-block Chimera grids. Comparison of pressure coefficient on the surface of the main wing and blades were made for 3cases of advance ratio(0.325, 0.350, 0.375) and thrust and power coefficients for the rotor were compared with experimental data.

Axial fan is used for the supplement of large amount of flows. Axial blowers show relatively high efficiency of the system. The present model of axial fan is for cooling a condenser in an air-conditioning unit that exhibits tendency toward compact size. In order to realize the compact model, the width of an axial blade should be cut down in axial distance. Main interest lies on the performance of the axial blowing system with blades having shorter chord length. One of the important design parameters for axial fan is the shape of the blades of it. Design of blades includes the cross-sectional shape and its dimension, including the chord length. We consider two types of blades; one is NACA airfoil with normal chord length and the other is with shortening chord length by $10\%$ of normal airfoil. Axial blower with the modified blades is essential for the compact model of an air-conditioner. The other design parameters are same in the two cases. Using a wind tunnel follows ASHRAE standards carries out evaluation of performance of the system. Detail of flows around the blades is prepared by velocity measurements using PIV. According to performance estimation, the axial blower with short chord blade show quite close to the performance results, including flow rate and pressure rise, of the standard one. The reason of the two similar results is that the flowpatterns depend on Reynolds number based on the chord length of a blade. In this investigation, the critical chord length is found, in which the flows near the airfoil are so unstable and the performance of the system is decreased. A series of figures is for the detail information on the flow.

고속 선박을 추진하는 한 방법으로 널리 사용되는 물분사 추진은 물을 내부 덕트로 빨아들여 임펠러로 물을 가속시켜 노즐을 통해 분사시킴으로써 입출구의 운동량차이에 의해 추력을 얻는 추진장치이다. 선박의 목적에 따라 사용되는 다양한 형태의 물분사 추진기의 개발을 위하여 모형실험을 통하여 그 성능을 검증하는 부분에서 로터의 피치각 변화에 따른 추진기의 성능 실험을 하는 것은 많은 비용과 시간이 따른다. 따라서 이러한 문제를 해결하기 위하여 본 연구에서는 추진기 내부의 유동장을 4가지 피치각에 따라 추진력을 3차원 비압축성 Navier-Stokes 방정식을 이용하여 해석하였다. 로터의 회전을 고려하여 슬라이딩 다중 격자기법을 적용하였고 추력계수, 토크계수, 그리고 모멘텀을 해석 결과와 비교 분석을 통하여 추진기의 성능과 효율을 추정하였다.

Up to recently, with the improvement of a computer power and high speed of network technology, advanced countries have researched a construction of the e-Science environment. As a major application part, a construction for environment of CFD, also, have studied together. During the research, people realize that not sharing hardware but also appropriate software development is really important to realize the environment. This paper describes about a construction of a collaborative research environment in the KISTI: Clients can connect to the computing resources through the web portal, run the Cactus simulation.: According to the computing resources, the simulation can migrate to some site to find better computing power.: Result of the calculation visualize at the web portal directly so that researchers of remote site can be share and analyze the result collaborative ways.

A new PC cluster was designed and constructed based on Windows XP Operating system. Primary target of the present design was the high node density per rack by using the general PC parts those are cost-effective and readily available in the market. Other major design points were system cooling and the convenient maintenance using standard PC parts. Presently 24 nodes per rack seems to be optimum considering the specification of the network switching device, system cooling and power supply, but 40 nodes can be accommodated within a single rack at maximum. Windows XP was selected as a high-performance computing environment considering the cost and the convenience in acquisition, maintenance and education. Both fast-Ethernet and Gigabit Ethernet network connection were tested and compared with previous data, especially for Linux doter using Myrinet. The result shows that there is no significant difference between the operating systems and the Fast-Ethernet and/or Gigabit Ethernet are good solution for the high-performance PC cluster considering the cost and performance.

The Grid'[1] means the collaboration of computing and experimental resources in dispersed organizations by high-speed network. It has been paid much attention for an unlimited number of potential resources available and the easiness to build collaborative environments among multiple disciplines. However, the difficulty in establishing the environments and accessing and utilizing the resources has prevented application scientists from conducting Grid computing. Thus, the present study focuses on building PSE(Problem Solving Environment) which assists application researchers to easily access and utilize the Grid. The Cactus toolkit, originally developed by astrophysicists, is used as a base frame for Grid PSE. Some modules are newly developed and modified for CFD(Computational Fluid Dynamics) analysis. Simultaneously, a web portal, Grid-One portal, is built for remote monitoring/control and job migration. Cactus frame through the web portal service has been applied to various CFD problems, demonstrating that the developed PSE is valuable for large-scaled applications on the Grid.

A new feedback control system based on system identification is proposed and preliminarily tested on Van der Pol equation which has a similar characteristic to circular cylinder. The same principle is applicable to circular cylinder in a uniform flow for suppresing the vortex shedding. The feedback controller is designed to impose feedback signal at the phase which is located outside the range of lock-on. The lift coefficient (CL) is employed as a feedback signal and the control forcing is given by a rotational oscillation of the cylinder. By applying the feedback control system, the lift coefficient is reduced.

Aerodynamic characteristics of three-dimensional wings in ground effect for Aero-levitation Electric Vehicle(AEV) are numerically investigated for various fairing shapes at the junctions of 3D Wings. Numerical results show that a sizeable three-dimensional comer flow separation occurs with formation of an arch vortex at the junction of main and vertical wings, and also that this is predicted the main cause of the high lift-to-drag(L/D) reduction rate of the main wing. To avoid the comer flow separation, the main idea of this study is to reduce the cross section gradient of the comer flow tube near the trailing edge for various fairing shapes. Improvements on L/D ratios of the wings are pursued by breaking the coherence of superimposed adverse pressure gradients at the wing junction when the cross section gradient is changed slowly at the trailing edge.

Curretly, the development of MEMS(Micro Electronic Mechanical System) technology awakes many research's interest for the aerodynamics. This work presents the development of a compact synthetic jet actuator for flow separation control at the flat plate. The formation and evolution of fluidic actuators based on synthetic jet technology are investigated using Reynolds-Averaged Navier-Stokes equations. Also, 2-Dimensional, unsteady, incompressible Navier-Stokes equation solver with single partitioning method for Multi-Block grid to analyze and a modeled boundary condition in developed fo. the synthetic jet actuator. Both laminar and turbulent jets are investigated. Results show very good agreement with experimental measurements. A jet flow develops, even though no net mass flow is introduced. Pair of counter-rotating vortices are observed near the jet exit as are observed in the experiments.

The HANARO, multi-purpose research reactor, 30 MWth open-tank-in-pool type, is planning to produce a fission moly-99 of radio isotopes, a mother nuclide of Tc-99m, a medical isotope and is under developing a target handling tool for loading and unloading it in a circular flow tube (OR-5). A guide tube is extended from the reactor core to the top of the reactor chimney for easily loading the target under the reactor normal operation. But active coolant through the core can be quickly raised up to the top of the chimney through the guide tube. The jet flow was suppressed in the guide tube after reducing the inner diameter of a flow restriction orifice installed in the OR-5 flow tube for adding the pressure difference in the flow tube after unloading the target. This paper describes an analytical analysis to calculate the flow distribution in the core of the HANARO after suppressing the jet flow of the guide tube. As results, it was confirmed through the analysis results that the flow distribution in the core of the HANARO were not adversely affected.

A numerical investigation has been performed to discuss the radiation-affected steady-laminar natural convection in an enclosure under a large temperature difference. Due to inherent nature of this study, the Boussinesq approximation is no longer valid. Therefore the radiating fluid in an enclosure is treated as a ideal gas. To examine the effects of thermal radiation on thermo-fluid dynamic behaviors in complex geometries, two incomplete partitions are introduced. Based on the results of this study, the dispositions of incomplete partitions with radiatively participating medium are found to incur a distinct difference in fluid-dynamic as well as thermal behavior.

The numerical solutions are examined on the effect of a centered heat conducting body on natural convection in a 2-D square cavity. The influences of the Rayleigh number, the dimensionless conducting body size, and the ratio of the thermal diffusivity of the body to that of the fluid have been investigated on the natural convection heat transfer in overall concerned region. The analysis reveals that the fluid flow and heat transfer processes are governed by all of them. Results for isotherms, vector plots and wall Nusselt numbers are reported for Pr = 0.71 and relatively wide ranges of the other parameters. Heat transfer across the cavity, in comparison to that in the absence of a body, are enhanced (reduced) in general by a body with a thermal diffusivity ratio less (greater) than unity. The heat transfer are also found to attain a minimum as the body size is increased.

The heat transfer enhancement by pulsatile flow in plate heat exchanger has been investigated numerically in the present study. The numerical study was performed ill the range of the Strouhal number from 0.04 to 2 and the Reynolds number from 370 to 730. The results showed that the pulsatile flow produces resonating vortex shedding at the groove sharp edges and a strong transient vortex rotation within the grooved channels. As a result, the mixing between the trapped volume in the grooved cavity and the main stream was enhanced. Good agreements between the predictions and measured data are obtained for the optimum frequency of pulsation and corresponding heat transfer enhancement

A heat dissipation modeling method of EEE parts is developed for thermal design and analysis of an satellite electronic equipment. The power consumption measurement value of each functional breadboard is used for the heat dissipation modeling method. For the purpose of conduction heat transfer modeling of EEE parts, surface heat model using very thin ignorable thermal plates is developed instead of conventional lumped capacity nodes. The thermal plates are projected to the printed circuit board and can be modeled and modified easily by numerically preprocessing programs according to design changes. These modeling methods are applied to the thermal design and analysis of CTU and verified by thermal cycling and vacuum tests.

The geostationary satellite propulsion system has thermistors which can measure liquid propellant temperature at tanks, pipes and etc. In the satellite propulsion system with several tanks, the propellant in the tanks is moved by temperature change and this temperature pattern is constant. In this paper, the temperature change pattern of KOREASAT 1 propulsion system is compared and the prediction study of pressurant inflow using temperature change of geostationary satellite propulsion system is described.

This work discusses the numerical analysis, the design and experimental test of the SMA (Shape Memory Alloy) actuator along with its capabilities and limitations. Convection heating and cooling using water actuate the SMA element of the actuator. The fuel such as propane, having a high energy density, is used as the energy source for the SMA actuator in order to increase power and energy density of the system, and thus in order to obviate the need for electrical power supplies such as batteries. The system is composed of a pump, valves, bellows, a heater (burner), control unit and a displacement amplification device. The actuation frequency is compared with the prediction obtained from numerical analysis. For the designed SMA actuator system, the results of numerical analysis were utilized in determining design parameters and operating conditions.

The Korean high speed train runs at 300 km/h, ballast-flying phenomenon often happens by strong train-wind. It is important to consider the prevention of ballast-flying phenomenon, because the train under-body and fares or walker around a rail might be damaged. In this study, Numerical analysis of the flow field of under-body of train and study of heighter-effect were conducted to decrease the speed of under-body. The shape of under-body was simplified for convenience of meshing and analysis. According to results of Taguchi's design by orthogonal arrays, a height of tie is dominant in the flow field, so if the heighter is installed on tie, the speed of under-body might be decreased. The result of this study is useful to build a new high-speed-line.

Recently, the weight of train is decreased by using the light material for improvement in energy efficiency. And the length of whole train is more increased for mass transportation of passengers and cargo. However, decrease of the weight and increase of the length of train can cause the train to be overturned or derailed by strong crosswind. In case of Korean Tilting Train eXpress (TTX), the situation can be more severe. TTX will be developed for a quasi-high speed train at 200km/h speed rate and operated on the existing tracks. Moreover, the weight of TTX will be much less than that of conventional train. It is supposed that TTX will be very sensitive to crosswind. In this paper, numerical analysis is used to investigate aerodynamic characteristics around TTX and obtain the induced lateral force by crosswind. After calculating derailment coefficient and overturning coefficient using numerical results, the crosswind safety of TTX is judged. This paper will be good data for judging crosswind safety of TTX.

CFD simulation was peformed for 2D and 3D flying and rotating frisbees. Multiple reference method(MRF) was utilized to consider the rotation of 3D model. Geometry change of 2D model shows dramatic increase of lift, but 3D simulation results for geometry change show decrease of lift and drag. Ground effect increases the lift of the frisbee being close to ground.

In this paper a method for analysing aerodynamic characteristics of a rotating parachute in steady descending motion is presented. Because of a complex geometric configuration of the parachute associated with side vents and discontinuous skirts, special procedure was adopted th handle the geometry in the analysis. A panel method was successfully applied to the present problem and yielded good results using above procedure. A CFD code using the full Navier-Stokes equations was also applied and provided good results. Parachute free drop and wind tunnel tests were performed to determine the fully developed canopy configuration and aerodynamic characteristics. The method can be used for optimizing the parachute size and side vent configurations in the design period.

The numerical analyses of the complicated flows are widely attempted in these days. Because of the enormous demanding memory and calculation time, parallel processing is used for these problems. In order to obtain calculation efficiency, it is important to choose proper domain decomposition technique and numerical algorithm. In this research we enhanced the efficiency of the CFDS code developed by ADD, using parallel computation and newly developed numerical algorithms. For the huge amount of data transfer between blocks non-blocking method is used, and newly developed data transfer algorithm is used for non-aligned block interface. Recently developed RoeM scheme is adpoted as a spatial difference method, and AF-ADI and LU-SGS methods are used as a time integration method to enhance the convergence of the code. Analyses of the flows around the ONERA M6 wing and the high angle of attack missile configuration are performed to show the efficiency improvement.

In the CFD area, the numerical analysis of high Mach number flow over a blunt-body poses many issues. Various numerical schemes have been developed to cover the issues, but the traditional schemes are still used widely due to the complexities of new schemes and intricacy of modifying the established codes. In the present study, the well-known Roe's FDS based on TVD-MUSCL scheme is used for the solution of very high Mach number three-dimensional flows posing carbuncle and non-physical phenomena in numerical analysis. A parametric study was carried out to account for the effects of the entropy fixing, grid configurations and initial condition. The carbuncle phenomena could be easily overcome by the entropy fixing, and the non-physical solution could be eliminated by the use of the modified initial condition regardless of entropy fixing and grid configurations.

A hybrid central-WENO scheme is proposed. The fifth order WENO-LF scheme is coupled with a central flux scheme at cell face. Two sub-schemes, the WENO-LF scheme and the central flux scheme, are switched by a weighting function. The efficiency and accuracy of the proposed hybrid central-WENO scheme is validated through several numerical experiments.

A numerical study of evaluation of turbulence models for thermal striping phenomenon is performed. The turbulence models chosen in the present study are the two-layer model, the shear stress transport (SST) model and the V2-f model. These three models are applied to the analysis of the triple jet flow with the same velocity but different temperature. The unsteady Reynolds-averaged Navier-Stokes (URANS) equation method is used together with the SIMPLE algorithm. The results of the present study show that the temporal oscillation of temperature is predicted only by the V2-f model, and the accuracy of the mean velocity, the turbulent shear stress and the mean temperature is a little dependent on the turbulence model used. The the two-layer model and the SST model shows nearly the same capability of predicting the thermal striping and the amplitude of the temperature fluctuation is predicted best by the V2-f model.

The three-dimensional incompressible flow past an open cavity in a channel is investigated using Detached Eddy Simulation(DES). The length to depth ratio of the cavity is 2 and the Reynolds number defined with the cavity depth is 3,360. The DES methods are based on the Mentor's SST model. In the present work, two types of inflow conditions are used; one is RANS profile, the other is LES inflow from another Large Eddy Simulation(LES) of fully developed channel flow. The results are compared with experimental data and LES results in terms of the mean statistics and temporal physics of the flow.

Numerical simulations of two-dimensional steady incompressible lid-driven flow in a square cavity are presented to verify the validity of a new solution code(PowerCFD) with unstructured grids. The code uses the non-staggered(collocated) grid approach which is very popular for incompressible flow analysis because of its numerical efficiency on the curvilinear or unstructured grids. Solutions are obtained for configurations with a Reynolds number as high as 10,000 with both rectangular and hybrid types of unstructured grid mesh. Interesting features of the flow are presented in detail and comparisons are made with benchmark solutions found in the literature. It is found that the code is capable of producing accurately the nature of the lid-driven cavity flow at high Reynolds numbers.

Noise generated by the blunt trailing edge of lifting surfaces is investigated in this study using fluid structure interaction theory. First, through the eddy modeling, noise generation doe to the flow instability on the rigid trailing edge is surveyed. Then the behavior of elastic cantileverd beam is investigated. Parametric study based on various material properties is employed to analyze the motion of the beam. Moreover, each eigenmode approach of cantilevered beam is used to find when flow induced vibration is resonant. To analyze elastic behavior of cantilever beam efficiently, moving grid generation technique based on non-conservative form of Navier-Stokes equation is used. Equation of the motion associated with the cantilever beam is discretized by the Galerkin procedure with forced vibration. As a consequence, behavior of the elastic cantilevered beam is stable when the first mode natural frequency of the material is relatively higher than that of flow induced pressure fluctuation.

A numerical analysis was made In investigate the pressure diminution of a silencer with baffles for high pressure blast flow fields. Reynolds-Averaged Navier-Stokes equations were solved for an axisymmetric computational domain constructed by multi block Chimera grids. A blast flow field was calculated for the silencers that are with one and three baffles. The evolution of high pressure blast flow fields was observed by depicting calculated contours of pressure. It was found that the tested silencer could achieve 97.7 percent pressure diminution.

A fluid transient analysis for the propellant flow in a monopropellant propulsion system is conducted using the method of characteristics (MOC). Algebraic simultaneous equations method and Clamor's rule method utilized to drive the compatible and characteristic equations are reviewed to understand MOC more extensively. The identification of fluid transient phenomena of propulsion system of Koreasat 1 is carried out through parametric studies. Also this work describes the reason that the propulsion system of Koreasat 1 has no orifice to control flow transients or to limit the initial hydrazine flow rate for the first-pulse firing.

In this study, the numerical and dynamic simulation on staging problem including forward jet mechanism is conducted. The forward jet plays a vital role in staging, which jets out from aftbody. This staging environment needs full dynamic characteristics study and flow analysis for securing staging safety. Present study performs dynamic simulation of prebody and aftbody with flow analysis using Chimera grid scheme which is usually used for moving simulations. As a result, separation mechanism using forward jet well works in staging for given initial conditions and reverse thrust, chamber pressure variation from experiments. Furthermore, it is found that the technique using forward jets for staging is excellent for securing the separation safety.

In this study, the reliablility based design optimization is peformed for an aircraft wing. The flexiblility of the wing was assumed by considering the interaction modeled by static aeroelasticity between aerodynamic forces and the structure. For a multidisciplinary design optimization the results of aerodynamic analysis and structural analysis were included in the optimization formulation. The First Order Reliability Method(FORM) was employed to consider the uncertainty of the designed points.

For the improvement of aerodynamic performance of the turbine blade in a turbopump for the liquid rocket engine, the optimization of turbine profile shape has been studied. The turbine in a turbopump in this study is a partial admission of impulse type, which has twelve nozzles and supersonic inflow. Due to the separated nozzles and supersonic expansion, the flow field becomes complicates and shows oblique shocks and flow separation. To increase the blade power, redesign of the blade shape using CFD and optimization method was attempted. The turbine cascade shape was represented by four design parameters. For optimization, genetic algorithm based upon non-gradient search has been selected as a optimizer. As a result, the final blade has about 4 percent more blade power than the initial shape.

According to operating conditions of each engine, a PCV valve has various flow rates and pressure characteristic. In a developed country, it has been developing by a computational design simulation. But, Korean companies have no ability of technical design for a PCV valve. So, they depend on their experiments and copy the designs of foreign companies when they need to design new PCV valves. These problems cause increase of error rate and take much time. Hence, optimal design for a PCV valve is needed to secure for continuous competition against foreign automobile companies. In this study, we used 4th order Runge-Kutta method for the prediction of spool movements and applied Bernoulli's equation for the determination of flow area. A spool geometry and spool displacement were predicted to be satisfied in comparison with their experiment.

A new numerical scheme for two-phase flows, the Hybrid VOF method has been developed for improved free surface capturing. The present new method is a volume capturing based VOF method coupled with a reinitialization procedure of a Level-set method. For validation, the proposed method is applied to two test cases: spherical bubble rising and dam breaking. The calculated results by using the Hybrid VOF method with the two previously applied VOF formulations are compared with available numerical and experimental data. It is found that the new method provides more accurate results than the two previous ones.

In this paper, we present the theoretical and numerical results of flow characteristics of a gas in a swirl injector. Proposed in this study are one-dimensional (theoretical) model and 2D/3D CFD models for use in the design of the injector. It was found that contrary to the classical theory about the compressible flow, the swirl gives a significant effect on the mass flow rate and the choking conditions. The one-dimensional model was found to Provide reasonably accurate results compared with the 2D/3D numerical results, so that it can be employed in th initial stage of the swirl-injector design process.

The purpose of this study is the characteristics of an innovative inlet system with shock/boundary layer interactions by using various types of bumps which are substituted for the conventional bleeding system in supersonic inlet. This study performs a comprehensive numerical effort that be directed at better understanding the three-dimensional flowfield includes shock/boundary layer interaction and growth of turbulent boundary layer that occur around a three-dimensional bump in a supersonic inlet. The characteristics of boundary layer seen in the current numerical simulations indicates the potential capability of the three-dimensional bump to control shock/boundary layer interaction in supersonic inlets.

A compressible two-fluid two-phase flow computation model using the stiffened-gas equation of state is formulated. Since the conservation equation system is of mixed type, it gives complex eigenvalues. The sonic speeds obtained from the individual single phase have been simply used in the literature for the fastest wave speeds necessary in the HLL scheme. This method has worked fine but proved to be quite diffusive according to our test. To improve the accuracy, we here propose to utilize the analytic eigenvalues evaluated from an approximate Jacobian matrix lot the fastest wave speeds. The interfacial transfer terms were dropped in constituting the Jacobian matrix for this purpose. The present scheme proved efficient, robust and accurate in comparison with other existing methods. We solved the cavitating flow problem using the present scheme. The result shows more detailed wave structure in the cavitating process caused by the strong expansion waves.

최근 전산 유체의 발달로 이상 유동해석의 캐비테이션 모델 적용 방법이 발전되어 왔으나 다양한 수력학적 시스템에서 발생하는 캐비테이션 유동은 난류이며 물과 공기 사이에서의 복잡한 상호 작용을 가지고 있으므로 그 적용 예가 아직은 미흡한 상태이다. 본 연구에서는 수중에서의 캐비테이션 해석과 이상 유동 해석을 위한 코드 개발 및 검증을 목적으로 3차원 회전체 주위의 캐비테이션 유동을 여러 가지 조건들의 변화를 적용하여 해석하였다. 또한 캐비테이션 발생과 관련한 다른 난류 모델에 적용하여 비교 분석을 수행하였다. 해석을 위한 모델의 지배방정식은 이상유동 Wavier-Stokes 방정식, 질량$\cdot$모멘텀 방정식의 혼합된 형태로 구성되어 있으며 방정식의 해를 구하기 위한 방법으로 유한차분법을 이용하였다. 해석결과의 신뢰성을 고려하여 반구형 실린더 주위의 캐비테이션 유동의 실험치와 비교 분석하였다. 그 결과, 본 연구의 수치 해석 방법과 실험적 방법의 결과가 강한 양의 상관관계를 가짐을 알 수 있었으며, 이러한 수치적 뒷받침은 본 연구의 전산수치 해석 방법이 앞으로의 여러 유동 해석으로의 적용 가능성을 보여준다.

Two-dimensional Stokes flows through a micro channel with a couple of symmetric vertical fins are investigated. At far up- and down-stream from the fins, the plane Poiseuille flow exists in the channel. The slip boundary conditions are applied to take account of the Knudsen number effects. For the analysis, the method of eigen function expansion and collocation method are employed. By the results, the streamline patterns and pressure distributions are shown and the force exerted on the fin and the excess pressure drop due to the fins are determined as functions of the length of the fin and Knudsen number. It may be conjectured that the force and the excess pressure drop are almost independent of the Knudsen number.