• Computers and Concrete
• /
• v.25 no.5
• /
• pp.411-425
• /
• 2020

In this paper, the shell material has been taken as functionally graded material and their material quantity is located by the exponential volume fraction law. Moreover, the impact of ring supports around the shell circumference has been examined for their various positions along the shell axial length. These rings support restraints the radial displacement in the transverse direction. While the axial modal deformation functions have been estimated by characteristic beam functions and nature of materials used for construction of cylindrical shells. The fundamental natural frequency of cylindrical shell of parameter versus ratios of length- and height-to-radius for a wide range has been reported and investigated through the study. In addition, by increasing height-to-radius ratio resulting frequencies also increase and frequencies decrease on ratio of length-to-radius. Though the trends of frequency values of both ratios are converse to each other with three different boundary conditions. Also it is examined the position of ring supports with length-to radius ratio, height-to-radius ratio and varying the exponent of volume fraction. MATLAB software package has been utilized for extracting shell frequency spectra. The obtained results are confirmed by comparing with available literature.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.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.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.9 no.3
• /
• pp.65-73
• /
• 1989

Continuum formulations for the expressions of dynamic energy release rates and computational methods for dynamic stress intensity factors are developed for the analysis of dynamic fracture problems subjected to stress wave loading. Explicit volume integral expressions for instantaneous dynamic energy release rates are derived by modeling virtual crack extensions with the dynamic Eulerian-Lagrangian kinematic description. In the finite element applications a finite region around a crack-tip is modeled by using quarter-point singular isoparametric elements, and the volume integrals are evaluated for each crack-tip element during virtual crack extensions while the singularity is maintained. It is shown that the use of the present method is more reliable and accurate for the dynamic fracture analysis than that of other path-independent integral methods when the effects of stress waves are significant.

• Journal of the Korean Society of Combustion
• /
• v.8 no.3
• /
• pp.15-23
• /
• 2003

Partial quenching structure of turbulent diffusion flames in a turbulent mixing layer is investigated by the method of flame hole dynamics in order to develop a prediction model for turbulent flame lift off. The essence of flame hole dynamics is derivation of the random walk mapping, from the flame-edge theory, which governs expansion or contraction of flame holes initially created by local quenching events. The numerical simulation for flame hole dynamics is carried out in two stages. First, a direct numerical simulation is performed for constant-density fuel-air channel mixing layer to obtain the turbulent flow and mixing fields, from which a time series of two dimensional scalar dissipation rate array is extracted at a fixed virtual flame surface horizontally extending from the end of split plate to the downstream. Then, the Lagrangian simulation of the flame hole random walk mapping projected to the scalar dissipation rate array yields temporally evolving turbulent extinction process and its statistics on partial quenching characteristics. The statistical results exhibit that the chance of partial quenching is strongly influenced by the crossover scalar dissipation rate while almost unaffected by the iteration number of the mapping that can be regarded as a flame-edge speed.

• Journal of the Korean Society for Precision Engineering
• /
• v.28 no.4
• /
• pp.391-397
• /
• 2011

Control of interior permanent magnet (IPMSM) is difficult because its nonlinearity and parameter uncertainty. In this paper, a fuzzy c-regression models clustering algorithm which is based on T-S fuzzy is used to model IPMSM with a series linear model and weight them by memberships. Lagrangian of constrained function is built for calculating clustering centers where training output data are considered. Based on these clustering centers, least square method is applied for T-S fuzzy linear model parameters. As a result, IPMSM can be modeled as T-S fuzzy model for T-S fuzzy control of them.

• Journal of Korean Society of Occupational and Environmental Hygiene
• /
• v.28 no.2
• /
• pp.190-199
• /
• 2018

Objectives: The purpose of this case study is to assess the current airflow and find the ideal ventilation conditions in tank reactors for minimizing the possibility of exposure respiratory dusts(size of $2.5{\mu}m$, $10{\mu}m$) when workers exchange catalysts in the tank reactors. Methods: A Numerical study was performed to determine ideal ventilation conditions, We considered two sizes of airborne respiratory particles($2.5{\mu}m$, $10{\mu}m$) at 12points from the bottom of tank reactor. We changed input & output ventilation conditions and analyzed the particle motion in the tank reactor. The star-ccm+, computational fluid dynamics tool was used to predict air & particle flow patterns in the tank reactor and a numerical simulation was achieved by applying the realized ${\kappa}-{\varepsilon}$ turbulence model and the Lagrangian particle tracking method. Results: From the results, the increase of recirculation air had a significant impact on removing dusts because they are removed by HEPA filter. To the contrary, Increasing the clean air quantity or changing the input position of clean air is not good for workers because it causes the exit of respiratory dusts through workers' entrance or cause it to remail suspended in the air in the workplace tank.

• Proceedings of the Korea Society for Energy Engineering kosee Conference
• /
• 1995.11a
• /
• pp.50-55
• /
• 1995

