• Wind and Structures
• /
• v.19 no.6
• /
• pp.687-708
• /
• 2014

The performance of the $k-{\varepsilon}$ and $k-{\omega}$ two-equation turbulence models was investigated in computational simulations of the neutrally stratified atmospheric boundary layer developing above various terrain types. This was achieved by using a proposed methodology that mimics the experimental setup in the boundary layer wind tunnel and accounts for a decrease in turbulence parameters with height, as observed in the atmosphere. An important feature of this approach is pressure regulation along the computational domain that is additionally supported by the nearly constant turbulent kinetic energy to Reynolds shear stress ratio at all heights. In addition to the mean velocity and turbulent kinetic energy commonly simulated in previous relevant studies, this approach focuses on the appropriate prediction of Reynolds shear stress as well. The computational results agree very well with experimental results. In particular, the difference between the calculated and measured mean velocity, turbulent kinetic energy and Reynolds shear stress profiles is less than ${\pm}10%$ in most parts of the computational domain.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.2
• /
• pp.141-149
• /
• 2005

This paper deals with formulations of the energy equilibrium equation by an introduction of the algebraic description, quarternion, which meets conservations of system energy for the equation of motion. Then the equation is discretized to analyze the dynamits analysis of flexible multibody systems in such a way that the work done by the constrained force completely is eliminated. Meanwhile, Rodrigues parameters we used to express the finite rotation lot the proposed method. This method lot the initial essential step to a guarantee of developments of the 3D dynamical problem provides unconditionally stable conditions for the nonlinear problems through the numerical examples.

• The Transactions of the Korean Institute of Electrical Engineers P
• /
• v.61 no.1
• /
• pp.13-17
• /
• 2012

Ionization and Attachment Coefficients 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-term approximation of the Boltzmann equation (BEq.) method and the Monte Carlo simulation (MCS). The calculations of electron swarm parameters require the knowledge of several collision cross-sections of electron beam. Thus, published momentum transfer, ionization, vibration, attachment, electronic excitation, and dissociation cross-sections of electrons for $CH_4$, $CF_4$ and Ar, were used. 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. The proposed theoretical simulation techniques in this work will be useful to predict the fundamental process of charged particles and the breakdown properties of gas mixtures.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1997.10a
• /
• pp.599-604
• /
• 1997

One of important inconsistencies in the six-equation model predictions has been found to be the force experienced by a single bubble placed in a convergent stream of liquid. Various sets of governing equations yield different amount of forces to hold the bubble stationary in a convergent nozzle. By using the first order potential flow theory, it is found that the six-equation model can not be used to estimate the force experienced by a deformed bubble. The theoretical value of the particle stress of a bubble in a convergent nozzle flow has been found to be a function of the Weber number when bubble distortion is allowed. This force has been calculated by using different sets of governing equations and compared with the theoretical value. It is suggested in this study that the bubble size distribution function can be used to remove the presented inconsistency by relating the interfacial variables with different moments of the bubble size distribution function. This study also shows that the inconsistencies in the thermal-hydraulic governing equation can be removed by mechanistic modeling of the phasic interface.

• The Proceedings of the Korean Institute of Illuminating and Electrical Installation Engineers
• /
• v.11 no.6
• /
• pp.37-43
• /
• 1997

Density of a dye or pigment teduces by radiant energym and dinsity change induces color change. In this paper, we habe worked out an equation to predict the color change by the radiant energy under the assumption of colorant density changing pattern by radiant energy and the color changing pattern by colorant density chane. And the damaging factors in the equation was established by experiments. To evaluate the efectiveness of the equation and the method, several samples were irradiated by common light source and color change was measured. Good coincidence between the ceasured and the calculated color change was found. And it was also found that there is a close relationhaip between the spectral damaging factor and the spectral absorptance in visible wavelength range.

• The Transactions of the Korean Institute of Electrical Engineers P
• /
• v.55 no.2
• /
• pp.57-61
• /
• 2006

