• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.1 s.244
• /
• pp.1-7
• /
• 2006

The characteristics of micro gaseous flows in microchannels have been analyzed in view of flow resistance using the direct simulation Monte Carlo (DSMC) method which is a molecule-based numerical modeling technique. For this purpose, a DSMC code where the pressure boundary condition was specified at the inlet and outlet, has been developed and the results of simulations showed satisfactory agreements with the analytic solution in the slip flow regime. (0.01 < Kn < 0.1) By varying the height and length of the microchannel, the effect of pressure difference between the inlet and outlet was examined. The present computation indicates that the curvature in pressure distribution along the channel increases due to the effect of compressibility when the pressure difference increases. To obtain the flow resistance regardless of the channel dimensions, a standard curve is devised in the present study by introducing the concept of unit mass flowrate and unit driving pressure force. From this curve, it is shown that in micro flows, a significant deviation from the laminar incompressible flow occurs by reducing the flow resistance.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.6
• /
• pp.860-869
• /
• 2000

Numerical and experimental investigations are performed for the molecular transition and slip flows in pumping channels of a disk-type drag pump. The flow occurring in the pumping channel develops from the molecular transition to the slip flow traveling downstream. Two different numerical methods are used in this analysis: the first one is a continuum approach in solving the Navier-Stokes equations with slip boundary conditions, and the second one is a stochastic approach through the use of the direct simulation Monte Carlo method. In the experimental study, the inlet pressures are measured for various outlet pressures in the range of 0.1{\sim}4Torr. From the present study, the numerical results of predicting the performance, obtained by both methods, agree well with the experimental data for the range of Knudsen number $Kn{\leq}0.1$ (i.e., the slip flow regime). But the results from the second method only agree with the experimental data for Kn>0.1(i.e., the molecular transition regime)

• Journal of computational fluids engineering
• /
• v.19 no.2
• /
• pp.40-48
• /
• 2014

For rarefied gas flow regimes, physical phenomena such as velocity slip and temperature jump occur on the solid body surface. To predict these phenomena accurately, either the Navier-Stokes solver with a slip boundary condition or the direct simulation Monte Carlo method should be used. In the present study, flow simulations of a wedge were conducted in Mach-10 flow of argon gas for several different flow regimes using a two-dimensional Navier-Stokes solver with the Maxwell slip boundary condition. The results of the simulations were compared with those of the direct simulation Monte Carlo method to assess the present method. It was found that the values of the velocity slip and the temperature jump predicted increase as the Knudsen number increases. Also, the results are comparatively reasonable up to the Knudsen number of 0.05.

• Journal of Power System Engineering
• /
• v.18 no.6
• /
• pp.207-214
• /
• 2014

Mass and flow resistance in a square ribbed microchannel have been studied numerically using the Lattice Boltzmann Method. It has been build up on two dimensional nine velocity vectors model with single relaxation time method called the Lattice Bhatnagor-Gross-Krook model. To analyze the roughness effect on the flow resistance namely the friction factor and mass flow has been discussed at the slip flow regime, $0.01{\leq}Kn{\leq}0.10$, where Kn is the Knudsen number. The wall roughness is considered by square microelements with a relative roughness height up to maximum 10% of channel height. The velocity profiles in terms of streamlines near the riblets are demonstrated to be responsible for the roughness effect. It is found that the roughness effect leads to increase the flow resistance with roughness height but it is decreased significantly with increasing the space between two roughness elements as well as the Knudsen number. In addition, the mass flow decreased linearly with increasing both roughness height and gap but significantly changed at the slip flow regime.

• Computers and Concrete
• /
• v.4 no.5
• /
• pp.377-402
• /
• 2007

The paper presents quasi-static plane strain FE-simulations of strain localization in reinforced concrete beams without stirrups. The material was modeled with two different isotropic continuum crack models: an elasto-plastic and a damage one. In case of elasto-plasticity, linear Drucker-Prager criterion with a non-associated flow rule was defined in the compressive regime and a Rankine criterion with an associated flow rule was adopted in the tensile regime. In the case of a damage model, the degradation of the material due to micro-cracking was described with a single scalar damage parameter. To ensure the mesh-independence and to capture size effects, both criteria were enhanced in a softening regime by nonlocal terms. Thus, a characteristic length of micro-structure was included. The effect of a characteristic length, reinforcement ratio, bond-slip stiffness, fracture energy and beam size on strain localization was investigated. The numerical results with reinforced concrete beams were quantitatively compared with corresponding laboratory tests by Walraven (1978).

• Journal of Advanced Marine Engineering and Technology
• /
• v.31 no.3
• /
• pp.272-281
• /
• 2007

