• 한국전산유체공학회지
• /
• 제5권1호
• /
• pp.1-7
• /
• 2000

In general, FDM(finite difference method) and FVM(finite volume method) are used for analyzing the fluid flow numerically. However it is difficult to apply them to problems involving complex geometries, multi-connected domains, and complex boundary conditions. On the contrary, FEM(finite element method) with coordinates transformation for the unstructured grid is effective for the complex geometries. Most of previous studies have used commercial codes such as KIVA or STAR-CD for the flow analyses in the engine cylinder, and these codes are mostly based on the FVM. In the present study, using the FEM for three-dimensional, unsteady, and incompressible Navier-Stokes equation, the velocity and pressure fields in the engine cylinder have been numerically analyzed. As a numerical algorithm, 4-step time-splitting method is used and ALE(arbitrary Lagrangian Eulerian) method is adopted for moving grids. In the Piston-Cylinder, the calculated results show good agreement in comparison with those by the FVM and the experimental results by the LDA.

• 대한기계학회논문집B
• /
• 제25권12호
• /
• pp.1802-1812
• /
• 2001

The purpose of this study is to analyze numerically the movement of particles included in turbulent fluid flow characteristics of metallic surfaces. To describe fluid flew, the incompressible Navier-Stokes equation discretized by the finite volume method were solved on the non-orthogonal coordinates with non-staggered variable arrangement, and the k-$\xi$ turbulence model was adapted. After fluid flow was calculated, particle movement was predicted from the Lagrangian approaches. Non-essential complexities were avoided by assuming that the particles had spherical shapes and the Stoke's drag formula only consisted of external farces acting upon them. In order to validate the numerical calculations, the results were compared with the experimental data reported in literature and agreed well with them. The drag force coefficient equation showed better agreement with the experimental data in the prediction of particle movement than the correction factor equation. Impact velocity and impact angle increased as inlet turbulence intensity decreased, relative jet height was lower. or the Reynolds number was larger.

• 한국소음진동공학회논문집
• /
• 제17권1호
• /
• pp.41-47
• /
• 2007

Scattering fields of two dimensional acoustic waves by a circular cylinder are investigated. The present numerical approach for the acoustic scattering problem has difficulties of numerical robustness, long-time stability and suitability of far-field boundary treatments. The time-dependent periodic acoustic source is used to analyze Interference patterns between incident waves and waves reflected by the cylinder. Characteristic boundary algorithms coupled with 4th order Modified-Flux-Approach ENO(essentially non-oscillatory) schemes are employed in generalized coordinates to examine the effect of the wane frequency on the interference patterns. Non-reflecting boundary conditions, which is crustal for accurate computations of aeroacoustic problems, are used not to contaminate scattering fields by reflected waves at the outer boundary. Computed scattering fields show the circumferential acoustic modes generated by interacting between acoustic sources and scattered waves. At a lower frequency, the wave passes almost straight through the cylinder without Interacting with circular cylinder. Simulation results are presented and compared with the analytic solution. Computed RMS-pressure distribution on the cylinder wall is good agreement with exact solution.

• Journal of the Optical Society of Korea
• /
• 제17권6호
• /
• pp.538-542
• /
• 2013

We propose a simple model to elucidate the dispersion behavior of spiraling modes on silver nanowire by finding correspondence parameters and building a simple equivalent relationship with the planar insulator-metal-insulator geometry. The characteristics approximated for the proposed structure are compared with the results from an exact solution obtained by solving Maxwell's equation in cylindrical coordinates. The effective refractive index for our proposed equivalent model is in good agreement with that for the exact solution in the 400-2000 nm wavelength range. In particular, when the radius of the silver nanowire is 100 nm, the calculated index shows typical improvements; the average percentage error for the real part of the effective refractive index is reduced to only 5% for the $0^{th}$ order mode (11.9% in previous results) and 1.5% for the $1^{st}$ order mode (24.8% in previous results) in the 400-800 nm wavelength range. This equivalent model approach is expected to provide further insight into understanding the important behavior of nanowire waveguides.

• 한국공간구조학회논문집
• /
• 제5권4호
• /
• pp.79-90
• /
• 2005

본 연구의 주 목적은 다양한 요인에 의해 케이블-막구조물에 요소이동이 발생할 때, 당초 해석 시 가정되었던 좌표나 응력상태의 변화에 대해 막과 케이블 사이에 발생하는 요소이동의 방향이나 크기를 산정하고, 요소이동이 발생한 후 응력상태의 변화를 규명하는 것이다. 먼저 케이블 보강 막구조물의 요소이동 문제를 해석하기 위한 이론적 배경인 ALE 유한요소법의 개념을 소개하고, ALE 개념이 도입된 케이블-막구조물에서의 요소이동을 고려한 강성매트릭스를 작성하여 해석 프로그램을 개발한다. 개발된 프로그램의 타당성을 검증하기 위해 다양한 예제 해석을 수행한다.

• Nuclear Engineering and Technology
• /
• 제41권10호
• /
• pp.1347-1360
• /
• 2009

