• Journal of the Institute of Electronics Engineers of Korea SD
• /
• v.39 no.2
• /
• pp.14-26
• /
• 2002

A new model for degenerate semiconductor quantum well devices was developed. In this model, the multi-subband Boltzmann transport equation was formulated by applying the Pauli exclusion principle and coupled to the Schrodinger and Poisson equations. For the solution of the resulted nonlinear system, the finite difference method and the Newton-Raphson method was used and carrier energy distribution function was obtained for each subband. The model was applied to a Si MOSFET inversion layer. The results of the simulation showed the changes of the distribution function from Boltzmann like to Fermi-Dirac like depending on the electron density in the quantum well, which presents the appropriateness of this modeling, the effectiveness of the solution method, and the importance of the Pauli -exclusion principle according to the reduced size of semiconductor devices.

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

Surface roughness is present in most of the microfluidic devices due to the microfabrication techniques. This paper presents lattice Boltzmann method (LBM) results for laminar flow in a microchannel with surface roughness. The surface roughness is modeled by an array of rectangular modules placed on top and bottom side of a parallel-plate channel. In this study, LBGK D2Q9 code in lattice Boltzmann Method is used to simulate flow field for low Reynolds number in a micro-channel. The effects of relative surface roughness, roughness distribution, roughness size and the results are presented in the form of the product of friction factor and Reynolds number. Finally, a significant increase in Poiseuille number is detected as the surface roughness is considered, while the effect of roughness on the microflow field depends on the surface roughness.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.54-57
• /
• 2011

In this paper we present an algorithm about how to simulate two dimensional dam breaking with lattice Boltzmann method (LBM). LBM considers a typical volume element of fluid to be composed of a collection of particles that represented by a particle velocity distribution function for each fluid component at each grid point. We use the modified Lattice Boltzmann Method for incompressible fluid. This paper will represent detailed information on single phase flow which considers only the water instead of both air and water. Interface treatment and conservation of mass are the most important things in simulating free surface where the Interface is treated by mass exchange with the water region. We consider the surface tension on the interface and also bounce back boundary condition for the treatment of solid obstacles. We will compare the results of the simulation with some methods and experimental results.

• The Transactions of the Korean Institute of Electrical Engineers P
• /
• v.60 no.1
• /
• pp.27-31
• /
• 2011

In this paper, the electron transport characteristics in $CF_4$ has been analysed over the E/N range 1~300[Td] by a two-term approximation Boltzmann equation method and by a Monte Carlo simulation. The motion has been calculated to give swarm parameters for the electron drift velocity, longitudinal diffusion coefficient, the ratio of the diffusion coefficient to the mobility, electron ionization and attachment coefficients, effective ionization coefficient, mean energy, collision frequency and the electron energy distribution function. The electron energy distribution function has been analysed in $CF_4$ at E/N=5, 10, 100, 200 and 300[Td] for a case of the equilibrium region in the mean electron energy and respective set of electron collision cross sections. The results of Boltzmann equation and Monte Carlo simulation have been compared with experimental data by Y. Nakamura and M. Hayashi. The swarm parameter from the swarm study are expected to serve as a critical test of current theories of low energy electron scattering by atoms and molecules, in particular, as well as crucial information for quantitative simulations of weakly ionized plasmas.

• Proceedings of the Korean Institute of Electrical and Electronic Material Engineers Conference
• /
• 1997.11a
• /
• pp.167-172
• /
• 1997

