• 대한기계학회논문집B
• /
• 제24권12호
• /
• pp.1555-1569
• /
• 2000

Numerical simulation of fluid flow with moving free surface has been carried out. For the free surface flow, a VOF(Volume of Fluid)-based algorithm utilizing a fixed grid system has been investigated. In order to reduce numerical smearing at the free surface represented on a fixed grid system, a new free surface tracking algorithm based on the donor-acceptor scheme has been presented. Novel features of the proposed algorithm are characterized as two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby-cell to determine the fluid volume flux at each cell boundary. The proposed algorithm can be easily implemented in any irregular non-uniform grid systems that are usual in finite element method (FEM). Moreover, the proposed algorithm can be extended and applied to the 3-D free surface flow problem without additional efforts. For computation of unsteady incompressible flow, a finite element approximation based on the explicit fractional step method has been adopted. In addition, the SUPG(streamline upwind/Petrov-Galerkin) method has been implemented to deal with convection dominated flows. Combination of the proposed free surface tracking scheme and explicit fractional step formulation resulted in an efficient solution algorithm. Validity of the present solution algorithm was demonstrated from its application to the broken dam and the solitary wave propagation problems.

• 대한조선학회논문집
• /
• 제29권4호
• /
• pp.114-131
• /
• 1992

본 논문은 2차원 자유표면파문제에서 시간영역해법을 이용하여 2차원 운동문제에 적용할 수 있는 수치해석을 하였다. 경계조건으로는 엄밀한 물체표면 경계조건과 비선형 자유표면경계 조건을 부과했다. 수치해를 구하는데는 코시이론을 이용하여 제2종 프레드흘름 경계적분방정식을 도출하고 이를 이산화시켜 처리하였다. 수치계산을 위해 전영역을 유한한 영역으로 제한하여야 한다. 제한된 영역에서 방사해의 부과를 위해 전영역을 수치해석영역과 외부영역으로 나누고, 외부영역의 해는 그린 제2정리를 이용하고, 선형자유표면조건을 만족하는 과도그린함수를 사용한다. 위의 그린 제2정리를 이용한 식으로 부터 초기조건, 선형 자유표면조건, 무한원방조건을 이용하여 단순화시킨 다음 포텐셜과 유동함수의 관계식으로 치환하면 비선형해와 정합할 수 있는 정합행렬을 구할 수 있다. 본 논문에서 개발한 정합방법을 이용하여 적용할 문제로서 첫째는 주상체가 상하동요, 수평동요하는 경우 계산이고 두번째로는 수면하에서 타원형실린더가 일정속도로 항진할 때 계산을 수행한 결과를 고차스펙트럴방법과 비교하였다.

• Journal of Korean Society for Atmospheric Environment
• /
• 제9권E호
• /
• pp.347-357
• /
• 1993

The Lagrangian grid-free numerical method, the discrete vortex method, was applied to solve the Navier-Stokes euqations. This method avoids the introduction of numerical viscosity swamping the real physical viscosity at high Reynolds number, unlike Eulerian method, e.g. finite difference and element methods. The boundary integral equation method for the potential flow solution was included to make the discrete vortex method more feasible for complex geometries. The fast adaptive multipole expansion method was incorporated to reduce the computational time from $O(N^2)$ to O(N) for the computations of vortex-vortex interactions. The test problems were air flow around one circular cylinder and two circular cylinders in tandem with various gaps. The numerical results were in excellent gareement with the experimental and other computational results. The applicabilty of the method was discussed with the indoor and the outdoor air pollution problems, especially the contaminant transport in the recirculation regions.

• 접착 및 계면
• /
• 제6권1호
• /
• pp.1-6
• /
• 2005

The high accurate simulation of very thin glue layers based on the finite element method is still connected to many problems which result from the necessity to construct a complicated mesh of essentially different sizes of elements. This can lead to a loss of accuracy, unstable calculations and even loss of convergence. However, the implementation of special transmission elements along the glue ling and special edge-elements in the near-edge region would lead to a dramatic decrease of number of finite elements in the mesh and thus, prevent unsatisfactory phenomena in numerical analysis and extensive computation time. The theoretical basis for such special elements is the knowledge about appropriate transmission conditions and the edge effects near the free boundary of the adhesive layer. Therefore, recently proposed so-called non-classical transmission conditions and the behavior near the free edge are investigated in the context of the single-lap tensile-shear test of adhesive technology.

• Korea-Australia Rheology Journal
• /
• 제12권2호
• /
• pp.107-117
• /
• 2000

Prediction of the product shape in sheet casting process is performed from the numerical simulation. A three-dimensional finite element method is used to investigate the flow behavior and to examine the effects of processing conditions on the sheet produced. Effects of inertia, gravity, surface tension and non-Newtonian viscosity on the thickness profile of the sheet are considered since the edge bead and the flow patterns in the chill roll region have great influence on the quality of the products. In the numerical simulation with free surface flows, the spine method is adopted to update the free surface, and the force-free boundary condition is imposed along the take-up plane to avoid severe singularity problems existing at the take-up plane. From the numerical results of steady isothermal flows of a generalized Newtonian fluid, it is shown that the draw ratio plays a major role in predicting the shape of the final sheet produced and the surface tension has considerable effect on the bead thickness ratio and the bead width fraction, while shear-thinning and/or tension-thickening viscosity affect the degree of neck-in.