This paper calculates and gives the analysis of mean energy in pure $SiH_4,\;Ar-SiH_4$ mixture gases ($SiH_4-0.5[%],\;5[%]$) over the range of $E/N =0.01{\sim}300[Td]$, p = 0.1, 1, 5.0 [Torr] by Monte Carlo the Backward prolongation method of the Boltzmann equation using computer simulation without using expensive equipment. The results have been obtained by using the electron collision cross sections by TOF, PT, SST sampling, compared with the experimental data determined by the other author. It also proved the reliability of the electron collision cross sections and shows the practical values of computer simulation. The calculations of electron swarm parameters require the knowledge of several collision cross-sections of electron beam. Thus, published momentum transfer, ionization, vibration, attachment, electronic excitation, and dissociation cross-sections of electrons for $SiH_4$ and Ar, were used. The differences of the transport coefficients of electrons in $SiH_4$, mixtures of $SiH_4$ and Ar, have been explained by the deduced energy distribution functions for electrons and the complete collision cross-sections for electrons. A two-term approximation of the Boltzmann equation analysis and Monte Carlo simulation have been used to study electron transport coefficients.

• 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).

• Structural Engineering and Mechanics
• /
• v.70 no.6
• /
• pp.737-750
• /
• 2019

This paper investigates the static and dynamic behaviors of imperfect single walled carbon nanotube (SWCNT) modeled as a beam structure by using energy-equivalent model (EEM), for the first time. Based on EEM Young's modulus and Poisson's ratio for zigzag (n, 0), and armchair (n, n) carbon nanotubes (CNTs) are presented as functions of orientation and force constants. Nonlinear Euler-Bernoulli assumptions are proposed considering mid-plane stretching to exhibit a large deformation and a small strain. To simulate the interaction of CNTs with the surrounding elastic medium, nonlinear elastic foundation with cubic nonlinearity and shearing layer are employed. The equation governed the motion of curved CNTs is a nonlinear integropartial-differential equation. It is derived in terms of only the lateral displacement. The nonlinear integro-differential equation that governs the buckling of CNT is numerically solved using the differential integral quadrature method (DIQM) and Newton's method. The linear vibration problem around the static configurations is discretized using DIQM and then is solved as a linear eigenvalue problem. Numerical results are depicted to illustrate the influence of chirality angle and imperfection amplitude on static response, buckling load and dynamic behaviors of armchair and zigzag CNTs. Both, clamped-clamped (C-C) and simply supported (SS-SS) boundary conditions are examined. This model is helpful especially in mechanical design of NEMS manufactured from CNTs.

• Nuclear Engineering and Technology
• /
• v.53 no.7
• /
• pp.2312-2322
• /
• 2021

In recent studies, the creep rate of Zircaloy-4, one of the basic property parameters of the nuclear fuel code, has been commonly used with the axial creep model proposed by Rosinger et al. However, in order to calculate the circumferential deformation of the fuel cladding, there is a limitation that a difference occurs depending on the anisotropic coefficients used in deriving the circumferential creep equation by using the axial creep equation. Therefore, in this study, the existing axial creep law and the derived circumferential creep results were analyzed through a circumferential creep test by the internal pressurization method in the isothermal conditions. The circumferential creep deformation was measured through the optical image analysis method, and the results of the experiment were investigated through constructed IDECA (In-situ DEformation Calculation Algorithm based on creep) code. First, preliminary tests were performed in the isotropic β-phase. Subsequently in the anisotropic α-phase, the correlations obtained from a series of circumferential creep tests were compared with the axial creep equation, and optimized anisotropic coefficients were proposed based on the performed circumferential creep results. Finally, the IDECA prediction results using optimized anisotropic coefficients based on creep tests were validated through tube burst tests in transient conditions.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.56 no.2
• /
• pp.367-373
• /
• 2007

The present study investigates numerically nonequilibrium energy transfer between electrons and phonons in metal thin films irradiated by ultrashort pulse lasers and it also provides the temporal and spatial variations of electron and phonon temperatures using the well-established two-temperature model(TTM) on the basis of the Boltzmann transport equation(BTE). This article predicts the crater shapes in gold film structures, and compares the results by using two-dimensional energy transport equation. From the results, it is found that nonequilibrium energy transfer between electrons and phonons takes place, and the equilibrium time increases with the increase of laser fluence. On the other hand, above threshold fluence the ablation time doesn't change nearly with increasing fluences. Compared with one-dimensional TTM, it also reveals that the temporal distributions of electron and phonon temperatures at the top surface estimated by using two-dimensional TTM have a similar tendency. The results show that two-dimensional TTM can simulate the crater shape of metals during the irradiation of femtosecond pulse lasers and the absorbed energy is propagated to z-direction faster than to r-direction.