• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.30 no.7 s.250
• /
• pp.833-840
• /
• 2006

The paper presents a numerical solution to the problem of a hot rigid indenter sliding over a thermoelastic Winkler foundation with a thermal contact resistance at constant speed. It is shown analytically that no steady-state solution can exist for sufficiently high temperature or sufficiently small normal load or speed, regardless of the thermal contact resistance. However the steady state solution may exist in the same situation if the thermal contact resistance is considered. This means that the effect of the large values of temperature difference and small value of force or velocity which occur at no steady state can be lessened due to the thermal contact resistance. When there is no steady state, the predicted transient behavior involves regions of transient stationary contact interspersed with regions of separation regardless of the thermal contact resistance. Initially, the system typically exhibits a small number of relatively large contact and separation regions, but after the initial transient, the trailing edge of the contact area is only established and the leading edge loses contact, reducing the total extent of contact considerably. As time progresses, larger and larger numbers of small contact areas are established, unlit eventually the accuracy of the algorithm is limited by the discretization used.

• Journal of the Korean earth science society
• /
• v.39 no.4
• /
• pp.317-326
• /
• 2018

Effect of the nonuniform grid on the two-dimensional transport equation was investigated in terms of theoretical analysis and finite difference method (FDM). The nonuniform grid having a typical structure of the numerical weather forecast model was incorporated in the vertical direction, while the uniform grid was used in the zonal direction. The staggered and non-staggered grid were placed in the vertical and zonal direction, respectively. Time stepping was performed with the third-order Runge Kutta scheme. An error analysis of the spatial discretization on the nonuniform grid was carried out, which indicated that the combined effect of the nonuniform grid and advection velocity produced either numerical diffusion or numerical adverse-diffusion. An analytic function is used for the quantitative evaluation of the errors associated with the discretized transport equation. Numerical experiments with the non-uniformity of vertical grid were found to support the analysis.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.11 no.4
• /
• pp.91-102
• /
• 1991

Traditionally, a dynamic programming(DP) technique has been used to the multireservoir control system. The algorithm has inherent problem to increase computational requirements exponentially due to discretization of variables and expanding the dimension of the system. To solve this problem, this paper describes transforming the optimal control system into a quadratic programming(QP), using Laguerre polynomials(LP) and its properties. The objective function of the proposed QP is independent of time variable. The solution of the QP is obtained by nonlinear programming(NLP) using augmented Lagrangian multiplier method. The numerical experiment shows that the water level of reservoirs is higher than Lee's and the evaluated benefit value is about the same as other researcher's.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2005.05a
• /
• pp.330-333
• /
• 2005

The paper presents a numerical solution to the problem of a hot rigid indenter siding over a thermoelastic Winkler foundation with a thermal contact resistance at constant speed. It is shown analytically that no steady-state solution can exist for sufficiently high temperature or sufficiently small normal load or speed regardless of the thermal contact resistance. However, the steady state solution may exist in the same situation if the thermal contact resistance is considered. This means that the effect of the large values of temperature difference and small value of force or velocity which occur at no steady state can be lessened due to the thermal contact resistance. When there is no steady-state the predicted transient behavior involves regions of transient stationary contact interspersed with regions of separation regardless of the thermal contact resistance. Initially, the system typically exhibits a small number of relatively large contact and separation regions, but after the initial transient the trailing edge of the contact area is only established and the leading edge loses contact, reducing the total extent of contact considerably. As time progresses, larger and larger number of small contact areas are established, until eventually the accuracy of the algorithm is limited by the discretization used.

• Journal of computational fluids engineering
• /
• v.17 no.2
• /
• pp.35-41
• /
• 2012

This paper presents incompressible Navier-Stokes solution algorithm for 2D Free-surface flow problems on the Cartesian mesh, which was implemented to run on Graphics Processing Units(GPU). The INS solver utilizes the variable arrangement on the Cartesian mesh, Finite Volume discretization along Constrained Interpolation Profile-Conservative Semi-Lagrangian(CIP-CSL). Solution procedure of incompressible Navier-Stokes equations for free-surface flow takes considerable amount of computation time and memory space even in modern multi-core computing architecture based on Central Processing Units(CPUs). By the recent development of computer architecture technology, Graphics Processing Unit(GPU)'s scientific computing performance outperforms that of CPU's. This paper focus on the utilization of GPU's high performance computing capability, and presents an efficient solution algorithm for free surface flow simulation. The performance of the GPU implementations with double precision accuracy is compared to that of the CPU code using an representative free-surface flow problem, namely. dam-break problem.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.2
• /
• pp.115-122
• /
• 2008

