• 연구논문집
• /
• s.27
• /
• pp.15-34
• /
• 1997

The 3-dimensional unsteady compressible flows around the high speed train have been simulated for the train entering a tunnel and for passing another train. The simulation method employs the implicit approximation-factorization finite difference algorithm for the inviscid Euler equations in general curvilinear coordinates. A moving grid scheme is applied in order to resolve the train movement relative to the tunnel and the other train. The velo-city and pressure fields and pressure drag are calculated to study the effects of tunnel and the other train. The side directional force which is time dependent is also computed for the passing train. Pressure distribution shows that the compression wave is generated in front of the train noise just after the tunnel entrance and proceeds along the inside of tunnel.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.04a
• /
• pp.5-29
• /
• 1995

A three-dimensional flow field induced by two trains passing by each other inside a tunnel is studied based on the numerical simulation of the three-dimensional compressible Euler/Navier-Stokes equations formulated in the finite difference approximation. Domain decomposition method with the FSA(fortified solution algorithm) interface scheme is used to treat this moving-body problem. The computed resluts show basic characteristic of the flow field created when two trains passing by each other. History of the pressure distributions and the aerodynamic forces acting on the trains are mailnly discussed. The results indicate that the phenomenon is complicated due to the interaction of the flow induced by two trains. Strong side force occurs between the two trains when the front portion of the opposite train passes by. It fluctuates rapidly and maximum suction force occurs when two trains are aligned side by side. The results also indicate the effectiveness of the present numerical method for moving boundary problems.

• Proceedings of the KIEE Conference
• /
• 1997.07b
• /
• pp.474-477
• /
• 1997

In this paper, an efficient attitude determination algorithm based on the Kalman Filter which combines earth/sun sensor data with gyro data in a mutually compensating manner is presented. Quaternion is used as the attitude state to save computation time and to prevent the gimbal-lock situation associated with Euler angles. Gyro data allows the use of the kinematic equation instead of space vehicle's dynamic equation which is usually based on approximation of the actual dynamics and inaccurate torque information. The gyro data are used to propagate the attitude through kinematic equation and the earth/sun sensor data are used to update the attitude and estimate the gyro bias. Simulation results for the KOMPSAT attitude determination system are presented.

• Journal of Ship and Ocean Technology
• /
• v.4 no.2
• /
• pp.1-12
• /
• 2000

The present study concentrates on the weak-scatterer hypothesis for the nonlinear wave-body interaction problems. In this method, the free surface boundary conditions are linearized on the incoming wave profile and the exact body motion is applied. The considered problems are the diffraction problem near a circular cylinder and the ship response in oblique waves. The numerical method of solution is a Rankine panel method. The Rankine panel method of this study adopts the higher-order B spline basis function for the approximation of physical variables. A modified Euler scheme is applied for the time stepping, which has neutral stability. The computational result shows some nonlinear behaviors of disturbance waves and wave forces. Moreover, the ship response shows very close results to experimental data.

• Proceedings of the KIEE Conference
• /
• 2003.07a
• /
• pp.107-109
• /
• 2003

전력계통의 과도 안정도 해석의 접근방법에는 SI(Simultaneous Implicit)법과 PE(Partitioned Explicit)법 두 가지방법을 사용해오고 있다. SI법에는 Trapezoidal법 등이 있고, PE법에는 Runge-Kutta법, Euler법등이 사용되고 있다. SI법인 Trapezoidal법은 PE법의 Runge-Kutta법 또는 Euler법에 비해 시간간격을 크게 해서 계산속도를 줄일 수 있다는 장점이 있지만, 정화도면에서는 신뢰한 수 없는 단점이 있다. 이 논문에서는 포물선 사법을 이용하여 Trapezoidal법의 정확도를 개선학 수 있는 방법을 제시하고 명확한 수학적 증명을 통해 타당성을 보여준다. 연속함수와 불연속함수에 대해서 Runge-Kutta법과 Trapezoidal법과 제안한 방법을 적용시켜서 제안한 방법의 정화함을 보여준다.