향상되어진 93-PCGC-2는 기존의 PCGC-2와 같이 미분탄 연소를 포함하는 다양한 반응성흐름과 비반응성 흐름을 설명하기 위해 2차원 정상상태 모델로 제시되어 졌다. 93-PCGC-2는 실린더형의 축 대칭계에 응용되어질 수 있고, 난류(Turbulence)는 유체역학식과 연소기구 양쪽을 위해 고려되어졌으며, 불연속 세로좌표 방법(Discrete Ordinates Method)을 이용하여 기체, 벽 및 입자들로부터의 복사열(Radiation)을 모사하였다. 입자상은 입자 무리들의 평균 경로들을 따라 해석하는 Lagrangian계의 해석법으로 모델화되어졌다. 석탄의 팽윤(Swelling)과 촤의 반응성에 관한 부모델과 더불어 새롭게 일반화된 석탄 탈휘발화 부모델 (FG-DVC)도 첨가되어졌다. 비균일 반응기구는 확산과 화학반응 둘 모두를 고려하였다. 주요 기상반응은 국부 순간 평형을 가정하여 모델화하였다. 그래서 반응속도는 혼합의 난류속도에 의해 제한되어진다. Thermal NOx과 Fuel NOx의 유한속도 화학론(Finite Rate Chemstry)에 대한 부모델은 화학반응속도론와 난류성의 통계치를 통합하여 만들어져 있다. 기상은 반복적인 line-by-line기교에 의해 풀려지는 elliptic partial differential equation으로 묘사되어진다. 수치적인 안정을 고려하기 위해 under-relaxation이 이용되어졌다. 이렇게 코드화된 93-PCGC-2는 연소를 위해 모사되어졌다. 또한 더 나아가 이 수치모델의 활용범위는 미분탄의 가스화에도 활용되어질 것으로 기대되어진다.

• Journal of Theoretical and Applied Mechanics
• /
• v.2 no.1
• /
• pp.31-50
• /
• 1996

In this paper, a finite element analysis based on the local approach concept to fracture in the continuum damage mechanics is performed to analyze ductile fracture in two dimensional quasi-static state. First an isotropic damage model based on the generalized concept of effective stress is proposed for structural materials in the context of large deformation. In this model, the stiffness degradation is taken as a measure of damage and so, the fracture phenomenon can be explained as the critical deterioration of stiffness at a material point. The modified Riks' continuation technique is used to solve incremental iterative equations. Crack propagation is achieved by removing critically damaged elements. The mesh size sensitivity analysis and the simulation of the well known shearing mode failure in plane strain state are carried out to verify the present formulation. As numerical examples, an edge cracked plate and the specimen with a circular hole under plane stress are taken. Load-displacement curves and successively fractured shapes are shown. From the results, it can be concluded that the proposed model based on the local approach concept in the continuum damage mechanics may be stated as a reasonable tool to explain ductile fracture initiation and crack propagation.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.86-97
• /
• 1993

A numerical method is developed for the nonlinear motion of two-dimensional wedges and axisymmetric-forced-heaving motion using Semi-Largrangian scheme under assumption of potential flows. In two-dimensional-problem Cauchy's integral theorem is applied to calculate the complex potential and its time derivative along boundary. In three-dimensional-problem Rankine ring sources are used in a Green's theorem boundary integral formulation to salve the field equation. The solution is stepped forward numerically in time by integrating the exact kinematic and dynamic free-surface boundary condition. Numerical computations are made for the entry of a wedge with a constant velocity and for the forced harmonic heaving motion from rest. The problem of the entry of wedge compared with the calculated results of Champan[4] and Kim[11]. By Fourier transform of forces in time domain, added mass coefficient, damping coefficient, second harmonic forces are obtained and compared with Yamashita's experiment[5].

• Bulletin of the Society of Naval Architects of Korea
• /
• v.23 no.2
• /
• pp.1-13
• /
• 1986

The nonlinear hydrodynamic forces acting on a two-dimensional circular cylinder, oscillating with large amplitude in the free surface, are calculated by using the Semi-Lagrangian Time-Step-ping Method used by O.M. Faltinsen. In present calculation the position and the potential value of free surface are calculated using the exact kinematic and dynamic free surface boundary condition. At each time step an integral equation is solved to obtain the value of potential and normal velocity along the boundaries, consisting of both the body surface and the free surface. Some effort was devoted to the elimination of instability arising in the range of high frequency. Numerical simulations were performed up to the 3rd or 4th period which seems to be enough for the transient effect to die out. Each harmonic component and time-mean force are obtained by the Fourier transform of forces in time domain. The results are compared with others' experimental and theoretical results. Particularly, the calculation shows the tendency that the acceleration-phase 1st-harmonic component(added mass) increases as the motion amplitude increases and a reverse tendency in the velocity-phase 1st-harmonic component(damping coefficient). The Yamashita's experimental result also shows the same tendency. In general, the present result show relatively good agreement with the Yamashita's experimental result except for the time-mean force.