• 한국입자에어로졸학회지
• /
• 제5권3호
• /
• pp.123-131
• /
• 2009

This paper presents simulation results of particle transport in low-pressure, low-temperature plasma environment. The size dependent transport of particles in the plasma is investigated with a two-dimensional simulation tool developed in-house for plasma chamber analysis and design. The plasma model consists of the first two and three moments of the Boltzmann equation for ion and electron fluids respectively, coupled to Poisson's equation for the self-consistent electric field. The particle transport model takes into account all important factors, such as gravitational, electrostatic, ion drag, neutral drag and Brownian forces, affecting the motion of particles in the plasma environment. The particle transport model coupled with both neutral fluid and plasma models is simulated through a Lagrangian approach tracking the individual trajectory of each particle by taking a force balance on the particle. The size dependant trap locations of particles ranging from a few nm to a few ${\mu}m$ are identified in both electropositive and electronegative plasmas. The simulation results show that particles are trapped at locations where the forces acting on them balance. While fine particles tend to be trapped in the bulk, large particles accumulate near bottom sheath boundaries and around material interfaces, such as wafer and electrode edges where a sudden change in electric field occurs. Overall, small particles form a "dome" shape around the center of the plasma reactor and are also trapped in a "ring" near the radial sheath boundaries, while larger particles accumulate only in the "ring". These simulation results are qualitatively in good agreement with experimental observation.

• Ocean Systems Engineering
• /
• 제12권1호
• /
• pp.63-86
• /
• 2022

A refined projection-based purely Lagrangian meshfree method is presented towards reliable numerical analysis of fluid flow interactions with saturated/unsaturated porous media of uniform/spatially-varying porosities. The governing equations are reformulated on the basis of two-phase mixture theory with incorporation of volume fraction. These principal equations of mixture are discretized in the context of Incompressible SPH (Smoothed Particle Hydrodynamics) method. Associated with the consideration of governing equations of mixture, a new term arises in the source term of PPE (Poisson Pressure Equation), resulting in modified source term. The linear and nonlinear force terms are included in momentum equation to represent the resistance from porous media. Volume increase of fluid particles are taken into consideration on account of the presence of porous media, and hence multi-resolution ISPH framework is also incorporated. The stability and accuracy of the proposed method are thoroughly examined by reproducing several numerical examples including the interactions between fluid flow and saturated/unsaturated porous media of uniform/spatially-varying porosities. The method shows continuous pressure field, smooth variations of particle volumes and regular distributions of particles at the interface between fluid and porous media.

• 콘크리트학회논문집
• /
• 제19권6호
• /
• pp.803-810
• /
• 2007

FRP 합성재료로 구속된 콘크리트의 응력-변형률 응답을 합리적으로 예측할 수 있는 해석 모델이 제시되었다. 제안된 모델은 하중이 증가함에 따라 점진적으로 발생하는 미세균열에 의한 부피팽창이 미세 재료 구조의 손상을 나타내는 중요한 척도이며, 이에 손상 정도에 따라 하중 지지 능력을 일관되게 산정할 수 있다는 기본 개념에 근거한다. 이를 위하여 제안 모델은 면적 변형률 및 공극의 함수로 표시된 탄성계수, 팽창 콘크리트와 구속 매체의 상호작용을 나타내는 에너지 평형식, 변화하는 구속력 및 점증 계산 논리를 포함한다. 따라서 실험으로부터 유도된 팽창비 관계식으로부터 횡방향 혹은 부피팽창변형률을 산정하는 기존의 해석 모델과는 달리 역학적 거동 및 에너지 평형식으로부터 연속적으로 변화하는 횡방향 변형률을 산정한다. 구속된 콘크리트의 전체 응답을 예측할 수 있는 기존의 여러 해석 모델에 대하여 검토하였으며, 특히 부피 팽창을 고려하는 방법에 초점을 맞추어 토의하였다. 제안된 모델 및 기존 Samaan의 2선식 모델을 사용하여 실험 결과를 예측한 결과, 만족할 만한 범위 내에서 일치를 나타냈으나, Samaan의 2선식 모델은 부피 팽창 거동을 위하여 단지 초기 포아송비와 최종 수렴 팽창변화율 만을 고려하기 때문에 횡방향 변형률 응답을 예측하는 데는 한계가 있는 것으로 판단된다. 제안된 모델은 역학적 거동에 근거하여 다양한 관련 응답을 산정하므로 다른 분야에도 쉽게 적용할 수 있다.

• 태양에너지
• /
• 제17권3호
• /
• pp.59-66
• /
• 1997

본 연구에서는 PTV기법에 의해 속도를 계측하였으며 계측시스템은 레이저를 이용한 시이트라이트와 이미지보오드를 내장한 퍼스널컴퓨터 그리고 관련 소프트웨어로 구성하였다. 속도벡터는 PTV 시스템에 의해 구하였으며 그 결과를 압력에 관한 프와송방정식에 적용하여 전유동장의 압력을 구하였다. 이때의 경계조건으로서는 벽면에서는 부착조건을 부여하였고 구동류가 흐르는 캐비티의 상부에서는 내부의 속도 값을 외삽하였다. 압력계산에 있어서 격자의 해상도는 $40{\times}40$이며 등간격의 엇갈림격자를 사용하였다. 계측결과는 유동장의 속도 및 압력분포를 잘 나타내었다.

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

