• 한국전산유체공학회지
• /
• 제16권4호
• /
• pp.72-83
• /
• 2011

A high-order discontinuous Galerkin method for the two-dimensional compressible Navier-Stokes equations was developed on unstructured triangular meshes. For this purpose, the BR2 methd(the second Bassi and Rebay discretization) was adopted for space discretization and an implicit Euler backward method was used for time integration. Numerical tests were conducted to estimate the convergence order of the numerical solutions of the Poiseuille flow for which analytic solutions are available for comparison. Also, the flows around a flat plate, a 2-D circular cylinder, and an NACA0012 airfoil were numerically simulated. The numerical results showed that the present implicit discontinuous Galerkin method is an efficient method to obtain very accurate numerical solutions of the compressible Navier-Stokes equations on unstructured meshes.

• Wind and Structures
• /
• 제31권5호
• /
• pp.473-482
• /
• 2020

The first-order method for estimating the extreme wind pressure on building envelopes with consideration of the directionality of wind speed and wind pressure is improved to enhance its computational efficiency. In this improved method, the result is obtained directly from the empirical distribution of a random selection of annual maximum wind pressure samples generated by a Monte Carlo method, rather than from the previously utilized extreme wind pressure probability distribution. A discussion of the relationship between the first- and full-order methods indicates that when extreme wind pressures in a non-typhoon climate with a high return period are estimated with consideration of directionality, using the relatively simple first-order method instead of the computationally intensive full-order method is reasonable. The validation of this reasonableness is equivalent to validating two assumptions to improve its computational efficiency: 1) The result obtained by the full-order method is conservative when the extreme wind pressure events among different sectors are independent. 2) The result obtained by the first-order method for a high return period is not significantly affected when the extreme wind speeds among the different sectors are assumed to be independent. These two assumptions are validated by examples in different regions and theoretical derivation.

• 충청수학회지
• /
• 제25권2호
• /
• pp.359-365
• /
• 2012

The purpose of this paper is to provide some corrections regarding algebraic flaws encountered in the paper entitled "Sixteenth-order method for nonlinear equations" which was published in January of 2010 by Li et al.[9]. Further detailed comments on their error equation are stated together with convergence analysis as well as high-precision numerical experiments.

• 한국해양공학회지
• /
• 제22권5호
• /
• pp.25-30
• /
• 2008

A B-spline based high order panel method was developed for the motion of bodies in an ideal fluid, either of infinite extent or with a free boundarysurface. In this method, both the geometry and the potential are represented by the B-spline, which guarantees more accurate results than most potential based low order methods. In the present work, we applied this B-spline based high order method to the radiation problem of floating bodies. The boundary condition on the free surface was satisfied by adopting a Kelvin-type Green function and irregular frequencies were removed by placing additional control points on the free surface surrounding the body. The numerical results were validated by comparison with existing numerical and experimental results.

• 대한조선학회논문집
• /
• 제57권2호
• /
• pp.114-123
• /
• 2020

A high-order potential-based panel method based on Green's theorem, with piecewise-linear dipole strength on triangular panels, is formulated for the analysis of potential flow around a three-dimensional wing. Previous low-order panel methods adopt square panels with piecewise-constant dipole strength, which results in inherent errors. Square panels can not represent a high curvature lifting body, such as propellers, since the four vertices of the square panel do not locate at the same flat plane. Moreover the piecewise-constant dipole strength induces inevitable errors due to the steps in dipole strength between adjacent panels. In this paper a high-order panel method is formulated to improve accuracy by adopting a piecewise linear dipole strength on triangular panels. Firstly, the square panels are replaced by triangular panels in order to increase the geometric accuracy in representing the shape of the object with large curvature. Next, the step difference of the dipole strength between adjacent panels is removed by adopting piecewise-linear dipole strength on the triangular panels. The calculated results by the present method is compared with analytical ones for simple non-lifting geometries, such as ellipsoid. The results for an elliptic wing with zero thickness at finite angle of attack are compared with Jordan's results. The comparison shows reasonable agrements for the both lifting and non-lifting bodies.

• Steel and Composite Structures
• /
• 제21권4호
• /
• pp.755-773
• /
• 2016

In the present paper, a new high-order simple-input analytical method is used to study thick laminated composite straight tubes subjected to combined axial force, torque and bending moment. The most general displacement field of elasticity for an arbitrary laminated composite straight tube is obtained to analytically calculate stresses under combined loadings based on a layerwise method. The accuracy of the proposed method is subsequently verified by comparing the numerical results obtained using the proposed method with finite element method (FEM) and experimental data. The results show good corresponded. The proposed method provides advantages in terms of computational time compared to FEM.

