• Journal of computational fluids engineering
• /
• v.16 no.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
• /
• v.31 no.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.

• Journal of the Chungcheong Mathematical Society
• /
• v.25 no.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.

• Journal of Ocean Engineering and Technology
• /
• v.22 no.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.

• Journal of the Society of Naval Architects of Korea
• /
• v.57 no.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
• /
• v.21 no.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.

• Journal of the Korean earth science society
• /
• v.38 no.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.

• Journal of Ocean Engineering and Technology
• /
• v.20 no.6 s.73
• /
• 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.

• Journal of Power System Engineering
• /
• v.14 no.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.

• Journal of Power System Engineering
• /
• v.16 no.5
• /
• pp.40-46
• /
• 2012

For free surface flow problem, a high-order spectral/boundary element method is 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. Using the combination of these two methods, the free surface flow problems of a submerged moving body are solved in time domain. In the present study, lifting surface theory is added to the former work to include effects of lift force. Therefore, a new formulation for the basic mathematical theory is introduced to contain the lift body in calculation.