• 한국콘크리트학회:학술대회논문집
• /
• 한국콘크리트학회 2008년도 추계 학술발표회 제20권2호
• /
• pp.201-204
• /
• 2008

FRP 합성재료로 구속된 콘크리트의 응력-변형률 응답을 합리적으로 예측할 수 있는 해석모델이 제시되었다. 제안된 모델은 하중이 증가함에 따라 점진적으로 발생하는 미세균열에 의한 부피팽창이 미세재료구조의 손상을 나타내는 중요한 척도이며, 이에 손상정도에 따라 하중지지능력을 일관되게 산정할 수 있다는 기본개념에 근거한다. 이를 위하여 제안모델은 면적변형률 및 공극의 함수로 표시된 탄성계수, 팽창콘크리트와 구속매체의 상호작용을 나타내는 에너지 평형식, 변화하는 구속력 및 점증계산논리를 포함한다. 따라서 실험으로부터 유도된 팽창비 관계식으로부터 횡방향 혹은 부피팽창변형률을 산정하는 기존의 해석모델과는 달리 역학적 거동 및 에너지 평형식으로부터 연속적으로 변화하는 횡방향 변형률을 산정한다.

• 한국방재학회 논문집
• /
• 제10권1호
• /
• pp.111-117
• /
• 2010

수중구조물에 의한 파랑의 변형을 예측하기 위해 3차원 수치모형을 도입하여 수치모형 실험을 수행하였다. 본 수치모형은 Navier-Stokes 방정식을 유한차분법을 이용하여 계산하는 동수압 모형으로서, 난류의 해석을 위해서 상대적으로 큰 에디(eddy)만을 고려하는 SANS(Spatially Averaged Navier-Stokes) 방정식의 해를 구하는 LES(large-eddy-simulation) 기반의 수치모형이다. 엇갈림 격자체계에서 유한차분법을 사용하여 지배방정식을 해석하는 모형으로서 수치기법으로 Two-step projection 기법을 사용하여 SANS 방정식을 계산하였으며, Bi-CGSTAB 기법을 이용하여 Poisson 방정식의 해를 구하고 압력장을 계산하였다. 또한, 자유수면의 추적을 위하여 2차 정확도의 VOF(volume-of-fluid) 기법을 사용하였다. 먼저 선형파를 일정 수심상에서 조파시켜 해석해와 비교한 후 수중구조물이 설치된 지형에 적용하여 파랑의 변형을 수치모의하여 수리모형 실험결과와 비교 및 분석하였다.

• 대한조선학회논문집
• /
• 제49권2호
• /
• pp.196-202
• /
• 2012

In this study, numerical analysis of viscous flows is carried out based on the unstructured grid. There exist some difficulties in expressing and computing numerical derivatives on the unstructured grid due to lack of the structured characteristics. The general computer algorithms are developed to perform numerical derivatives easily and extended to be applicable to various geometries composed of hybrid meshes. And the optimal method of strongly implicit procedure is newly contrived to accelerate the rate of convergence in solving the pressure Poisson equation. To verify numerical schemes, the driven cavity problems of 2 and 3 dimension are simulated. The numerical results are compared with others and our numerical schemes are shown to be valid.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제5권1호
• /
• pp.55-78
• /
• 2001

Given the anisotropic Poisson equation $-{\nabla}{\cdot}{\mathcal{K}}{\nabla}p=f$, one can convert it into a system of two first order PDEs: the Darcy law for the flux $u=-{\mathcal{K}{\nabla}p$ and conservation of mass ${\nabla}{\cdot}u=f$. A very natural mixed finite volume method for this system is to seek the pressure in the nonconforming P1 space and the Darcy velocity in the lowest order Raviart-Thomas space. The equations for these variables are obtained by integrating the two first order systems over the triangular volumes. In this paper we show that such a method is really a standard finite element method with local recovery of the flux in disguise. As a consequence, we compare two approaches in analyzing finite volume methods (FVM) and shed light on the proper way of analyzing non co-volume type of FVM. Numerical results for Dirichlet and Neumann problems are included.

• Environmental Engineering Research
• /
• 제11권3호
• /
• pp.164-171
• /
• 2006

Numerical modeling of the flow velocity fields for the near corona wire electrohydrodynamic (EHD) flow was conducted. The steady, two-dimensional momentum equations have been computed for a wire-plate type electrostatic precipitator (ESP). The equations were solved in the conservative finite-difference form on a fine uniform rectilinear grid of sufficient resolution to accurately capture the momentum boundary layers. The numerical procedure for the differential equations was used by SIMPLEST algorithm. The Phoenics (Version 3.5.1) CFD code, coupled with Poisson's electric field, ion transport equations and the momentum equation with electric body force were used for the numerical simulation and the Chen-Kim ${\kappa}-{\varepsilon}$ turbulent model numerical results that an EHD secondary flow was clearly visible in the downstream regions of the corona wire despite the low Reynolds number for the electrode ($Re_{cw}=12.4$). Secondary flow vortices caused by the EHD increases with increasing discharge current or EHD number, hence pressure drop of ESP increases.