• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2004년도 추계 학술대회논문집
• /
• pp.7-14
• /
• 2004

It is well known that preconditioning methods are efficient for convergence acceleration at compressible low Mach number flows. In this study, the original Euler equations and three preconditioners nondimensionalized differently are implemented in two dimensional inviscid bump flows using the 3rd order MUSCL and DADI schemes as flux discretization and time integration respectively. The multigrid and local time stepping methods are also used to accelerate the convergence. The test case indicates that a properly modified local preconditioning technique involving concepts of a global preconditioning one produces Mach number independent convergence. Besides, an asymptotic analysis for properties of preconditioning methods is added.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제20권3호
• /
• pp.243-259
• /
• 2016

The Richards equation for water movement in unsaturated soil is highly nonlinear partial differential equations which are not solvable analytically unless unrealistic and oversimplifying assumptions are made regarding the attributes, dynamics, and properties of the physical systems. Therefore, conventionally, numerical solutions are the only feasible procedures to model flow in partially saturated porous media. The standard Finite element numerical technique is usually coupled with an Euler time discretizations scheme. Except for the fully explicit forward method, any other Euler time-marching algorithm generates nonlinear algebraic equations which should be solved using iterative procedures such as Newton and Picard iterations. In this study, lumped mass and distributed mass in the frame of Picard and Newton iterative techniques were evaluated to determine the most efficient method to solve the Richards equation with finite element model. The accuracy and computational efficiency of the scheme and of the Picard and Newton models are assessed for three test problems simulating one-dimensional flow processes in unsaturated porous media. Results demonstrated that, the conventional mass distributed finite element method suffers from numerical oscillations at the wetting front, especially for very dry initial conditions. Even though small mesh sizes are applied for all the test problems, it is shown that the traditional mass-distributed scheme can still generate an incorrect response due to the highly nonlinear properties of water flow in unsaturated soil and cause numerical oscillation. On the other hand, non oscillatory solutions are obtained and non-physics solutions for these problems are evaded by using the mass-lumped finite element method.

• 한국전산구조공학회논문집
• /
• 제21권5호
• /
• pp.475-482
• /
• 2008

본 논문은 로터 시스템의 디스크 회전운동을 표현하는데 있어 운동방정식을 통합하는 과정에서 기존 연구자들이 채택한 오일러 각 사용법이 일관성이 없음을 지적하였다. 기존 연구자들은 오일러 각 순서가 달라서 속도와 운동에너지도 달리 산정하였음은 물론, 운동방정식은 오직 선형 시스템만 취급해 왔다 이러한 오일러 각 사용법의 단점을 극복하기 위하여 회전운동을 더욱 단순하게 매개화할 수 있는 4원법(quaternion)과 구 좌표계를 적용하여 비선형 시스템을 도출하였다. 이를 바탕으로 수치해석을 통하여 기존 방법과 비교하여 제안한 방법의 신뢰성과 우수성을 보였다.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2002년도 춘계학술대회논문집
• /
• pp.435-440
• /
• 2002

The free vibration of partially embedded piles is investigated. The pile model is based on the Bernoulli-Euler beam theory and the soil is idealized as a Winkler model for mathematical simplicity. The governing differential equation for the free vibrations of such members is solved numerically The piles with one typical end constraint (clamped/hinged/free) and the other hinged end with rotational spring are applied in numerical examples. The lowest three natural frequencies are calculated over a range of non-dimensional system parameters: the rotational spring parameter, the relative stiffness and the embedded ratio.

• 대한조선학회지
• /
• 제25권2호
• /
• pp.1-10
• /
• 1988

An axially marching numerical method is developed for the simulation of the internal waves produced by translation of a submersed vehicle in a density-stratified ocean. The method provides for the direct solution of the primitive variables [$\upsilon,\;p,\;\rho$] for the nonlinear and steady state three-dimensional Euler's equation with a non-constant density term in the vehicle-fixed cartesian co-ordinate system. By utilizing a known potential flow around the vehicle for an estimate of the axial velocity gradient, the present parabolic algorithm local upstreamwise disturbances and axial velocity variation.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2006년도 제27회 추계학술대회논문집
• /
• pp.277-280
• /
• 2006

