• Journal of applied mathematics & informatics
• /
• 제24권1_2호
• /
• pp.17-30
• /
• 2007

A new scheme for solving the Vlasov equation using a compactly supported wavelets basis is proposed. We use a numerical method which minimizes the numerical diffusion and conserves a reasonable time computing cost. So we introduce a representation in a compactly supported wavelet of the derivative operator. This method makes easy and simple the computation of the coefficients of the matrix representing the operator. This allows us to solve the two equations which result from the splitting technique of the main Vlasov equation. Some numerical results are exposed using different numbers of wavelets.

• Journal of applied mathematics & informatics
• /
• 제9권1호
• /
• pp.1-13
• /
• 2002

In this paper we present numerical methods fur computations of nonequilibrium hypersonic flow of air around bodies including chemical reaction effects and present numerical result of the flow over concave corners. We developed implicit finite difference method to overcome numerical difficulties with the lack of resolution behind the shock and near the body. Using our method we were able to find details of the flow properties near the shock and body and were able to continue the computation of the flow for a long distance from the corner of the body.

• 대한조선학회논문집
• /
• 제33권1호
• /
• pp.1-8
• /
• 1996

3차원 날개면이 부착된 복합지지형 고속선의 정상상태에서의 조파저항을 포텐셜 기저 판요소법으로 해석하였다. 계산에 사용된 고속선은 물 속에 잠겨있는 몸체와 몸체의 중간 및 후방에 붙어 있는 3차원 날개면과 수면을 관통하는 앞, 뒤의 스트럿트로 구성되었다. 물체 표면(몸체, 날개면 및 스트럿트)에는 쏘오스와 다이폴을, 자유표면에는 쏘오스를 분포하였고, 선형화된 자유표면 조건과 방사조건을 만족시키기 위해 4점 유한차분을 이용하였으며, 날개면의 유동해석을 위해 준압력 Kutta조건(semi-linear pressure Kutta condition)을 적용하였다. 패널사이의 틈새 문제를 개선하기 위해 비 평면성이 고려된 쌍곡면 판요소법을 적용하여 각 선체 표면에서의 수치 계산 정도를 높이고자 하였다. 수치 계산 결과는 회류수조에서 모형 시험을 수행한 결과와 비교하였으며, 이로부터 본 연구에서 개발된 수치 계산법은 고속선의 최적 선형 개발에 이용 가능한 도구가 될 수 있는 것으로 나타났다.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 춘계학술대회논문집
• /
• pp.75-80
• /
• 2009

The present paper focuses on the analysis of aero-acoustics characteristic by several boundary conditions. In this simulation, a high-order and high-resolution numerical schemes are used for the accurate computation of compressible flow with several boundary conditions including characteristic boundary conditions as well as extrapolation and zonal characteristic boundary condition. These boundary conditions are applied to the computation of two dimensional circular cylinder flows with Mach number of 0.3 and Reynolds number of 400. The computation results are validated with measurement datum and other computation results for the Strouhal frequency of vortex shedding, the mean drag coefficient and root-mean-square lift for the unsteady periodic flow regime. Secondary frequency is predicted by three kinds of boundary conditions characteristic.

• 한국전산유체공학회지
• /
• 제14권3호
• /
• pp.25-32
• /
• 2009

The present paper focuses on the analysis of aero-acoustics characteristic by appling different four boundary conditions. The high-order and high-resolution numerical schemes are used for discrete accurate computation of compressible flow. The four boundary conditions include extrapolation, characteristic boundary condition, zonal characteristic boundary condition. These boundary conditions are applied to the computation of two dimensional circular cylinder flows with Mach number of 0.3 and Reynolds number of 400. The computation results are validated against measurement data and other computation results for the Strouhal frequency of vortex shedding, the mean drag coefficient and root-mean-square lift for the unsteady periodic flow regime. The characteristics of secondary frequency is predicted by three kinds of boundary conditions.

• 한국소음진동공학회논문집
• /
• 제12권9호
• /
• pp.674-684
• /
• 2002

In order to decrease remarkably the computation time and storage used in the direct integration method without the loss of accuracy, authors suggest a new transient analysis algorithm. This algorithm is derived from the combination of three techniques, that is, the transfer technique of the transfer stiffness coefficient method, the modeling technique of the finite element method, and the numerical integration technique of the Newmark method. In this paper, the transient analysis algorithm of a frame structure is formulated by the proposed method. The accuracy and computation efficiency of the proposed method are demonstrated through the comparing with the computation results by the direct integration method for three computation models under various excitations.

