• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2005년도 학술발표회 논문집
• /
• pp.409-413
• /
• 2005

Even though the relative importance of length scale of flow system allow us to simplify three dimensional flow problem to one or two dimensional representation, many systems still require three dimensional analysis. The objective of this study is to develop an efficient and accurate finite element model for analyzing and predicting three dimensional flow features in natural rivers and to offend to model spreading of pollutants and transport of sediments in the future. Firstly, three dimensional Reynolds averaged Navier-Stokes equations with the hydrostatic pressure assumption in generalized curvilinear coordinates were combined with the kinematic free-surface condition. Secondly. to simulate realistic high Reynolds number flow, the model employed the Streamline Upwind/Petrov-Galerkin(SU/PG) scheme as a weighting function for the finite element method in conjunction with an appropriate turbulence model(Smagorinsky scheme for the horizontal plain and Mellor-Yamada scheme for the vertical direction). Several tests is performed for the purpose of validation and verification of the developed model. A simple rectangular channel, 5-shaped and U-shaped channel are used for tests and comparisons are made with RMA-10 model. Runs for each case is converged stably without a oscillation and calculated water-surface deformation, longitudinal and transversal velocities, and velocity vector fields are in good agreement with the results of RMA-10 model.

• Journal of Theoretical and Applied Mechanics
• /
• 제1권1호
• /
• pp.97-109
• /
• 1995

In this paper, a finite viscoelastic continuum model for rubber and its finite element analysis are presented. This finite viscoelatic model based on continuum mechanics is an extended model of Johnson and Wuigley's 1-D model. In this extended model, continuum based kinematic measures are rigorously defied and by using this kinematic measures, elastic stage law and flow rule are introduced. In kinematics, three configuration are introduced. In kinematics, three configuration are introduced. They are reference, current and virtual visco configurations. In elastic state law, it is assumed that at a certain time, there exists an elastic potential which describes the recoverable elastic energy. From this elastic potential, elastic state law is derived. The proposed flow rule is based on phenomenological observation. The flow rule gives precise relaxation response. In finite element approximation, mixed Lagrangian description is used, where total and similar method of updated Lagrangian descriptions are used together. This approach reduces the numerical job and gives simple nonlinear syatem of equations. To satisfy the incompressible condition, penalty-type modified Mooney-Rivlin energy function is adopted. By this method nearly incompressible condition is obtain the virtual visco configuration. For verification, uniaxial stretch tests are simulated for various stretch rates. The simulated results show good agreement with experiments. As a practical experiments. As a preactical example, pressurized rubber plate is simulated. The result shows finite viscoelastic effects clearly.

• Fibers and Polymers
• /
• 제5권3호
• /
• pp.209-212
• /
• 2004

A nonsteady axi-symmetric ideal flow solution is obtained here. It is based on the rigid perfect-plastic constitutive law with the Tresca yield condition and its associated flow rule. The process is to deform a circular solid disk into a spherical shell of prescribed geometry. It is assumed that there are no rigid zones and friction stresses. The solution obtained provides the distribution of kinematic variables and involves one undetermined function of the time. This function can be in general found by superimposing an optimality criterion.

• 한국수자원학회논문집
• /
• 제33권2호
• /
• pp.195-205
• /
• 2000

도달시간 산정에 있어 고려되어야 할 요소에는 유역형상과 흐름상태, 유출특성 등이 있다. 한편 Singh(1976)은 지표면흐름을 Kinematic Wave 이론으로 해석하였으며 유역형상을 평면형상과 수렴형상, 발산형상으로 분류하여 유역형상에 따른 도달시간 산정식을 제시하였다. 본 논문에서는 Singh이 제시한 평면형상에서의 도달시간 산정식에 기초하여 불투수 지표면 유출에서의 도달시간 산정을 다루었다. 이를 위한 이론식 유도는 흐름을 개수로 흐름(층류, 완난류, 전난류)으로 파악하여 수행하였고 각 개수로 흐름상태에 대한 이론식을 제시하였다. 또한 강우강도를 주요 인자로 고려하여 재현기간별 강우강도를 포함한 근사식 즉 복합형 도달시간 산정식을 개발하였다.

• 대한기계학회논문집A
• /
• 제27권9호
• /
• pp.1562-1570
• /
• 2003

The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. The complexity of these constitutive equations generally requires a stable and accurate numerical method. The radial return mapping is one of the most robust integration scheme currently used. Nonlinear kinematic hardening model of Armstrong-Fredrick type has recovery term and the direction of kinematic hardening increment is not parallel to that of plastic strain increment. In this case, The conventional radial return mapping method cannot be applied directly. In this investigation, we expanded the radial return mapping method to consider the nonlinear kinematic hardening model and implemented this integration scheme into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using Newton method and bisection method. Using dynamic yield condition derived from linearization of flow rule, the integration scheme for elastoplastic and viscoplastic constitutive model was unified. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 춘계 학술대회논문집
• /
• pp.43-48
• /
• 1998

