• Structural Engineering and Mechanics
• /
• v.22 no.1
• /
• pp.1-16
• /
• 2006

This paper presents a virtual boundary element-equivalent collocation method (VBEM) for the plane magnetoelectroelastic solids, which is based on the fundamental solutions of the plane magnetoelectroelastic solids and the basic idea of the virtual boundary element method for elasticity. Besides all the advantages of the conventional boundary element method (BEM) over domain discretization methods, this method avoids the computation of singular integral on the boundary by introducing the virtual boundary. In the end, several numerical examples are performed to demonstrate the performance of this method, and the results show that they agree well with the exact solutions. So the method is one of the efficient numerical methods used to analyze megnatoelectroelastic solids.

• The Bulletin of The Korean Astronomical Society
• /
• v.44 no.1
• /
• pp.46.1-46.1
• /
• 2019

We present an accurate and efficient method to calculate the gravitational potential of an isolated system in three-dimensional Cartesian and cylindrical coordinates subject to vacuum (open) boundary conditions. Our method consists of two parts: an interior solver and a boundary solver. The interior solver adopts an eigenfunction expansion method together with a tridiagonal matrix solver to solve the Poisson equation subject to the zero boundary condition. The boundary solver employs James's method to calculate the boundary potential due to the screening charges required to keep the zero boundary condition for the interior solver. A full computation of gravitational potential requires running the interior solver twice and the boundary solver once. We develop a method to compute the discrete Green's function in cylindrical coordinates, which is an integral part of the James algorithm to maintain second-order accuracy. We implement our method in the {\tt Athena++} magnetohydrodynamics code, and perform various tests to check that our solver is second-order accurate and exhibits good parallel performance.

• Journal of Advanced Research in Ocean Engineering
• /
• v.3 no.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.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.430-435
• /
• 2001

This paper introduces a pressure correction method for microflow computation. Conventional CFD methods with no slip boundary condition fail to predict the rarefaction effect of the wall when simulating gas microflows in the slip-flow regime. Pressure correction method with an appropriate slip boundary condition is an efficient tool in analyzing microscale flows. The present unstructured SIMPLE algorithm adopts both the classical Maxwell boundary condition and Langmuir boundary condition proposed by Myong. The simulation results of microchannel flows show that the proposed method has an effective predictive capability for microscale flows.

• Journal of computational fluids engineering
• /
• v.18 no.2
• /
• pp.56-65
• /
• 2013

Two compression ramp problems and an impinging shock problem are computed to investigate influence of turbulence models and eddy viscosity on the shock-wave / boundary layer interaction. A Navier-Stokes boundary layer generation code was applied to the generation of inflow boundary conditions. Computational results are validated well with the experimental data and effects of turbulence models are investigated. It is shown that the behavior of turbulence (eddy) viscosity directly affects both the extent of the separation and shock-wave positions over the separation.

• Journal of computational fluids engineering
• /
• v.5 no.1
• /
• pp.43-59
• /
• 2000

Characteristic boundary conditions are discussed in conjunction with a flux-difference splitting formulation as modified from Roe's linearization. Details of how one can implement the characteristic boundary conditions which are made compatible with the interior point formulation are described for different types of boundaries including subsonic outflow and adiabatic wall. The validity of boundary conditions are demonstrated through computation of transonic airfoil, supersonic ogive-cylinder, hypersonic cylinder, and S-duct internal flows. The computed wall pressure distributions are compared with published experimental and computed data. Objectives of this paper are thus to give insight of formulation procedure of a flux-difference splitting method and to pave ways for other users to adopt present boundary procedure on their numerical methods.

• Journal of applied mathematics & informatics
• /
• v.28 no.3_4
• /
• pp.813-821
• /
• 2010

In this paper, the tenth-order linear boundary value problems are solved using reproducing kernel method. The algorithm developed approximates the solutions, and their higher-order derivatives, of differential equations and it avoids the complexity provided by other numerical approaches. First a new reproducing kernel space is constructed to solve this class of tenth-order linear boundary value problems; then the approximate solutions of such problems are given in the form of series using the present method. Three examples compared with those considered by Siddiqi, Twizell and Akram [S.S. Siddiqi, E.H. Twizell, Spline solutions of linear tenth order boundary value problems, Int. J. Comput. Math. 68 (1998) 345-362; S.S.Siddiqi, G.Akram, Solutions of tenth-order boundary value problems using eleventh degree spline, Applied Mathematics and Computation 185 (1)(2007) 115-127] show that the method developed in this paper is more efficient.

• Journal of applied mathematics & informatics
• /
• v.26 no.1_2
• /
• pp.45-60
• /
• 2008

Here we consider the evaluation of the the dynamic component of the second order force due to wave diffraction by a circular cylinder analytically and numerically. The cylinder is fixed, vertical, surface piercing in water of finite uniform depth. The formulation of the wave-structure interaction is based on the assumption of a homogeneous, ideal, incompressible, and inviscid fluid. The nonlinearity in the wave-structure interaction problem arises from the free surface boundary conditions, namely, dynamic and kinematic free surface boundary conditions. We expand the velocity potential and free surface elevation functions in terms of a small parameter and then consider the second order diffraction problem. After deriving the pressure using Bernoulli's equation, we obtain the analytical expression for the dynamic component of the second order force on the cylinder by integrating the pressure over the wetted surface. The computation of the dynamic force component requires only the first order velocity potential. Numerical results for the dynamic force component are presented.

• Journal of computational fluids engineering
• /
• v.6 no.4
• /
• pp.15-25
• /
• 2001

Numerical prediction of the diffusion controlled transition in a turbine gas pass is important because it can change the local heat transfer rate over a turbine blade as much as three times. In this study, the gas flow over turbine blade is simplified to the flat plate boundary layer, and an adaptive grid scheme redistributing grid points within the computation domain is proposed with a great emphasis on the construction of the grid control function. The function is sensitized to the second invariant of the mean strain tensor, its spatial gradient, and the interaction of pressure gradient and flow deformation. The transition process is assumed to be described with a κ-ε turbulence model. An elliptic solver is employed to integrate governing equations. Numerical results show that the proposed adaptive grid scheme is very effective in obtaining grid independent numerical solution with a very low grid number. It is expected that present scheme is helpful in predicting actual flow within a turbine to improve computation efficiency.

• Journal of Institute of Control, Robotics and Systems
• /
• v.22 no.11
• /
• pp.932-936
• /
• 2016

In this paper, we present a method for the inspection of diamond wheels. In total, six items, including height, radius, and angle, need to be checked during the manufacturing of a diamond wheel. Automatic inspection through image processing is presented in this paper. First, a contour corresponding to the boundary of the diamond wheel is extracted from an image. Next, control points are selected by processing the contour. Seven control points are detected and used for the computation of the required item. Detailed procedures for the computation of the height, radius, and angle using control points are presented in this paper. Experimental results show the feasibility of the presented method.