• Journal of the Korean Mathematical Society
• /
• v.59 no.5
• /
• pp.891-909
• /
• 2022

The main objective of this paper is to develop a concrete inverse formula of the system induced by the fourth-order finite difference method for two-point boundary value problems with Robin boundary conditions. This inverse formula facilitates to make a fast algorithm for solving the problems. Our numerical results show the efficiency and accuracy of the proposed method, which is implemented by the Thomas algorithm.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.34 no.2
• /
• pp.174-184
• /
• 1998

In order to accumulate fundamental. data for control of fishes’ behavior at the real fishing ground, the fitness of the numerical model for describing the behavior of fishes was examined by the marine fish. Mullet, Mugil cephalus were used as experimental fishes. The numerical model of fishes’ behavior presented in our earlier paper was modified on the vertical movement of fish school. For the comparision of parameters of the modified numerical model between mullet and rainbow trout, the estimated values of parameters were identified with dimension. The fitness of the modified numerical model was examined by the comparision between experiment and simulation on the several indexes represented by fishes’ swimming characteristics. The obtained result are summarized a follows : 1. The non-dimensional parameter a’ of propulsive force and kb’ of interactive force by the experiment without model net showed a similarity, but the non-dimensional parameter k sub(c’) of schooling force for rainbow trout was lager than one for mullet and the non-dimensional parameter k sub(w’) of repulsive force for mullet was lager than one for rainbow trout. 2. The non-dimensional parameter a’ and k sub(b’) for rainbow trout by the experiment with model net were a little lager than ones for mullet, but non-dimensional parameter k sub(c’) and k sub(w’) for mullet were lager than ones for rainbow trout. 3. The non-dimensional parameter k sub(c’) and k sub(b’) showed the largest and the smallest value among the non-dimensional parameters for rainbow trout and mullet, respectively. 4. The fitness of the modified numerical model was confirmed by means of the compulsion between experiment and simulation on the swimming trajectory of fishes, the mean distance of individual from wall, the mean swimming speed, the mean swimming depth and the mean distance between the nearest individuals. Especially, the similarity of mean swimming depth was improved by using the modified numerical model.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2261-2266
• /
• 2007

In this study, the numerical analysis to estimate condensation heat and mass transfer of the condenser was carried out using the PMA (porous medium approach). In the PMA, the details of tube bundle in the condenser are replaced by the porous medium, and the flow resistance term is added in the momentum equation. In this regard, the PMA is quite helpful for the study of tube bundle in the large condenser. The pressure loss through tube bundle can be compensated by viscous and inertial momentum sink terms, which was validated numerically. Value of the pressure drop was compared to that of Butterworth correlation. Three dimensional analysis of condensation for McAllister condenser with the PMA was conducted using Fluent 6.2 and UDFs (use-defined functions). The result of condensation rate was analogous to previous results (experimental and numerical data).

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2050-2054
• /
• 2007

This study numerically and experimentally investigated on thermal strain analysis of aluminum alloy casting mold using metal foundry. To predict the numerical result of thermal strain in Al alloy casting mold during the cooling process, it is performed the investigation of temperature distribution, stress and displacement based on the physical properties of Al alloy. In results of this study, Al alloy casting mold represented rapidly cooling graph during initial 20minutes after beginning cooling process, therefore value of stress and displacement is rapidly changed during initial 20minutes after beginning cooling process. In addition to, temperature distribution obtained by experiment confirmed corresponding pattern then compared numerical analysis with experiment. These results are distribute to make the effective and the high precision casting mold.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2684-2689
• /
• 2007

Non-equilibrium first order extrapolation boundary condition proposed by Guo et $al.^{(9)}$ proposed has a good application for complex geometries, a second order accuracy and a treatment on non-slip wall boundary condition easily. However it has a lack of the numerical stability from high Reynolds number. Guo et $al.^{(9)}$ substituted the density value of adjacent nodes for the density of boundary nodes. This procedure causes the numerical instability on the boundary. In this paper, we derived a procedure of density extrapolation and compared to previous results.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.2 no.1
• /
• pp.39-44
• /
• 2010

Recently underwater systems moving at high speed such as a super-cavitating torpedo have been studied for their practical advantage of the dramatic drag reduction. In this study we are focusing our attention on super-cavitating flows around axisymmetric cavitators. A numerical method based on inviscid flow is developed and the results for several shapes of the cavitator are presented. First using a potential based boundary element method, we find the shape of the cavtiator yielding a sufficiently large enough cavity to surround the body. Second, numerical predictions of supercavity are validated by comparing, with experimental observations carried out in a high speed cavitation tunnel at Chungnam National University (CNU CT).

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.9
• /
• pp.1148-1156
• /
• 2000

The QUICK scheme for convection terms is modified for unstructured finite volume method by using linear reconstruction technique and validated through the computation of two well defined laminar flows. It uses two upstream grid points and one downstream grid point in approximating the convection terms. The most upstream grid point is generated by considering both the direction of flow and local grid line. Its value is calculated from surrounding grid points by using a linear construction method. Numerical error by the modified QUICK scheme is shown to decrease about 2.5 times faster than first order upwind scheme as grid size decreases. Computations are also carried out to see effects of the skewness and irregularity of grid on numerical solution. All numerical solutions show that the modified QUICK scheme is insensitive to both the skewness and irregularity of grid in terms of the accuracy of solution.

• East Asian mathematical journal
• /
• v.33 no.5
• /
• pp.495-509
• /
• 2017

With the accuracy of the nonlinearity guaranteed, plenty of time and large memory space are needed when we solve the finite element numerical solution of nonlinear partial differential equations. In this paper, we use the Group Element Method (GEM) to deal with the non-linearity of the BBM-Burgers Equation with Conservation form and perform a numerical analysis for two particular initial-boundary value (the Dirichlet boundary conditions and Neumann-Dirichlet boundary conditions) problems with the Finite Element Method (FEM). Some numerical experiments are performed to analyze the error between the exact solution and the FEM solution in MATLAB.

• 한국연소학회:학술대회논문집
• /
• 2014.11a
• /
• pp.297-299
• /
• 2014

Petroleum coke has high heating value and low price. Due to the steadily increasing demand for heavy oil processing, the production volume of petroleum coke tends to be expanded. The high availability and low price of petroleum coke have been strongly considered as candidate fuel for power generation facilities. However the high carbon content, high sulfur content and nitrogen content of petroleum fuel are known to produce relatively large quantity of CO2, high NOx and SO2 emission. In this work, a series of numerical simulations have been carried out in order to investigate the effects of swirl flow intensity on combustion furnace, which is most important operating condition. Results show that the temperature distribution was spatially uniform at about 1600K but high temperature region are located quite differently depending on swirl number. In addition, numerical temperature data was compared with experimental temperature data and its temperature difference shows less than 10%. On the other hand, discrepancy between numerical and experimental emission data were slightly large with necessities of improved emission model.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.3
• /
• pp.624-634
• /
• 1994

A numerical method which combines equal-order velocity-pressure formulation originated from SIMPLE algorithm and streamline upwinding method has been developed. To verify the proposed numerical method, we considered the lid-driven cavity flow and backward facing step flow. The trend of convergence history is stable up to the error criterion beyond which the maximum value of error is oscillatory due4 to the round-off error. In the present study, all results were obtained with the single precision calculation up to the given error criterion and it was found to be sufficient for our purpose. The present results were then compared with existing experimental results using laser doppler velocimetry and numerical results using finite difference method and mixed interpolation finite element method. It has been shown that the present method gives accurate results with less memories and execution time than the coventional finite element method.