The development of unsteady three-dimensional incompressible viscous solver based on unsteady physical curvilinear coordinate system is presented. A 12-point finite analytic scheme based on local uniform grid spacing is extended for nonuniform grid spacing. The formulation of a condition is suggested to avoid the oscillation of the series summations produced by the application of the method of separation of variables. SIMPLER and pressure Poisson equation techniques are used for solving a velocity-pressure coupled problem. The matrix is solved using the Generalized Minimal RESidual (GMRES) method to enhance the convergence rate of unsteady flow solver and the Kinematic boundary condition of a free surface flow. It is demonstrated that the numerical solutions of these equations are less mesh sensitive.

• Tribology and Lubricants
• /
• 제25권4호
• /
• pp.236-242
• /
• 2009

Oil viscosity is one of the important parameters for machinery condition monitoring. Basically, it is expressed as kinematic viscosity measured by capillary flow and dynamic or absolute viscosity measured by rotary shear viscometry. Recently, acoustic wave techniques appear in the market, measuring viscosity as the product of dynamic viscosity and density. For Newtonian fluids, knowledge of density allows conversion from one viscosity parameter to the other at a specific shear rate and temperature. In this work, oil samples with different chain lengths of viscosity index (VI) improvers and concentrations were examined by different viscometric techniques. Results showed that acoustic viscosity measurements give misleading results for oil samples with high molecular weight VI improvers and at low temperatures ${\leq}40^{\circ}C$.

• Journal of Ship and Ocean Technology
• /
• 제5권1호
• /
• pp.38-54
• /
• 2001

A finite difference method for calculating turbulent flow and surface wave fields around a ship model is evaluated through the comparison with the experimental data of a Series 60 $C_B$=0.6 ship model. The method solves the Reynolds-averaged Navior-Stokes Equations using the non-staggered grid system, the four-stage Runge-Kutta scheme for the temporal integration of governing equations and the Bladwin-Lomax model for the turbulence closure. The free surface waves are captured by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and free-surface conforming grids are generated at each time step so that one of the grid surfaces coincides always with the free surface. The computational results show an overall close agreement with the experimental data and verify that the present method can simulate well the turbulent boundary layers and wakes as well as the free-surface waves.

• 물과 미래
• /
• 제27권2호
• /
• pp.147-154
• /
• 1994

지체시간 내의 누가강우량이 특정 사면경사를 지나게 될 때 발생하는 토석류의 생기조건은 경사면에서의 흐름에 대한 식을 사용하여 구할 수 있으며, 이때의 토석류가 발생하는 한계 누가강우량을 산정하기 위하여 유역의 시스템적 분석 기법이 필요하다. 따라서 본 연구에서는 토석류가 발생하는 한계 누가강유량의 산정 과정을 제시하였으며, 이 방법을 최근 폭발한 운젠 화산의 토석류에 적용하여 신뢰성을 검토하였다. 한편 토석류의 유량은 수리학적 방법인 Kinematic Wave방법을 이용하여 계산하였으며, 이 때의 단면적은 유량에 대한 함수관계에 있다고 가정하였다. 이 방법에 의한 계산치와 사쿠라지마 및 운젠 화산 지역의 관측치의 수문곡선은 대체로 잘 일치하였으며, 토석류의 유출강도를 강우강도와 누가강우량의 조합에 따른 식으로 산정하여 다른 지역에서의 토석류 계산에 이용할 수 있게 하였다.

• Journal of Mechanical Science and Technology
• /
• 제19권12호
• /
• pp.2263-2269
• /
• 2005

The ship wave phenomena in the restricted waterway were investigated by a numerical analysis. The Euler and continuity equations were employed for the present study. The boundary fitted and moving grid system was adopted to enhance the computational efficiency. The convective terms in the governing equations and the kinematic free surface boundary condition were solved by the Constrained Interpolated Profile (CIP) algorithm in order to solve accurately wave heights in far field as well as near field. The advantage of the CIP method was verified by the comparison of the computed results by the CIP and the Maker and Cell (MAC) method. The free surface flow simulation around Wigley hull was performed and compared with the experiment for the sake of the validation of the numerical method. The present numerical scheme was applied to the free surface simulation for various canal sections in order to understand the effect of the sectional shape of waterways on the ship waves. The wave heights on the side wall and the shape of the wave patterns with their characteristics of flow are discussed.