• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.4
• /
• pp.1000-1013
• /
• 2014

This paper presents some of the efforts by the authors towards numerical prediction of springing of ships. A time-domain Higher Order Boundary Element Method (HOBEM) based on cubic shape function is first presented to solve a complete second-order problem in terms of wave steepness and ship motions in a consistent manner. In order to avoid high order derivatives on the body surfaces, e.g. mj-terms, a new formulation of the Boundary Value Problem in a body-fixed coordinate system has been proposed instead of traditional formulation in inertial coordinate system. The local steady flow effects on the unsteady waves are taken into account. Double-body flow is used as the basis flow which is an appropriate approximation for ships with moderate forward speed. This numerical model was used to estimate the complete second order wave excitation of springing of a displacement ship at constant forward speeds.

• Journal of applied mathematics & informatics
• /
• v.7 no.2
• /
• pp.453-465
• /
• 2000

In this paper, a new method for finding the approximate solution of a second order nonlinear partial differential equation is introduced. In this method the problem is transformed to an equivalent optimization problem. them , by considering it as a distributed parameter control system the theory of measure is used for obtaining the approximate solution of the original problem.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.9
• /
• pp.1423-1431
• /
• 2001

Based on the exponential intervening variable, a new two-point approximation method is presented. This introduces the shifting level into each exponential intervening variable to avoid the lack of def inition of the conventional exponential intervening variables due to zero-or negative-valued design variables. Then a new quadratic approximation whose Hessian matrix has only diagonal elements of different values is proposed in terms of these intervening variables. These diagonal elements are determined in a closed form that corrects the typical error in the approximate gradient of the TANA series due to the lack of definition of exponential type intervening variables and their incomplete second-order terms. Also, a correction coefficient is multiplied to the pre-determined quadratic term to match the value of approximate function with that of the previous point. Finally, in order to show the numerical performance of the proposed method, a sequential approximate optimizer is developed and applied to solve six typical design problems. These optimization results are compared with those of TANA-3. These comparisons show that the proposed method gives more efficient and reliable results than TANA-3.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.9 no.2
• /
• pp.59-65
• /
• 2000

In this paper the comparison of the first order approximation schemes such as SLP(sequential linear programming) CONLIN(convex linearization) MMA(method of moving asymptotes) and the second order approximation scheme SQP(sequential quadratic programming) was accomplished for optimization of nonlinear structures. It was found that MMA and SQP are the most efficient methods for optimization. But the number of function call of SQP is much more than that of MMA. Therefore when it is considered with the expense of computation MMA is more efficient than SQP. In order to examine the efficiency of MMA for complex optimization problem it was applied to the helicopter tail boom con-sidering column buckling and local wall buckling constraints. it is concluded that MMA can be a very efficient approxima-tion scheme from simple problems to complex problems.

• The Journal of the Acoustical Society of Korea
• /
• v.33 no.4
• /
• pp.232-237
• /
• 2014

A passive sonar of warship is composed of several directional or omni-directional sensors. In order to model the acoustic signal received into a warship sonar, the wave propagation modeling is usually required from arbitrary noise source to all sensors equipped to the sonar. However, the full calculation for all sensors is time-consuming and the performance of sonar simulator deteriorates. In this study, we suggest an asymptotic method to estimate the sonar signal arrived to sensors adjacent to the reference sensor, where it is assumed that all information of eigenrays is known. This method is developed using Taylor series for the time delay of eigenray and similar to Fraunhofer and Fresnel approximation for sonar aperture. To validate the proposed method, some numerical experiments are performed for the passive sonar. The approximation when the second-order term is kept is vastly superior. In addition, the error criterion for each approximation is provided with a practical example.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.7 no.4
• /
• pp.149-154
• /
• 2008

The concept of robust design was introduced by Dr. G. Taguchi in the late 1940s, and his technique has become commonly known as the Taguchi method or the robust design. In this research, a robust design procedure for microgyroscope is suggested based on the kriging and optimization approaches. The kriging interpolation method is introduced to obtain the surrogate approximation model of true function. Robustness is calculated by the kriging model to reduce real function calculations. For this, objective function is represented by the probability of success, thus facilitating robust optimization. The statistics such as mean and variance are obtained based on the reliable kriging model and the second-order statistical approximation method.