A numerical simulation on nitrogen gas flow in a long parallel plate microchannel was performed to obtain the effect of compressibility and rarefaction on gaseous flow in microchannels. The simulation was based on steady. two dimensional compressible Navier-Stokes and energy equations with noslip and first order slip boundary conditions. The channel was $1.2{\mu}m$ deep and $3000{\mu}m$ long. The Reynolds numbers were in the range of order from $10^{-2}$ to $10^{-1}$. So the flow was assumed to be laminar. The computations were performed on various pressure ratios. The outlet pressure was fixed to atmospheric pressure. The outlet Knudsen number was 0.0585, consequently the flow was in the slip flow regime. The computations were performed with the assumption of isothermal channel walls. The results were compared with the experimental data. The agreement was good.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.4
• /
• pp.509-517
• /
• 1997

Results of flat plate compressible boundary layer calculation, based on discrete formulation of DSMC method, are presented in low Mach number and low Knudsen number range. The free stream is a uniform flow of pure nitrogen at various Mach numbers in low pressures (i.e. rarefied gas). Complete thermal accommodation and diffuse molecular reflections are used as the wall boundary condition, replacing unreal no-slip condition used in continuum calculations. In the discrete formulation of DSMC method, there is no need to use ad hoc assumptions on transport properties like viscosity and thermal conductivity, instead viscosity is calculated from values of other field variables (velocity and shear stress). Also the results are compared with existing self-similar continuum solutions. In all Mach number cases computed, velocity slip is most pronounced in regions near the leading edge where continuum formulation renders the solution singular. As the boundary layer develops further downstream, velocity slips asymptote to values that are between 10 to 20% of the magnitude of free stream velocity. When the free stream number density is reduced, so the gas more rarefied, the velocity slip increases as expected.

• 유체기계공업학회:학술대회논문집
• /
• 2006.08a
• /
• pp.391-394
• /
• 2006

In this paper, the thermal lattice Boltzmann method(TLBM) proposed by Guo et al.(2002) is applied to analyze the forced convective flow and heat transfer of 2-D micro channel. Nonequilibrium extrapolation boundary condition is adopted to simulate the velocity and temperature behavior at wall boundaries. Numerical results obtained by the present study give a good prediction of the micro fluidic characteristics with thermal effects.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.430-435
• /
• 2001

This paper introduces a pressure correction method for microflow computation. Conventional CFD methods with no slip boundary condition fail to predict the rarefaction effect of the wall when simulating gas microflows in the slip-flow regime. Pressure correction method with an appropriate slip boundary condition is an efficient tool in analyzing microscale flows. The present unstructured SIMPLE algorithm adopts both the classical Maxwell boundary condition and Langmuir boundary condition proposed by Myong. The simulation results of microchannel flows show that the proposed method has an effective predictive capability for microscale flows.

• KSTLE International Journal
• /
• v.7 no.2
• /
• pp.27-34
• /
• 2006

Knudsen number is the ratio of molecular mean free path versus mm thickness and the criterion to determine the flow form. When its value is lower than 0.01, the flow can be assumed to has no slip boundary condition. And in the case that the value is between 0.01 and 10, then the flow has slip boundary condition at both the adjacent walls. The condition of the air flow between the rotating journal and top foil in the air foil bearing is determined by the rotating speed and load, and the Knudsen number is also varied by those values. Because the molecular mean free path is variable to the pressure and temperature, more exact formulation is necessary to understand and analyze the flow regime. In this study, the analysis considering Knudsen number formulated with those variables (pressure, temperature and mm thickness) was executed. The approximate value was examined using the equation to confirm whether the flow has the slip or no-slip boundary condition. From the analytic investigation, it was decided to range approximately 0.01 to 1.0 and the flow can be supposed to have the slip boundary condition. Under the condition of the slip flow, the static characteristics of the air foil bearing were examined using modified Reynolds equations. The results were compared with those considering no slip condition. It shows that the slip condition makes the flow decelerates and the load carrying capacity decreases compared with no slip condition. And as the bearing number and eccentricity ratio increase, the load carrying capacity also increased at both the cases. From this result, it can be supposed that the bearing torque also increases. In the analysis of the dynamic characteristics, the perturbed Knudsen number was taken into consideration. Because the Knudsen number is expressed as the terms of each variable, the perturbed equation can be simply derived. The results of both cases considering and not considering Knudsen number were compared each other. In the case of the direct terms of the stiffness and damping coefficients, the difference between both cases was little and increased as the bearing number and eccentricity ratio increased. And the cross terms have less or more differences.