The electron swarm parameters and Energy distribution function have been calculated for electrons motion through CH$_4$ pure gas under the action of uniform electric field for 0.1$\leq$E/N(Td)$\leq$300, at the 300( $^{\circ}$K), using MCS method and Boltzmann transport equation. And then the resulting values of electron drift velocity were compared to experimental data and adjustment made in assumed cross sections until good agreement was obtained. The electron drift velocity is very useful in the fields of study relating to the conductive and dielectric phenomena of gas medium. The electron energy distribution in gas discharge are generally nonmaxwellian , and must be calculated by a numerical solution of the Boltzmann equation which takes in the elastic and inelastic collisions. To analyze the physical phenomena and properties (or electron swarm motion in a gas under the influence of an electric field, the energy distribution function of electrons and the theoretical deriveration of the electron drift velocity are calculated by the Backward Prolongation with respect to the Boltzmann transport equation as a parameter of E/N(Td).

• The KSFM Journal of Fluid Machinery
• /
• v.15 no.3
• /
• pp.12-18
• /
• 2012

The Lattice Boltzmann Method has developed for solving the Boltzmann equation in Cartesian domains containing immersed boundaries of arbitrary geometrical complexity moving with prescribed kinematics. When a immersed boundaries are sweeping the fixed fluid node, refilling the node information in a vicinity of fluid nodes is one of the important issues in Lattice Boltzmann Method. In this study, we propose a simple refill algorithm for the particle distribution function based on a proper velocity, density and strain rate to enhance accuracy and stability of the method. The refill scheme based on a asymptotic analysis of LBGK model has improved accuracy than interpolation schemes. The proposed scheme in this study is validated by the simulations of an impulsively started rotating circular cylinder to investigate adaptability for fluid-structure interaction (FSI) problem. This refill scheme has improved stability and accuracy especially at high Reynolds number region.

• Journal of the Korean Vacuum Society
• /
• v.10 no.4
• /
• pp.424-430
• /
• 2001

An electron Boltzmann equation is solved numerically to calculate the electron energy distribution functions in plasma discharge which is generated by radio-frequency (RF) and microwave frequency electric field. The maintenance field strengths are determined self-consistently by solving the homogeneous electron Boltzmann equation in the Lorentz approximation expressed by 2nd order differential equation and an additional particle balance equation expressed by integro-differential equation. By using this numerical code, the electron energy distribution functions in argon discharge are calculated in the range from RF to microwave frequency. The influence of frequency of the HF electric field on the electron energy distribution functions and ionization rate are investigated.

• Proceedings of the Korean Institute of Electrical and Electronic Material Engineers Conference
• /
• 2002.05c
• /
• pp.194-198
• /
• 2002

Energy Distribution Function in pure $CH_4$, $CF_4$ and mixtures of $CF_4$ and Ar, have been analyzed over a range of the reduced electric field strength between 0.1 and 350[Td] by the two-tenn approximation of the Boltzmann equation (BEq.) method and the Monte Carlo simulation (MCS). The results of the Boltzmann equation and the Monte Carlo simulation have been compared with the data presented by several workers. The deduced transport coefficients for electrons agree reasonably well with the experimental and simulation data obtained by Nakamura and Hayashi. The energy distribution function of electrons in $CF_4-Ar$ mixtures shows the Maxwellian distribution for energy. That is, f(${\varepsilon}$) has the symmetrical shape whose axis of symmetry is a most probably energy

• Proceedings of the KIEE Conference
• /
• 2004.07e
• /
• pp.37-40
• /
• 2004

Electron swarm parameters in pure $CF_4$ and mixtures of $CF_4$ and Ar, have been analyzed over a range of the reduced electric field strength between 0.1 and 350[Td] by the two-term approximation of the Boltzmann equation(BEq.) method and the Monte Carlo simulation(MCS) The results of the Boltzmann equation and the Monte Carlo simulation have been compared with the data presented by several workers. The deduced transport coefficients for electrons agree reasonably well with the experimental and simulation data obtained by Nakamura and Hayashi. The energy distribution function of electrons in $CF_4$-Ar mixtures shows the Maxwellian distribution for energy. That is, f(${\varepsilon}$) has the symmetrical shape whose axis of symmetry is a most probably energy

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2797-2802
• /
• 2007

It has been confirmed that implementation of the no-slip boundary conditions for the lattice-Boltzmann method play an important role in the overall accuracy of the numerical solutions as well as the stability of the solution procedure. We in this paper propose a new algorithm, i.e. the method of the dynamic boundary condition for no-slip boundary condition. The distribution functions on the wall along each of the links across the physical boundary are assumed to be composed of equilibrium and nonequilibrium parts which inherit the idea of Guo's extrapolation method. In the proposed algorithm, we apply a dynamic equation to reflect the computational slip velocity error occurred on the actual wall boundary to the correction; the calculated slip velocity error dynamically corrects the fictitious velocity on the wall nodes which are subsequently employed to the computation of equilibrium distribution functions on the wall nodes. Along with the dynamic selfcorrecting process, the calculation efficiently approaches the steady state. Numerical results show that the dynamic boundary method is featured with high accuracy and simplicity.