• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.194-199
• /
• 2007

IB (immersed boundary) method is one of the prominent tool in computational fluid dynamics for the analysis of flows over complex geometries. The IB technique simplyfies the solution procedure by eliminating the requirement of complex body fitted grids and it is also superior in terms of memory requirement. In this study we have developed numerical code (FOTRAN) for the analysis of 2D flow over a cylinder using IB technique. The code is validated by comparing the wake lengths and separation angles given by Guo et. al. We employed fractional-step procedure for solving the Navier-Stokes equations governing the flow and discrete forcing IB technique for imposing boundary conditions. Also we have developed a 3D code for the backward-facing-step flow and flow over a sphere. The reattachment length in backward-facing-step flow was compared with the one given by Nie and Armaly, which has proven the validity of our code.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.17 no.4
• /
• pp.237-247
• /
• 1997

A geometrical inverse heat conduction problem is solved for the infrared scanning cavity detection by the boundary element method using minimal energy technique. By minimizing the kinetic energy of temperature field, boundary element equations are converted to the quadratic programming problem. A hypothetical inner boundary is defined such that the actual cavity is located interior to the domain. Temperatures at hypothetical inner boundary are determined to meet the constraints of mea- surement error of surface temperature obtained by infrared scanning, and then boundary element analysis is peformed for the position of an unknown boundary (cavity). Cavity detection algorithm is provided, and the effects of minimal energy technique on the inverse solution method are investigated by means of numerical analysis.

• Communications of the Korean Mathematical Society
• /
• v.24 no.4
• /
• pp.605-615
• /
• 2009

In this paper, we develop a reliable algorithm which is called the variation of parameters method for solving sixth-order boundary value problems. The proposed technique is quite efficient and is practically well suited for use in these problems. The suggested iterative scheme finds the solution without any perturbation, discritization, linearization or restrictive assumptions. Moreover, the method is free from the identification of Lagrange multipliers. The fact that the proposed technique solves nonlinear problems without using the Adomian's polynomials can be considered as a clear advantage of this technique over the decomposition method. Several examples are given to verify the reliability and efficiency of the proposed method. Comparisons are made to reconfirm the efficiency and accuracy of the suggested technique.

• Communications of the Korean Mathematical Society
• /
• v.22 no.1
• /
• pp.53-64
• /
• 2007

In this paper, we apply the generalized quasilinearization technique to a forced Duffing equation with three-point mixed nonlinear nonlocal boundary conditions and obtain sequences of upper and lower solutions converging monotonically and quadratically to the unique solution of the problem.

• Journal of Mechanical Science and Technology
• /
• v.14 no.1
• /
• pp.84-92
• /
• 2000

A new bounce back boundary method of the second order in error is proposed for the lattice Boltzmann fluid simulation. This new method can be used for the arbitrarily irregular lattice geometry of a non-slip boundary. The traditional bounce back boundary condition for the lattice Boltzmann simulation is of the first order in error. Since the lattice Boltzmann method is the second order scheme by itself, a boundary technique of the second order has been desired to replace the first order bounce back method. This study shows that, contrary to the common belief that the bounce back boundary condition is unilaterally of the first order, the second order bounce back boundary condition can be realized. This study also shows that there exists a generalized bounce back technique that can be characterized by a single interpolation parameter. The second order bounce back method can be obtained by proper selection of this parameter in accordance with the detailed lattice geometry of the boundary. For an illustrative purpose, the transient Couette and the plane Poiseuille flows are solved by the lattice Boltzmann simulation with various boundary conditions. The results show that the generalized bounce back method yields the second order behavior in the error of the solution, provided that the interpolation parameter is properly selected. Coupled with its intuitive nature and the ease of implementation, the bounce back method can be as good as any second order boundary method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.37 no.1
• /
• pp.41-47
• /
• 2024

This paper utilizes computational asymptotic analysis to compute the boundary layer solution for composite beams and validates the findings through a comparison with ANSYS results. The boundary layer solution, presented as a sum of the interior solution and pure boundary layer effects, necessitates a mathematically rigorous formalization for both interior and boundary layer aspects. Computational asymptotic analysis emerges as a robust technique for addressing such problems. However, the challenge lies in connecting the boundary layer and interior solutions. In this study, we systematically separate the principles of virtual work and the principles of Saint-Venant to tackle internal and boundary layer issues. The boundary layer solution is articulated by calculating the Papkovich-Fadle eigenfunctions, representing them as linear combinations of real and imaginary vectors. To address warping functions in the interior solutions, we employed a least squares method. The computed solutions exhibit excellent agreement with 2D finite element analysis results, both quantitatively and qualitatively. This validates the effectiveness and accuracy of the proposed approach in capturing the behavior of composite beams.

• Journal of applied mathematics & informatics
• /
• v.31 no.1_2
• /
• pp.55-63
• /
• 2013

In this paper, a general technique is proposed for solving a system of fourth-order boundary value problems. The solution is given in the form of series and its approximate solution is obtained by truncating the series. Advantages of the method are that the representation of exact solution is obtained in a new reproducing kernel Hilbert space and accuracy of numerical computation is higher. Numerical results show that the method employed in the paper is valid. Numerical evidence is presented to show the applicability and superiority of the new method.

• 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.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.2
• /
• pp.304-311
• /
• 2001

The present work examines the influence of surface coating on the temperature and the thermo-mechanical stress field produced by friction due to sliding contact. A two-dimensional transient model of a layered medium submitted to a moving heat flux is prsented. A solution technique based on the boundary element method employing the multiregion technique is utilized. Results are presented showing the influence of coating thickness, thermal properties, Peclet number, and mechanical properties. It has been shown that the mechanical properties and thickness of coating have a significant influence on the stress field, even for low temperature increase. The effects of the ratios of shear modulus become more important for low temperature increase than the effects of the ratios of other mechanical properties.

• Journal of applied mathematics & informatics
• /
• v.31 no.1_2
• /
• pp.221-228
• /
• 2013

The method of upper and lower solutions and the generalized quasilinearization technique is developed for the existence and approximation of solutions to boundary value problems for higher order fractional differential equations of the type $^c\mathcal{D}^qu(t)+f(t,u(t))=0$, $t{\in}(0,1),q{\in}(n-1,n],n{\geq}2$ $u^{\prime}(0)=0,u^{\prime\prime}(0)=0,{\ldots},u^{n-1}(0)=0,u(1)={\xi}u({\eta})$, where ${\xi},{\eta}{\in}(0,1)$, the nonlinear function f is assumed to be continuous and $^c\mathcal{D}^q$ is the fractional derivative in the sense of Caputo. Existence of solution is established via the upper and lower solutions method and approximation of solutions uses the generalized quasilinearization technique.