A multi-dimensional component for the thermal-hydraulic system analysis code, MARS, was developed for a more realistic three-dimensional analysis of nuclear systems. A three-dimensional and two-fluid model for a two-phase flow in Cartesian and cylindrical coordinates was employed. The governing equations and physical constitutive relationships were extended from those of a one-dimensional version. The numerical solution method adopted a semi-implicit and finite-difference method based on a staggered-grid mesh and a donor-cell scheme. The relevant length scale was very coarse compared to commercial computational fluid dynamics tools. Thus a simple Prandtl's mixing length turbulence model was applied to interpret the turbulent induced momentum and energy diffusivity. Non drag interfacial forces were not considered as in the general nuclear system codes. Several conceptual cases with analytic solutions were chosen and analyzed to assess the fundamental terms. RPI air-water and UPTF 7 tests were simulated and compared to the experimental data. The simulation results for the RPI air-water two-phase flow experiment showed good agreement with the measured void fraction. The simulation results for the UPTF downcomer test 7 were compared to the experiment data and the results from other multi-dimensional system codes for the ECC delivery flow.

• JSTS:Journal of Semiconductor Technology and Science
• /
• 제11권4호
• /
• pp.278-286
• /
• 2011

In this paper, we present a compact model of gate-voltage-dependent quantum effects in short-channel surrounding-gate (SG) metal-oxide-semiconductor field-effect transistors (MOSFETs). We based the model on a two-dimensional (2-D) analytical solution of Poisson's equation using cylindrical coordinates. We used the model to investigate the electrostatic potential and current sensitivities of various gate lengths ($L_g$) and radii (R). Schr$\ddot{o}$dinger's equation was solved analytically for a one-dimensional (1-D) quantum well to include quantum effects in the model. The model takes into account quantum effects in the inversion region of the SG MOSFET using a triangular well. We show that the new model is in excellent agreement with the device simulation results in all regions of operation.

• Structural Engineering and Mechanics
• /
• 제13권2호
• /
• pp.135-154
• /
• 2002

This paper presents the convected material frame approach to study the nonlinear behavior of inelastic frame structures. The convected material frame approach is a modification of the co-rotational approximation by incorporating an adaptive convected material frame in the basic definition of the displacement vector and strain tensor. In the formulation, each discrete element is associated with a local coordinate system that rotates and translates with the element. For each load increment, the corresponding strain-displacement and nodal force-stress relationships are defined in the updated local coordinates, and based on the updated element geometry. The rigid body motion and deformation displacements are decoupled for each increment. This modified approach incorporates the geometrical nonlinearities through the continuous updating of the material frame geometry. A generalized nonlinear function is used to derive the inelastic constitutive relation and the kinematic hardening is considered. The equation of motion is integrated by an explicit procedure and it involves only vector assemblage and vector storage in the analysis by assuming a lumped mass matrix of diagonal form. Several numerical examples are demonstrated in close agreement with the solutions obtained by the ANSYS code. Numerical studies show that the proposed approach is capable of investigating large deflection of inelastic planar structures and providing an excellent numerical performance.

• Bulletin of the Korean Chemical Society
• /
• 제13권3호
• /
• pp.317-324
• /
• 1992

Structure and dynamics of $Na^+$ ions are investigated by molecular dynamics simulations of rigid dehydrated zeolite-A at several temperatures using a simple Lennard-Jones potential plus Coulomb potential. A best-fitted set of electrostatic charges is chosen from the results of simulation at 298.15 K and Ewald summation technique is used for the long-ranged character of Coulomb interaction. The calculated x, y, and z coordinates of $Na^+$ ions are in good agreement with the positions determined by X-ray crystallography within statistical errors, their random movings in different types of closed cages are well described by time-correlation functions, and $Na_Ⅰ$ type ions are found to be less diffusive than $Na_Ⅱ$ and $Na_{III}$. At 600.0 K, the unstable $Na_{III}$ type ion pushes down one of nearest $Na_{I}$ ions into the $\beta-cage$ and sits on the stable site Ⅰ, and the captured ion in the $\beta-cage$ wanders over and attacks one of 8 $Na_{I}$ type ions.

• Steel and Composite Structures
• /
• 제38권2호
• /
• pp.189-211
• /
• 2021

An analytical solution is presented to analyze the thermoelastoplastic response of a rotating thick-walled cylindrical pressure vessel made of functionally graded material (FGM). The analysis is based on Tresca's yield condition, its associated flow rule and linear strain hardening material behaviour. The uncoupled theory of thermoelasticity is used, and the plane strain condition is assumed. The material properties except for Poisson's ratio, are assumed to vary nonlinearly in the radial direction. Elastic, partially plastic, fully plastic, and residual stress states are investigated. The heat conduction equation for the one-dimensional problem in cylindrical coordinates is used to obtain temperature distribution in the vessel. It is assumed that the inner surface is exposed to an airstream and that the outer surface is exposed to a uniform heat flux. Tresca's yield criterion and its associated flow rule are used to formulate six different plastic regions for a linearly hardening condition. All these stages are studied in detail. It is shown that the thermoelastoplastic stress response of a rotating FGM pressure vessel is affected significantly by the nonhomogeneity of the material and temperature gradient. The results are validated with those of other researchers for appropriate values of the system parameters and excellent agreement is observed.