• 한국정밀공학회지
• /
• 제12권9호
• /
• pp.137-147
• /
• 1995

The unbalance response analysis is one of the essential area in the forced vibration analysis of rotor-bearing systems. Local bearing parameters in rotor-bearing systems are the major sources which give rise to a difficulty in unbalance response computation due to the complicated dynamic properties such as rotational speed dependency and anisotropy. In the present paper, an efficient method for unbalance responses is proposed so as to easily take into account bearing parameters in computation. An exact matrix condensation procedure is proposed which enables the present method to compute unbalance responses by dealing with condensed, small matrices. The proposed method causes no errors even though the computation procedure is based on the small matrices condensed from the full matrices. The present method is illustrated through a numerical example and compared with the conventional method.

• Journal of Advanced Research in Ocean Engineering
• /
• 제3권4호
• /
• pp.193-203
• /
• 2017

This study investigated the steady-flow effects present in the numerical computation of the resistance added to a ship in waves. For a ship advancing in the forward direction, a time-domain 3D Rankine panel method is applied to solve the ship motion problem, and the added resistance due to waves is calculated using a near-field method, with the direct integration of the second-order pressure on the hull surface. In the linear potential theory, the steady flow is approximated by the basis potential of a uniform flow or double-body flow in order to linearize the boundary conditions. By applying these two different linearization schemes, the coupling effects between steady and unsteady solutions were examined. Furthermore, in order to analyze the steady-flow effects on the hull geometry, the computation results for two realistic hull forms, a KVLCC2 tanker and DTC containership, were compared. In particular, the mj term, which represents the coupling effects under the body boundary condition, was evaluated considering the geometry of a non-wall-sided ship. Lastly, the characteristics of the linearization schemes were examined in relation to the disturbed waves around a ship and the components of added resistance.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 1996년도 자동차부품 제작기술의 진보
• /
• pp.97-109
• /
• 1996

In the present work computations are carried out for analysis of complicated sheet metal forming process such as forming of a rear hinge. Finite element formulation using dynamic explicit time integration scheme and step-wise combined Implicit/Explicit scheme are introduced for numerical analysis of sheet metal forming process. The rigid-plastic finite element method based on membrane elements has long been employed as a useful numerical technique for the analysis of sheet metal forming because of its time effectiveness. The explicit scheme in general use is based on the elastic-plastic modelling of material requiring large computation time. In finite element simulation of sheet metal forming processes, the robustness and stability of computation are important requirements since the computation time and convergency become major points of consideration besides the solution accuracy due to the complexity of geometry and boundary conditions. The implicit scheme employs a more reliable and rigorous scheme in considering the equilibrium at each step of deformation, while in the explicit scheme the problem of convergency is eliminated at the cost of solution accuracy. The explicit approach and the implicit approach have merits and demerits, respectively. In order to combine the merits of these two methods a step-wise combined implicit/explicit scheme has been developed.

• Structural Engineering and Mechanics
• /
• 제38권6호
• /
• pp.803-823
• /
• 2011

In this paper, the three dimensional Parallel Realistic Failure Process Analysis ($RFPA^{3D}$-Parallel) code based on micromechanical model is employed to investigate the bonding behavior in FRP sheet bonded to concrete in single shear test. In the model, the heterogeneity of brittle disordered material at a meso-scale was taken into consideration in order to realistically demonstrate the mechanical characteristics of FRP-to-concrete. Modified Mohr-coulomb strength criterion with tension cut-off, where a stressed element can damage in shear or in tension, was adopted and a stiffness degradation approach was used to simulate the initiation, propagation and growth of microcracks in the model. In addition, a Master-Slave parallel operation control technique was adopted to implement the parallel computation of a large numerical model. Parallel computational results of debonding of FRP-concrete visually reproduce the spatial and temporal debonding failure progression of microcracks in FRP sheet bonded to concrete, which agrees well with the existing testing results in laboratory. The numerical approach in this study provides a useful tool for enhancing our understanding of cracking and debonding failure process and mechanism of FRP-concrete and our ability to predict mechanical performance and reliability of these FRP sheet bonded to concrete structures.