• Journal of applied mathematics & informatics
• /
• v.15 no.1_2
• /
• pp.173-184
• /
• 2004

In this paper, we find the approximate solution of a second order nonlinear partial differential equation on a simple connected region in $R^2$. We transfer this problem to a new problem of second order nonlinear partial differential equation on a rectangle. Then, we transformed the later one to an equivalent optimization problem. Then we consider the optimization problem as a distributed parameter system with artificial controls. Finally, by using the theory of measure, we obtain the approximate solution of the original problem. In this paper also the global error in $L_1$ is controlled.

• Journal of applied mathematics & informatics
• /
• v.35 no.3_4
• /
• pp.277-302
• /
• 2017

In this paper, we have proposed a numerical method for Singularly Perturbed Boundary Value Problems (SPBVPs) of reaction-diffusion type of third order Ordinary Differential Equations (ODEs). The SPBVP is reduced into a weakly coupled system of one first order and one second order ODEs, one without the parameter and the other with the parameter ${\varepsilon}$ multiplying the highest derivative subject to suitable initial and boundary conditions, respectively. The numerical method combines boundary value technique, asymptotic expansion approximation, shooting method and finite difference scheme. The weakly coupled system is decoupled by replacing one of the unknowns by its zero-order asymptotic expansion. Finally the present numerical method is applied to the decoupled system. In order to get a numerical solution for the derivative of the solution, the domain is divided into three regions namely two inner regions and one outer region. The Shooting method is applied to two inner regions whereas for the outer region, standard finite difference (FD) scheme is applied. Necessary error estimates are derived for the method. Computational efficiency and accuracy are verified through numerical examples. The method is easy to implement and suitable for parallel computing. The main advantage of this method is that due to decoupling the system, the computation time is very much reduced.

• Structural Engineering and Mechanics
• /
• v.34 no.6
• /
• pp.779-799
• /
• 2010

In order to consider high-order effects on the actual limit state function, a new response surface method is proposed for structural reliability analysis by the use of high-order approximation concept in this study. Hermite polynomials are used to determine the highest orders of input random variables, and the sampling points for the determination of highest orders are located on Gaussian points of Gauss-Hermite integration. The cross terms between two random variables, only in case that their corresponding percent contributions to the total variation of limit state function are significant, will be added to the response surface function to improve the approximation accuracy. As a result, significant reduction in computational cost is achieved with this strategy. Due to the addition of cross terms, the additional sampling points, laid on two-dimensional Gaussian points off axis on the plane of two significant variables, are required to determine the coefficients of the approximated limit state function. All available sampling points are employed to construct the final response surface function. Then, Monte Carlo Simulation is carried out on the final approximation response surface function to estimate the failure probability. Due to the use of high order polynomial, the proposed method is more accurate than the traditional second-order or linear response surface method. It also provides much more efficient solutions than the available high-order response surface method with less loss in accuracy. The efficiency and the accuracy of the proposed method compared with those of various response surface methods available are illustrated by five numerical examples.

• Journal of the Chungcheong Mathematical Society
• /
• v.20 no.4
• /
• pp.419-432
• /
• 2007

Double-diffusive convection inside a triangular porous cavity is studied numerically. Galerkin finite element method is adopted to derive the discrete form of the governing differential equations. The first-order backward Euler scheme is used for temporal discretization with the second-order Adams-Bashforth scheme for the convection terms in the energy and species conservation equations. The Boussinesq-Oberbeck approximation is used to calculate the density dependence on the temperature and concentration fields. A parametric study is performed with the Lewis number, the Rayleigh number, the buoyancy ratio, and the shape of the triangle. The effect of gravity orientation is considered also. Results obtained include the flow, temperature, and concentration fields. The differences induced by varying physical parameters are analyzed and discussed. It is found that the heat transfer rate is sensitive to the shape of the triangles. For the given geometries, buoyancy ratio and Rayleigh numbers are the dominating parameters controlling the heat transfer.