천음속 영역에서 순항 속도로 비행하는 상용 여객기에 대해 공력을 해석하였다. 실제 여객기의 축소모형을 사진으로 계측하여 역설계하여 형상을 모델링하였다. 사면체로 전체 격자계를 구성하였고 정상 상태에서 압축성 비점성 3차원 Euler 방정식을 풀었다. 날개에 엔진이 장착되었을 때와 제거되었을 때에 대해 계산하였는데 엔진이 장착되면 양력과 항력이 각각 1.49%와 3.9%씩 증가하였다. 또한 날개 표면에서 압력 분포와 엔진 주위의 유동 특성을 밝혔다.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 추계학술대회논문집B
• /
• pp.570-577
• /
• 2001

Three-dimensional time-dependent flow past a circular cylinder is numerically investigated using direct numerical simulation for Reynolds number 280 and 300. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme and the pressure and viscosity terms are applied by the 4th order central difference. And in spanwise, Navier-Stokes equation is distributed using of Spectral Method. At Reynolds number 259 the two-dimensional wake becomes linearly unstable to a second branch of modes with wavelength about 1.0 diameters at onset (B-mode). Present results of three-dimensional effects of in wake of a circular cylinder is represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• Journal of Mechanical Science and Technology
• /
• 제15권5호
• /
• pp.671-680
• /
• 2001

Boundary condition is one of the major factors to influence the numerical stability and solution accuracy in numerical analysis. One of the most important physical boundary conditions in the flowfield analysis is the wall boundary condition imposed on the body surface. To solve a two-dimensional Euler equation, totally four numerical wall boundary conditions should be prescribed. Two of them are supplied by the flow tangency condition. The other two conditions, therefore, should be prepared additionally in a suitable way. In this paper, four different sets of wall boundary conditions are proposed and then applied to solve high-speed flowfields around a quarter circle geometry. A two-dimensional compressible Euler solver is prepared based on the finite volume method. This solver hires three different upwind schemes; Steger-Warmings flux vector splitting, Roes flux difference splitting, and Lious advection upstream splitting method. It is found that the way to specify the additional numerical wall boundary conditions strongly affects the overall stability and accuracy of the upwind schemes in high-speed flow calculation. The optimal wall boundary conditions should be also chosen very carefully depending on the numerical schemes used to solve the problem.

• 한국운동역학회지
• /
• 제16권4호
• /
• pp.49-59
• /
• 2006

The purpose of this study was to investigate the relations between the segments of the body, the three dimensional anatomical angle according to sweep shot in ice hockey. The subjects of this study were five professional ice hockey players. The reflective makers were attached on anatomical boundary line of body. For the movement analysis three dimensional cinematographical method(APAS) was used and for the calculation of the kinematic variables a self developed program was used with the LabVIEW 6.1 graphical programming(Johnson, 1999) program. By using Eular's equations the three dimensional anatomical Cardan angles of the joint and ice hockey stick were defined. 1. In three dimensional linear velocity of blade the Y axis showed maximum linear velocity almost impact, the X axis(horizontal direction) and the Z axis(vertical direction) maximum linear velocity of blade did not show at impact but after impact this will resulted influence upon hitting puck. 2. The resultant linear velocity of each segment of right arm showed maximum resultant linear velocity at impact. It could be suggest that the right arm swing patterns is kind of push-like movement. therefore the upper arm is the most important role in the right arm swing. 3. The three dimensional anatomical angular displacement of trunk in flexion-extension showed flexion all around the wrist shot. The angular displacement of trunk in internal-external rotation showed internal rotation angle at the backswing top and and increased the angle after the impact. while there is no significant adduction-abduction. 4. The three dimensional anatomical angular displacement of trunk showed most important role in wrist shot. and is follwed by shoulder joints, in addition the movement of elbow/wrist joints showed least to the shot. this study result showed upperlimb of left is more important role than upperlimb of right.

• Structural Engineering and Mechanics
• /
• 제58권1호
• /
• pp.103-119
• /
• 2016

Multi-storey frame structures are frequently exposed to static and dynamic forces. Therefore analyses of static (buckling) and dynamic stability come into prominence for these structures. In this study, the effects of number of storey, static and dynamic load parameters, crack depth and crack location on the in-plane static and dynamic stability of cracked multi-storey frame structures subjected to periodic loading have been investigated numerically by using the Finite Element Method. A crack element based on the Euler beam theory is developed by using the principles of fracture mechanics. The equation of motion for the cracked multi-storey frame subjected to periodic loading is achieved by Lagrange's equation. The results obtained from the stability analysis are presented in three dimensional graphs and tables.