The effects of characteristic condition number on the convergence of preconditioned Euler equations were investigated. The two-dimensional preconditioned Euler equations adopting Choi and Merkle's preconditioning and the temperature preconditioning are considered. Preconditioned Roe's FDS scheme was adopted for spatial discretization and preconditioned LU-SGS scheme was used for time integration. It is shown that the convergence characteristics of the Euler equations are strongly affected by the characteristic condition number, and there is an optimal characteristic condition number for a problem. The optimal characteristic condition numbers for the Choi and Merkle's preconditioning and temperature preconditioning are different.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.2
• /
• pp.123-130
• /
• 2008

The effects of characteristic condition number on the convergence of preconditioned Navier-Stokes equations were investigated. The two-dimensional preconditioned Navier-Stokes adopting Choi and Merkle's preconditioning and the temperature preconditioning are considered. Preconditioned Roe's FDS scheme was adopted for spatial discretization and preconditioned LU-SGS scheme was used for time integration. It is shown that the convergence characteristics of the Navier-Stokes equations are strongly affected by the characteristic condition number. Also it is shown that the optimal characteristic condition numbers for viscous flows are larger than that in inviscid flows.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.2
• /
• pp.450-460
• /
• 1996

To investigate the effect of stress wave propagation for crack tip, impact responses of two-dimensional plates with oblique cracks are investigated by a numerical method. In the numerical analysis, the finite element method is used in space domain discretization and the Newmark constant acceleration algorithm is used in time integration. According to the numerical results from the impact response analysis. it is found that the stress fields are bisected at the crack surface and the parts of stress intensity are moved along the crack face. The crack tip stress fields are yaried rapidly. The magnitude of crack tip stress fields are converted to dynamic stress intensity factor. Dynamic sress intensity factor appears when the stress wave has reached at the crack tip and the aspect of change of dynamic stress intensity factor is shown to be the same as the part of the flow of stress intensity.

• Interaction and multiscale mechanics
• /
• v.1 no.1
• /
• pp.1-31
• /
• 2008

The present paper presents multiscale modelling via coupling of the discrete and finite element methods. Theoretical formulation of the discrete element method using spherical or cylindrical particles has been briefly reviewed. Basic equations of the finite element method using the explicit time integration have been given. The micr-macro transition for the discrete element method has been discussed. Theoretical formulations for macroscopic stress and strain tensors have been given. Determination of macroscopic constitutive properties using dimensionless micro-macro relationships has been proposed. The formulation of the multiscale DEM/FEM model employing the DEM and FEM in different subdomains of the same body has been presented. The coupling allows the use of partially overlapping DEM and FEM subdomains. The overlap zone in the two coupling algorithms is introduced in order to provide a smooth transition from one discretization method to the other. Coupling between the DEM and FEM subdomains is provided by additional kinematic constraints imposed by means of either the Lagrange multipliers or penalty function method. The coupled DEM/FEM formulation has been implemented in the authors' own numerical program. Good performance of the numerical algorithms has been demonstrated in a number of examples.

• Structural Engineering and Mechanics
• /
• v.68 no.1
• /
• pp.1-15
• /
• 2018

In ocean industry, free surface type ART (Anti Roll tank) system has been widely used to suppress the roll motion of floating structures. In those, various obstacles have been devised to obtain the sufficient damping and to enhance the controllability of freely rushing water inside the tank. Most of previous researches have paid on the development of simple mathematical formula for coupled ship-ARTs analysis although other numerical and experimental approaches exist. Little attention has been focused on the use of 3D panel method for preliminary design of free surface type ART despite its advantages in computational time and general capacity for hydrodynamic damping estimation. This study aims at developing a potential theory based hydrodynamic code for the analysis of floating structure with baffled ARTs. The sloshing in baffled tanks is modeled through the linear potential theory with FE discretization and it coupled with hydrodynamic equations of floating structures discretized by BEM and FEM, resulting in direct coupled FE-BE formulation. The general capacity of proposed formulation is emphasized through the coupled hydrodynamic analysis of floating structure and sloshing inside baffled ARTs. In addition, the numerical methods for natural sloshing frequency tuning and estimation of hydrodynamic damping ratio of liquid sloshing in baffled tanks undergoing wave exiting loads are developed through the proposed formulation. In numerical examples, effects of natural frequency tuning and baffle ratios on the maximum and significant roll motions are investigated.