• Geomechanics and Engineering
• /
• 제2권2호
• /
• pp.141-156
• /
• 2010

In the framework of meshfree methods, a new methodology is developed based on radial point interpolation method (RPIM). This methodology is applied to a one-dimensional contaminant transport modelling in the saturated porous media. The one-dimensional form of advection-dispersion equation involving reactive contaminant is considered in the analysis. The Galerkin weak form of the governing equation is formulated using 1D meshfree shape functions constructed using thin plate spline radial basis functions. MATLAB code is developed to obtain the numerical solution. Numerical examples representing various phenomena, which occur during migration of contaminants, are presented to illustrate the applicability of the proposed method and the results are compared with those obtained from the analytical and finite element solutions. The proposed RPIM has generated results with no oscillations and they are insensitive to Peclet constraints. In order to test the practical applicability and performance of the RPIM, three case studies of contaminant transport through the landfill liners are presented. A good agreement is obtained between the results of the RPIM and the field investigation data.

• 대한기계학회논문집B
• /
• 제26권12호
• /
• pp.1738-1746
• /
• 2002

In the FSI (Fluid-Structure Interaction) problems, two different governing equations are to be solved together. One is fur the fluid and the other for the structure. Furthermore, a kinematic constraint should be imposed along the boundary between the fluid and the structure. We use the combined formulation, which incorporates both the fluid and structure equations of motion into a single coupled variational equation so that it is not necessary to calculate the fluid force on the surface of structure explicitly when solving the equations of motion of the structure. A two-dimensional channel flow divided by a Bernoulli-Euler beam is considered and the dynamic response of the beam under the influence of channel flow is studied. The Navier-Stokes equations are solved using a P2P1 Galerkin finite element method with ALE (Arbitrary Lagrangian-Eulerian) algorithm. The internal structural damping effect is not considered in this study and numerical results are compared with a previous work fer steady case. In addition to the Reynolds number, two non-dimensional parameters, which govern this fluid-structure system, are proposed. It is found that the larger the dynamic viscosity and density of the fluid are, the larger the damping of the beam is. Also, the added mass is found to be linearly proportional to the density of the fluid.

• Korean Chemical Engineering Research
• /
• 제47권1호
• /
• pp.72-78
• /
• 2009

전기전도성 이방성 복합재료의 방전가공에 대하여 비정상상태 수식모델을 세우고 갤러킨의 유한요소법으로 해를 구하였다. 피삭재의 온도 분포와 분화구의 모양 및 공작물 제거 속도를 공정 매개변수에 관하여 구득하였다. 계산의 정확도와 효율을 위하여 앞선 연구에서 최적치로 선정된 $12{\times}12$ 요소의 비규칙 체눈을 사용하였다. 알루미나/티타늄 카바이드 복합재료의 물성을 재료의 물성으로 선정하였고 51.4 V의 전압과 7 A의 전류를 갖는 전력을 적용하였으며 제거 효율을 10%로 전열 이방성 계수를 2와 3으로 가정하였다. 불꽃이 일어나면서 피삭재는 즉시 녹기 시작하였고 열적 손상 영역이 형성되었다. 또한 시간이 흘러감에 따라서 분화구의 경계가 이동하는 것이 확인되었다. 반경 방향과 축 방향의 열전도도가 독립적으로 커지면 온도분포와 분화구의 모양이 각각 반경 방향과 축 방향으로 이동하였다. 공작물 제거 속도는 축 방향의 열전도도보다 반경 방향의 열전도도가 증가할 때 더욱 커지는 것으로 나타났다.

• 대한기계학회논문집B
• /
• 제25권7호
• /
• pp.989-996
• /
• 2001

High-speed and low-speed flows are simulated numerically by flowfield-dependent mixed explicit-implicit (FDMEI) method. This algorithm depends on implicitness parameters of convection, diffusion, diffusion gradients, and source terms which are calculated from the changes of local Mach, Reynolds, Peclet, and Damkohler numbers between adjacent nodes. Convection phenomena or shock waves are resolved from Mach number-dependent implicitness parameters whereas diffusion or viscous actions are simulated by Reynolds number or Peclet number-dependent implicitness parameters. Fluctuation components of all variables are properly accommodated spatially and temporally in the FDMEI procedure. To illustrate, some benchmark example problems are presented for comparisons of the FDMEI results with other available data. These results appear to be encouraging and point toward the need for further investigations of the FDMEI theory.

• 한국철도학회:학술대회논문집
• /
• 한국철도학회 2007년도 추계학술대회 논문집
• /
• pp.1605-1609
• /
• 2007

In the case of Linear induction motor(LIM), numerical analysis method like Finite Element Method(FEM) has been mainly used to analyze the travelling magnetic field problem which includes the velocity-induced electromotive force. If the problem including the velocity-induced electromotive force is analyzed by FEM using Galerkin method, the solution can be oscillated according to Peclet Number, which is determined by conductivity, permeability, moving velocity and size of mesh. Consequently, the accuracy of the solution can be low and the Loops of flux can be occurred at the secondary back-iron. These loops of flux occurred at the secondary back-iron does not exist physically, but it can be occurred in the analysis. In this case, the loops of flux can be generally removed by using Up-Wind method which is impossible to apply a conventional S/W tool(Maxwell 2D). Therefore, in this paper, authors examined the Loops of flux occurred at the secondary back-iron of LIM according to variations of Peclet Number, and analyzed whether these loops of flux affect on the thrust force characteristics of LIM or not.

