• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.2558-2561
• /
• 2008

In the present study, a parallel Taylor-Galerkin/level set based two-phase flow code was developed using finite element discretization and domain decomposition method based on MPI (Message Passing Interface). The proposed method can be utilized for the analysis of a large scale free surface problem in a complex geometry due to the feature of FEM and domain decomposition method. Four-step fractional step method was used for the solution of the incompressible Navier-Stokes equations and Taylor-Galerkin method was adopted for the discretization of hyperbolic type redistancing and advection equations. A Parallel ILU(0) type preconditioner was chosen to accelerate the convergence of a conjugate gradient type iterative solvers. From the present parallel numerical experiments, it has been shown that the proposed method is applicable to the simulation of large scale free surface flows.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.26 no.6_spc
• /
• pp.721-728
• /
• 2016

The combined rotational and translational dynamic vibration absorbers (DVA) with no dampers for the beam vibration control can effectively isolate the vibration within the external excitation force region. This paper investigates the damping efficacy for the combined rotational and translational dynamic vibration absorbers to impose some robustness to the DVA system for the excitation force frequency variation. The beam is assumed to be subjected to a concentrated harmonic excitation force. The solution to the problem is found based on Galerkin method.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.11 no.6
• /
• pp.208-215
• /
• 2001

The paper presents both theoretical and experimental study fur dynamic instabilities of a vortical cantilevered pipe with two attached lumped masses conveying fluid. The two attached lumped masses can be considered as valves or some mechanical paras in real pipe systems. Eigenvalue behaviors depending on the flow velocity are investigated for the change of Positions and magnitudes of an attached lumped mass and a tip mass. In order to verify appropriate of numerical solutions, experiments were accomplished. Theoretical predictions have a good agreement with experimental ones.

• Journal of the Korea Concrete Institute
• /
• v.24 no.5
• /
• pp.567-575
• /
• 2012

A research was conducted to develop a 2-D nonlinear Galerkin finite element analysis of reinforced concrete structures subjected to high temperature with experiments. Algorithms for calculating the closed-form element stiffness for a triangular element with a fully populated material conductance are developed. The validity of the numerical model used in the program is established by comparing the prediction from the computer program with results from full-scale fire resistance tests. Details of fire resistance experiments carried out on reinforced concrete slabs, together with results, are presented. The results obtained from experimental test indicated in that the proposed numerical model and the implemented codes are accurate and reliable. The changes in thermal parameters are discussed from the point of view of changes of structure and chemical composition due to the high temperature exposure. The proposed numerical model takes into account time-varying thermal loads, convection and radiation affected heat fluctuation, and temperature-dependent material properties. Although, this study considered standard fire scenario for reinforced concrete slabs, other time versus temperature relationship can be easily incorporated.

• Journal of computational fluids engineering
• /
• v.19 no.3
• /
• pp.52-61
• /
• 2014

The present paper deals with the efficient computation of higher-order CFD methods for compressible flow using graphics processing units (GPU). The higher-order CFD methods, such as discontinuous Galerkin (DG) methods and correction procedure via reconstruction (CPR) methods, can realize arbitrary higher-order accuracy with compact stencil on unstructured mesh. However, they require much more computational costs compared to the widely used finite volume methods (FVM). Graphics processing unit, consisting of hundreds or thousands small cores, is apt to massive parallel computations of compressible flow based on the higher-order CFD methods and can reduce computational time greatly. Higher-order multi-dimensional limiting process (MLP) is applied for the robust control of numerical oscillations around shock discontinuity and implemented efficiently on GPU. The program is written and optimized in CUDA library offered from NVIDIA. The whole algorithms are implemented to guarantee accurate and efficient computations for parallel programming on shared-memory model of GPU. The extensive numerical experiments validates that the GPU successfully accelerates computing compressible flow using higher-order method.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.14 no.5
• /
• pp.458-465
• /
• 2003

In this paper, we analyze the transient electromagnetic response from three-dimensional(3-D) dielectric bodies using a time-domain electric field integral equation formulation. The solution method in this paper is based on the Galerkin＇s method that involves separate spatial and temporal testing procedures. Triangular patch basis functions are used for spatial expansion and testing functions for arbitrarily shaped 3-D dielectric structures. The time-domain unknown coefficients of the equivalent electric and magnetic currents are approximated as an orthonormal basis function set that is derived from the Laguerre functions. These basis functions are also used as the temporal testing. Numerical results involving equivalent currents and far fields computed by the proposed method are presented.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.14 no.10
• /
• pp.1096-1103
• /
• 2003

In this paper, we analyze the transient electromagnetic response from three-dimensional(3-D) dielectric bodies using a time-domain PMCHW(Poggio, Miller, Chang, Harrington, Wu) formulation. The solution method in this paper is based on the Galerkin's method that involves separate spatial and temporal testing procedures. Triangular patch basis functions are used for spatial expansion and testing functions for arbitrarily shaped 3-D dielectric structures. The time-domain unknown coefficients of the equivalent currents are approximated by a set of orthonormal basis functions that are derived from the Laguerre polynomials. These basis functions are also used as the temporal testing. Numerical results involving equivalent currents and far fields computed by the proposed method are presented.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.19 no.5
• /
• pp.857-871
• /
• 1994

In this paper, six kinds of spiral antenna, a combination of two types of spiral arm-width and three types of spiral curvature are analyzed by using moment method. Dividing spiral arms into N sections, the current distribution is calculated by Galerkin`s method. The radiation pattern and the antenna gain are derived from antenna currents. All os the six spiral antenna have amni-dirctional and wide-band characteristics, although the antenna gain changes within +_ 5dB bound for operating range(600MHz-2GHz). The variation of antenna`s gain is caused by the return loss in connection the Balun to the antenna. Simulation and experimental results on the radiation pattern also show spiral antennas have omni-directional and wide-band characteristics.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.3
• /
• pp.322-329
• /
• 2001

A study of the shock wave turbulent boundary layer interaction is presented. The focus of the study is the interactions of the shock waves with the turbulent boundary layer on the falt plate. Three examples are investigated. The computations are performed, using mixed explicit-implicit generalized Galerkin finite element method. The linear equations at each time step are solved by a preconditioned GMRES algorithm. Numerical results indicate that the implicit scheme converges to the asymptotic steady state much faster than the explicit counterpart. The computed surface pressures and skin friction coefficients display good agreement with experimental data. The flowfield manifests a complex shock wave system and a pair of counter-rotating vortices.

• Nuclear Engineering and Technology
• /
• v.26 no.2
• /
• pp.212-224
• /
• 1994

A numerical method has been developed for studying the dynamics of a flexible cylinder in a coaxial cylindrical duct, immersed in inviscid flow. The unsteady inviscid fluid-dynamic force acting on the oscillating cylinder has been estimated more rigorously by means of a spectral collocation method without simplification of governing equations. This numerical approach is applicable to the system haying wider annular gap and/or shorter length of cylinder as compared to existing potential theory. The governing equation of the unsteady flow was obtained from Laplace equation. The equation of cylinder motion coupled with the fluid motion was discretized by Galerkin's method, from which the dynamic behaviour of the system has been evaluated. The effect of the length of the cylinder and the annular gap on the critical flour velocity, where the system loses stability by buckling, was investigated. To validate the numerical method, the potential flow theory developed by Hobson based on thin film approximation has been improved. Typical results of the present numerical theory on the dynamics and stability of the system are compared with those of available existing theory and the present approximate results. Good agreement was found between the results. It was also found that a nondimensional critical flow velocity becomes larger as increasing the annular gap and decreasing the length of cylinder.