• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2005.04a
• /
• pp.90-93
• /
• 2005

A thermal model of carbon spun yarn and its composite is presented. Based on voxelization method, the unit cells of spun carbon yam and its composite are divided into a number of volume elements and the local material properties have been given to each element. By using Finite Difference Method, temperature distribution in the unit cell can be obtained. Effective thermal conductivity of unit cell is calculated using the temperature distribution and thermal conductivities of local elements.

• Journal of the Korean Wood Science and Technology
• /
• v.22 no.2
• /
• pp.61-70
• /
• 1994

This study was executed to find the applicability of press drying of tree disk by investigating the shrinkage and drying defect and to form appropriate model by comparing the actual moisture content(MC) and internal temperature in respect of drying time with calculated values based computer simulation to which was applied finite difference method. In press drying disk, heating period, constant drying rate period maintained plateau temperature at 100$^{\circ}C$ and falling drying rate period were significantly distinguished. Actual MC and internal temperature were analogous to those calculated at comparing points. Heat transfer model formed by Fourier's law using specific heat of moist wood and conduction coefficient considering fractional volume of each element of wood cell wall, bound water, free water and air showed applicability as basic data to developing heat expansion, shrinkage and drying stress during press drying. Also mass transfer model formed by Fick's diffusion law using water vapor diffusion coefficient showed applicability. Longitudinal shrinkage was developed by pressure of hot press and tangential shrinkage was restrained by hygrothermal recovery. The heart check, surface check and ring failure were occurred differently in species, but V-shaped crack didn't develop.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.55-63
• /
• 1996

Numerical computations are carried out to analyze the characteristics of flow fields around Air Supported Ships. The computations are performed in a rectangular grid system based on MAC(Marker And Cell) method. The governing equations are represented in finite difference forms by forward differencing in time and centered differencing in space except for its convection terms. For the certification of this numerical analysis method, the computations of flow fields around a Catamaran, an ACV(Air Cushion Vehicle) modeled with pressure distribution on free surface and two SES(Surface Effect Ship)'s are carried out, The results of the present computations are compared with the previously presented computational and experimental results in the same condition.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.2
• /
• pp.77-94
• /
• 1991

A finite difference method based on MAC method is used to simulate free-surface waves around a ship. Euler equations and continuity equation are differentiated using the forward time and central space, and solved by time marching scheme. By the employment of variable mesh system in horizontal and vertical direction, the numerical accuracy of wave simulation results is grossly improved. To verify the improvement of numerical accuracy, some numerical simulations are accomplished for Wigley, Series 60($C_{b}$=0.6) and a bulk carrier model. The computational results are compared to the various experimental data and show good agreements.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.4
• /
• pp.48-55
• /
• 2000

A basic research of air bubble locus around the ship hull is performed to know the bubble's behavior which is generated around the hull surface. In this paper. bubble's behavior around the ship hull is calculated and experimented with the variation of bubble size and the location of bubble generation. For the use of the equation for bubble locus, Kawakita's equation which include the effect of bubble buoyancy is adopted. For the numerical simulation. a finite difference method based on the MAC method is used. Also experiments are performed about series 60 hull form at CWC. The results of calculation are compared with those of experiment.

• Nuclear Engineering and Technology
• /
• v.21 no.1
• /
• pp.40-47
• /
• 1989

For the effective design of a beam forming structure of the negative-ion-based neutral beam injector, a computer program based on a particle simulation model is developed for the calculation of charged particle motions in the electrostatic fields. The motions of negative ions inside the acceleration tube of a multiple-aperture triode are computed at finite time steps. The electrostatic potentials are obtained from the Poisson's equation by the finite difference method. The successive overrelaxation method is used to solve the matrix equation. The particle and force weighting methods are used on a cloud-in-cell model. The optimum design of the beam forming structure has been studied by using this computer code for the various conditions of elctrodes. The effects of the acceleration-deceleration gap distance, the thickness of the deceleration electrode and the shape of the acceleration electrode on beam trajectories are exmined to find the minimum beam divergence. Some numerical illustrations are presented for the particle movements at finite time steps in the beam forming tubes. It is found in this particle simulation modelling that the shape of the acceleration electrode is the most significant factor of beam divergence.

• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.39-46
• /
• 2011

A numerical study of a turbulent natural convection in an enclosure with the lattice Boltzmann method (LBM) is presented. The primary emphasis of the present study is placed on investigation of accuracy and numerical stability of the LBM for the turbulent natural convection flow. A HYBRID method in which the thermal equation is solved by the conventional Reynolds averaged Navier-Stokes equation method while the conservation of mass and momentum equations are resolved by the LBM is employed in the present study. The elliptic-relaxation model is employed for the turbulence model and the turbulent heat fluxes are treated by the algebraic flux model. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with the deferred correction way to ensure accuracy and stability of solutions. The present LBM is applied to the prediction of a turbulent natural convection in a rectangular cavity and the computed results are compared with the experimental data commonly used for the validation of turbulence models and those by the conventional finite-volume method. It is shown that the LBM with the present HYBRID thermal model predicts the mean velocity components and turbulent quantities which are as good as those by the conventional finite-volume method. It is also found that the accuracy and stability of the solution is significantly affected by the treatment of the convection term, especially near the wall.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.2
• /
• pp.93-114
• /
• 2019

In this article, we introduce a finite difference method for solving the Navier-Stokes equations in rectangular domains. The method is proved to be energy stable and shown to be second-order accurate in several benchmark problems. Due to the guaranteed stability and the second order accuracy, the method can be a reliable tool in real-time simulations and physics-based animations with very dynamic fluid motion. We first discuss a simple convection equation, on which many standard explicit methods fail to be energy stable. Our method is an implicit Runge-Kutta method that preserves the energy for inviscid fluid and does not increase the energy for viscous fluid. Integration-by-parts in space is essential to achieve the energy stability, and we could achieve the integration-by-parts in discrete level by using the Marker-And-Cell configuration and central finite differences. The method, which is implicit and second-order accurate, extends our previous method [1] that was explicit and first-order accurate. It satisfies the energy stability and assumes rectangular domains. We acknowledge that the assumption on domains is restrictive, but the method is one of the few methods that are fully stable and second-order accurate.

• Journal of Ship and Ocean Technology
• /
• v.10 no.4
• /
• pp.1-11
• /
• 2006

The finite volume based multi-block RANS code, WAVIS developed at MOERI is applied to the numerical self-propulsion test. WAVIS uses the cell-centered finite volume method for discretization of the governing equations. The realizable $k-{\epsilon}$ turbulence model with a wall function is employed for the turbulence closure. The free surface is captured with the two-phase level set method and body forces are used to model the effects of a propeller without resolving the detail blade flow. The propeller forces are obtained using an unsteady lifting surface method based on potential flow theory. The numerical procedure followed the self-propulsion model experiment based on the 1978 ITTC performance prediction method. The self-propulsion point is obtained iteratively through balancing the propeller thrust, the ship hull resistance and towing force that is correction for Reynolds number difference between the model and full scale. The unsteady lifting surface code is also iterated until the propeller induced velocity is converged in order to obtain the propeller force. The self-propulsion characteristics such as thrust deduction, wake fraction, propeller efficiency, and hull efficiency are compared with the experimental data of the practical container ship. The present paper shows that hybrid RANS and potential flow based numerical method is promising to predict the self-propulsion parameters of practical ships as a useful tool for the hull form and propeller design.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.13 no.4
• /
• pp.267-280
• /
• 2009

In this paper we investigate a simple two-level additive Schwarz preconditioner for the P1 symmetric interior penalty Galerkin method of the Poisson equation on rectangular meshes. The construction is based on the decomposition of the global space of piecewise linear polynomials into the sum of local subspaces, each of which corresponds to an element of the underlying mesh, and the global coarse subspace consisting of piecewise constants. This preconditioner is a direct combination of the block Jacobi iteration and the cell-centered finite difference method, and thus very easy to implement. Explicit upper and lower bounds for the maximum and minimum eigenvalues of the preconditioned matrix system are derived and confirmed by some numerical experiments.