• Journal of the Korean Mathematical Society
• /
• v.60 no.2
• /
• pp.273-302
• /
• 2023

In this paper, a fully discrete numerical scheme for the viscoelastic Oldroyd flow is considered with an introduced auxiliary variable. Our scheme is based on the finite element approximation for the spatial discretization and the backward Euler scheme for the time discretization. The integral term is discretized by the right trapezoidal rule. Firstly, we present the corresponding equivalent form of the considered model, and show the relationship between the origin problem and its equivalent system in finite element discretization. Secondly, unconditional stability and optimal error estimates of fully discrete numerical solutions in various norms are established. Finally, some numerical results are provided to confirm the established theoretical analysis and show the performances of the considered numerical scheme.

• Journal of Ocean Engineering and Technology
• /
• v.37 no.6
• /
• pp.256-265
• /
• 2023

Many numerical methods have been developed since 1961, but unresolved issues remain. This study developed a numerical method to address these issues and determine the coefficients and properties of rotational waves with a shear current in a finite water depth. The number of unknown constants was reduced significantly by introducing a wavelength-independent coordinate system. The reference depth was calculated independently using the shooting method. Therefore, there was no need for partial derivatives with respect to the wavelength and the reference depth, which simplified the numerical formulation. This method had less than half of the unknown constants of the other method because Newton's method only determines the coefficients. The breaking limit was calculated for verification, and the result agreed with the Miche formula. The water particle velocities were calculated, and the results were consistent with the experimental data. Dispersion relations were calculated, and the results are consistent with other numerical findings. The convergence of this method was examined. Although the required series order was reduced significantly, the total error was smaller, with a faster convergence speed.

• Structural Engineering and Mechanics
• /
• v.75 no.3
• /
• pp.283-294
• /
• 2020

In general, the study of a high-rise building's behaviour when subjected to a horizontal load (wind or earthquake) is carried out through numerical modelling with finite elements method. This paper proposes a new, original approach based on the use of a multi-beams model. By redistributing bending and axial stiffness of horizontal elements (beams and slabs) along vertical elements, it becomes possible to produce a system of differential equations able to represent the structural behaviour of the whole building. In this paper this approach is applied to the study of bending behaviour in a 37-storey building (Torre Pontina, Latina, Italy) with a regular reinforced concrete structure. The load considered is the wind, estimated in accordance with Italian national technical rules and regulations. To simplify the explanation of the approach, the wind load was considered uniform on the height of building with a value equal to the average value of the wind load distribution. The system of differential equations' is assessed numerically, using Matlab, and compared with the obtainable solution from a finite elements model along with the obtainable solutions via classical Euler-Bernoulli beam theory. The comparison carried out demonstrates, in the case study examined, an excellent approximation of structural behaviour.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.249-257
• /
• 2007

A high resolution numerical method aimed at solving gas-liquid two-phase flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media and some data related to computational efficiency are made. Comparisons of predicted results and exact solutions are provided and discussed.

• Journal of computational fluids engineering
• /
• v.12 no.2
• /
• pp.53-62
• /
• 2007

A high resolution numerical method aimed at solving cavitating flow is proposed and applied to gas-liquid two-phase shock tube problem. The present method employs a finite-difference 4th-order Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL TVD scheme. By applying the homogeneous equilibrium cavitation model, the present density-based numerical method permits simple treatment of the whole gas-liquid two-phase flow field, including wave propagation and large density changes. The speed of sound for gas-liquid two-phase media is derived on the basis of thermodynamic relations and compared with that by eigenvalues. By this method, a Riemann problem for Euler equations of one dimensional shock tube was computed. Numerical results such as detailed observations of shock and expansion wave propagations through the gas-liquid two-phase media at isothermal condition and some data related to computational efficiency are made. Comparisons of predicted results and exact solutions are provided and discussed.