• Earthquakes and Structures
• /
• 제14권2호
• /
• pp.117-128
• /
• 2018

In this research work, free vibrations of simply supported functionally graded plate resting on a Winkler-Pasternak elastic foundation are investigated by a new shear deformation theory. The influence of alternative micromechanical models on the macroscopic behavior of a functionally graded plate based on shear-deformation plate theories is examined. Several micromechanical models are tested to obtain the effective material properties of a two-phase particle composite as a function of the volume fraction of particles which continuously varies through the thickness of a functionally graded plate. Present theory exactly satisfies stress boundary conditions on the top and the bottom of the plate. The energy functional of the system is obtained using Hamilton's principle. The closed form solutions are obtained by using Navier technique, and then fundamental frequencies are found by solving the results of eigenvalue problems. Finally, the numerical results are provided to reveal the effect of explicit micromechanical models on natural fundamental frequencies.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2004년도 추계학술대회논문집
• /
• pp.898-901
• /
• 2004

Since soil-structure interactions are one of the most important subjects in the structural/foundation engineering, much study concerning the soil-structure interactions had been carried out. One of typical structures related to the soil-structure interactions is the strip foundation which is basically defined as the beam or strip rested on or supported by the soils. At the present time, lack of studies on dynamic problems related to the strip foundations is still found in the literature. From these viewpoint, this paper aims to theoretically investigate dynamics of the circular strip foundations and also to present the practical engineering data for the design purpose. Differential equations governing the free, out-of-plane vibrations of such strip foundations are derived, in which effects of the rotatory and torsional inertias and also shear deformation are included although the warping of the cross-section is excluded. Governing differential equations subjected to the boundary conditions of corresponding end constraints are numerically solved for obtaining the natural frequencies and mode shapes by using the numerical integration technique and the numerical method of non-linear equation.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제10권5호
• /
• pp.545-565
• /
• 2018

A numerical code is developed based on potential flow theory to investigate nonlinear sloshing in rectangular Liquefied Natural Gas (LNG) tanks under forced excitation. Using this code, internal free-surface elevation and sloshing loads on liquid tanks can be obtained both in time domain and frequency domain. In the mathematical model, acceleration potential is solved in the calculation of pressure on tanks and the artificial damping model is adopted to account for energy dissipation during sloshing. The Boundary Element Method (BEM) is used to solve boundary value problems of both velocity potential and acceleration potential. Numerical calculation results are compared with published results to determine the efficiency and accuracy of the numerical code. Sloshing properties in partially filled rectangular and membrane tank under translational and rotational excitations are investigated. It is found that sloshing under horizontal and rotational excitations share similar properties. The first resonant mode and excitation frequency are the dominant response frequencies. Resonant sloshing will be excited when vertical excitation lies in the instability region. For liquid tank under rotational excitation, sloshing responses including amplitude and phase are sensitive to the location of the center of rotation. Moreover, experimental tests were conducted to analyze viscous effects on sloshing and to validate the feasibility of artificial damping models. The results show that the artificial damping model with modifying wall boundary conditions has better applicability in simulating sloshing under different fill levels and excitations.

• Structural Engineering and Mechanics
• /
• 제72권1호
• /
• pp.113-129
• /
• 2019

A four-variable shear deformation refined plate theory has been proposed for dynamic characteristics of smart plates made of porous magneto-electro-elastic functionally graded (MEE-FG) materials with various boundary conditions by using an analytical method. Magneto-electro-elastic properties of FGM plate are supposed to vary through the thickness direction and are estimated through the modified power-law rule in which the porosities with even and uneven type are approximated. Pores possibly occur inside functionally graded materials (FGMs) due the result of technical problems that lead to creation of micro-voids in these materials. The variation of pores along the thickness direction influences the mechanical properties. The governing differential equations and boundary conditions of embedded porous FGM plate under magneto-electrical field are derived through Hamilton's principle based on a four-variable tangential-exponential refined theory which avoids the use of shear correction factors. An analytical solution procedure is used to achieve the natural frequencies of embedded porous FG plate supposed to magneto-electrical field with various boundary condition. A parametric study is led to carry out the effects of material graduation exponent, coefficient of porosity, magnetic potential, electric voltage, elastic foundation parameters, various boundary conditions and plate side-to-thickness ratio on natural frequencies of the porous MEE-FG plate. It is concluded that these parameters play significant roles on the dynamic behavior of porous MEE-FG plates. Presented numerical results can serve as benchmarks for future analyses of MEE-FG plates with porosity phases.

• 한국전산유체공학회지
• /
• 제20권2호
• /
• pp.103-108
• /
• 2015

Fluid-body interaction analysis of floating body with six degree-of-freedom motion is presented. In this study, three-dimensional incompressible Navier-Stokes equations are employed as a governing equation. The numerical method is based on a finite-volume approach on a cartesian grid together with a fractional-step method. To represent the body motion, the immersed boundary method for direct forcing is employed. In order to simulate the coupled six degree-of-freedom motion, Euler's equations based on rigid body dynamics are utilized. To represent the complex body shape, level-set based algorithm is utilized. In order to describe the free surface motion, the volume of fluid method utilizing the tangent of hyperbola for interface capturing scheme is employed. This study showed three different continuums(air, water and body) are simultaneously simulated by newly developed code. To demonstrate the applicability of the current approach, two different problems(dam-breaking with stationary obstacle and water entry) are simulated and all results are validated.