• Structural Engineering and Mechanics
• /
• 제47권5호
• /
• pp.599-622
• /
• 2013

Seismic response of two dimensional liquid tanks is numerically simulated using fully nonlinear velocity potential theory. Galerkin-weighted-residual based finite element method is used for solving the governing Laplace equation with fully nonlinear free surface boundary conditions and also for velocity recovery. Based on mixed Eulerian-Lagrangian (MEL) method, fourth order explicit Runge-Kutta scheme is used for time integration of free surface boundary conditions. A cubic-spline fitted regridding technique is used at every time step to eliminate possible numerical instabilities on account of Lagrangian node induced mesh distortion. An artificial surface damping term is used which mimics the viscosity induced damping and brings in numerical stability. Four earthquake motions have been suitably selected to study the effect of frequency content on the dynamic response of tank-liquid system. The nonlinear seismic response vis-a-vis linear response of rectangular liquid tank has been studied. The impulsive and convective components of hydrodynamic forces, e.g., base shear, overturning base moment and pressure distribution on tank-wall are quantified. It is observed that the convective response of tank-liquid system is very much sensitive to the frequency content of the ground motion. Such sensitivity is more pronounced in shallow tanks.

• 대한기계학회논문집
• /
• 제17권9호
• /
• pp.2304-2314
• /
• 1993

Natural convective flow and heat transfer in a two-dimensional square enclosure fitted with a horizontal partition are investigated numerically. The enclosure was composed of the lower hot and the upper cold horizontal walls and the adiabatic vertical walls, and a partition was situated perpendicularly at the one vertical insulated wall. The governing equations are solved by using the finite element method with Galerkin method. The computations were carried out with the variations of length, position and thermal conductivity of the partition, and Rayleigh number based on the temperature difference between two horizontal walls and the enclosure height with water(Pr=4.95). As the results, an oscillatory motion of natural convection is resulted in a sudden rise of overall heat transfer, but the increase of length of partition is significantly restrained the increase of Nusselt number. The maximum heat transfer was shown just before the transition of the direction of oscillating flow. An oscillatory motion of flow was perfectly shown the stability with the decrease of the length of partition and Rayleigh number. Also, the heat transfer was raised with the increase of the thermal conductivity in proportion to the increase of the length of partition. The stability and oscillation of flow are affected by the position of partition.

• 설비공학논문집
• /
• 제5권4호
• /
• pp.285-294
• /
• 1993

An oscillatory motion of natural convection in a two-dimensional square enclosure fitted with a horizontal partition is investigated numerically. The enclosure was composed of the lower hot and the upper cold horizontal walls and the adiabatic vertical walls, and a partition was positioned perpendicularly at the mid-height of one vertical insulated wall. The governing equations are solved by using the finite element method with Galerkin method. The computations were carried out with the variations of the partition length and Rayleigh number based on the temperature difference between two horizontal walls and the enclosure height with water(Pr=4.95). As the results, an oscillatory motion of natural convection has perfectly shown the periodicity with the decrease of Rayleigh number, and the stability was reduced to a chaotic state with the increase of Rayleigh number. The period of oscillation gets shorten with the decrease of the partition length and the increase of Rayleigh number. The frequency of oscillation obtained by the variations of stream function is more similar to the experimental results than that of the average Nusselt number. The stability of oscillation grows worse with the increase of Rayleigh number. The transition Rayleigh number for the chaos is gradually decreased with the increase of the partition length.

• 설비공학논문집
• /
• 제9권3호
• /
• pp.312-322
• /
• 1997

The natural convective flow in a two-dimensional square enclosure with horizontal partitions is investigated numerically. The enclosure was composed of the lower hot and the upper cold horizontal walls and the adiabatic vertical walls, and two identical partitions were positioned perpendicularly at the mid-height of the right and left walls, respectively. The governing equations are solved by using the finite element method with Galerkin method. Calculations are made for different partition lengths, partition conductivites, and Rayleigh numbers based on the temperature difference between two horizontal walls and the enclosure height with water(Pr=4.95). An oscillatory motion of the natural convective flow is affected significantly by the variation of the gap width and Rayleigh number. When the gap width is comparatively short, the heat transfer rate is raised with the increase of the thermal conductivity of partitions. However, for sufficiently large gap widths at higher Rayleigh numbers, the average Nusselt numbers of the conductive partitions are smaller than those of the adiabatic partitions.

• Steel and Composite Structures
• /
• 제26권1호
• /
• pp.17-29
• /
• 2018

The thermal effects on the buckling, postbuckling and nonlinear vibration behaviors of composite laminated trapezoidal plates are studied. Aiming at the complex plate structure and to simulate the temperature distribution of the plate, a finite element method (FEM) is applied in this paper. In the temperature model, based on the thermal diffusion equation, the Galerkin's method is employed to establish the temperature equation of the composite laminated trapezoidal plate. The geometrical nonlinearity of the plate is considered by using the von Karman large deformation theory, and combining the thermal model and aeroelastic model, Hamilton's principle is employed to establish the thermoelastic equation of motion of the composite laminated trapezoidal plate. The thermal buckling and postbuckling of the composite laminated rectangular plate are analyzed to verify the validity and correctness of the present methodology by comparing with the results reported in the literature. Moreover, the effects of the temperature with the ply-angle on the thermal buckling and postbuckling of the composite laminated trapezoidal plates are studied, the thermal effects on the nonlinear vibration behaviors of the composite laminated trapezoidal plates are discussed, and the frequency-response curves are also presented for the different temperatures and ply angles.