• 한국지구과학회지
• /
• 제38권2호
• /
• pp.105-114
• /
• 2017

Two dimensional finite element method with quadrilateral basis functions was applied to the spherical high order filter on the spherical surface limited area domain. The basis function consists of four shape functions which are defined on separate four grid boxes sharing the same gridpoint. With the basis functions, the first order derivative was expressed as an algebraic equation associated with nine point stencil. As the theory depicts, the convergence rate of the error for the spherical Laplacian operator was found to be fourth order, while it was the second order for the spherical Laplacian operator. The accuracy of the new high order filter was shown to be almost the same as those of Fourier finite element high order filter. The two-dimension finite element high order filter was incorporated in the weather research and forecasting (WRF) model as a hyper viscosity. The effect of the high order filter was compared with the built-in viscosity scheme of the WRF model. It was revealed that the high order filter performed better than the built in viscosity scheme did in providing a sharper cutoff of small scale disturbances without affecting the large scale field. Simulation of the tropical cyclone track and intensity with the high order filter showed a forecast performance comparable to the built in viscosity scheme. However, the predicted amount and spatial distribution of the rainfall for the simulation with the high order filter was closer to the observed values than the case of built in viscosity scheme.

• 한국해양공학회지
• /
• 제20권6호
• /
• pp.67-74
• /
• 2006

A high-order spectral/boundary-element method is newly adapted as an efficient numerical tool. This method is one of the most efficient numerical methods by which the nonlinear gravity waves can be simulated and hydrodynamic forces also can be calculated in time-domain. In this method, the velocity potential is expressed as the sum of surface potential and body potential. Then, surface potential is solved by using the high-order spectral method and body potential is solved by using the high-order boundary element method. By the combination of these two methods, the wave-making problems by a submerged sphere oscillating with large amplitude under the free~surface are solved in time-domain. Through the example calculations, nonlinear effects on free-surface profiles and hydrodynamic forces are shown and discussed.

• 동력기계공학회지
• /
• 제14권6호
• /
• pp.75-82
• /
• 2010

In this study, the topics for free-surface wave simulation, nonlinear hydrodynamic force, and the critical resonance frequency of so-called ${\tau}=U{\omega}/g$=1/4 are discussed. A high-order spectral/boundary element method is newly adapted as an efficient numerical tool. This method is one of the most efficient numerical methods by which the nonlinear gravity waves can be simulated and hydrodynamic forces also can be calculated in time domain. In this method, the velocity potential is expressed as the sum of surface potential and body potential. Then, surface potential is solved by using the high-order spectral method and body potential is solved by using the high-order boundary element method. By the combination of these two methods, the wave-making problems by a submerged sphere oscillating with forward speed under the free-surface are solved in time domain.

• 동력기계공학회지
• /
• 제16권5호
• /
• pp.40-46
• /
• 2012

본 연구에서는 자유표면 유동문제를 효율적으로 계산하기 위한 방법으로서, 수치 파수조를 구현하여 잠수체에 의한 조파현상을 시간영역에서 다룰 수 있는 수치해법에 대하여 소개하였다. 그리고 이를 이용하여 수행된 연구내용들을 검토하고 양력물체의 경우를 포함하기 위하여 개선된 수학적 정식화 및 수치해법의 개요를 다루었다. 임의의 운동을 하는 양력물체에 의한 조파현상을 전산기로 구현하는 수치 파수조는 과중한 계산시간이 문제가 되는데, 이는 수치 Kutta 조건의 구현과 양력표면 후방의 wake 영역을 계산에서 고려해야 하기 때문이다. 따라서 한층 더 수치계산의 효율성이 중요하다고 판단되므로, 본 연구에서 소개된 3차원 고차 스펙트럴/경계요소법(High-Order Spectral/Boundary Element Method)은 자유표면 요소수를 N이라 할 때 그 계산량이 NlogN에 비례(N이 클 때는 거의 선형적으로 비례)하여 증가하므로 기존의 방법들 보다 매우 효율적인 수치해법이라 할 수 있다. 향후, 본 연구의 타당성 검증을 위한 수치코드의 개선과 여러 가지 수치계산결과 비교 등의 노력이 더 필요하다고 생각된다.