• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.3B
• /
• pp.293-303
• /
• 2011

An efficient diagonalized approximate factorization algorithm (DAF) is developed for the solution of three-dimensional incompressible viscous flows. The pressure-based, artificial compressibility (AC) method is used for calculating steady incompressible Navier-Stokes equations. The AC form of the governing equations is discretized in space using a second-order-accurate finite volume method. The present DAF method is applied to derive a second-order accurate splitting of the discrete system of equations. The primary objective of this study is to investigate the computational efficiency of the present DAF method. The solutions of the DAF method are evaluated relative to those of well-known four-stage Runge-Kutta (RK4) method for fully developed and developing laminar flows in curved square ducts and a laminar flow in a cavity. While converged solutions obtained by DAF and RK4 methods on the same computational meshes are essentially identical because of employing the same discrete schemes in space, both algorithms shows significant discrepancy in the computing efficiency. The results reveal that the DAF method requires substantially at least two times less computational time than RK4 to solve all applied flow fields. The increase in computational efficiency of the DAF methods is achieved with no increase in computational resources and coding complexity.

• Journal of the Korea Computer Graphics Society
• /
• v.9 no.2
• /
• pp.11-18
• /
• 2003

For extracting geometric features in 3D meshes according to user-steering with effective interactions. this paper generalizes the 2D algorithms of snapping and wrapping that. respectively. moves a cursor to a nearby feature and constructs feature boundaries. First. we define approximate curvatures and move cost functions that are the numerical values measuring the geometric characteristics of the meshes, By exploiting the measuring values. the algorithms of geometric snapping and geometric wrapping are developed and implemented. We also visualize the results from applying the algorithms to extracting geometric features of general 3D mesh models such as a face model and a tooth model.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.22 no.4
• /
• pp.411-421
• /
• 2002

The accurate analysis of ultrasonic wave propagation and scattering plays an important role in many aspects of nondestructive evaluation. A numerical analysis makes it possible to perform parametric studies, and in this way the probability of detection and reliability of test results can be improved. In this study, a finite element method was developed for the analysis of ultrasonic fields, the accuracy of results was checked by solving several representative problems. The size of element and the integral time step, which are the critical components for the convergence of numerical results, were determined in a commercial finite element code. Several propagation and scattering problems in 2-D isotropic and anisotropic materials were solved and their results were compared with known analytical or experimental results.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.4
• /
• pp.213-221
• /
• 2013

This paper introduces a mesh continuation scheme for a one-dimensional inverse medium problem to reconstruct the spatial distribution of elastic wave velocities in heterogeneous semi-infinite solid domains. To formulate the inverse problem, perfectly-matched-layers(PMLs) are introduced as wave-absorbing boundaries that surround the finite computational domain truncated from the originally semi-infinite extent. To tackle the inverse problem in the PML-truncated domain, a partial-differential-equations(PDE)-constrained optimization approach is utilized, where a least-squares misfit between calculated and measured surface responses is minimized under the constraint of PML-endowed wave equations. The optimization problem iteratively solves for the unknown wave velocities with their updates calculated by Fletcher-Reeves conjugate gradient algorithms. The optimization is performed using a mesh continuation scheme through which the wave velocity profile is reconstructed in successively denser mesh conditions. Numerical results showed the robust performance of the mesh continuation scheme in reconstructing target wave velocity profile in a layered heterogeneous solid domain.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.4_1
• /
• pp.67-76
• /
• 1992

A hierarchical formulation based on p-version of the finite element method for linear elastic axisymmetric stress analysis is presented. This is accomplished by introducing additional nodal variables in the element displacement approximation on the basis of integrals of Legendre polynomials. Since the displacement approximation is hierarchical, the resulting element stiffness matrix and equivalent nodal load vectors are hierarchical also. The merits of the propoosed element are as follow: i) improved conditioning, ii) ease of joining finite elements of different polynomial order, and iii) utilizing previous solutions and computation when attempting a refinement. Numerical examples are presented to demonstrate the accuracy, efficiency, modeling convenience, robustness and overall superiority of the present formulation. The results obtained from the present formulation are also compared with those available in the literature as well as with the analytical solutions.

• Solar Energy
• /
• v.12 no.3
• /
• pp.47-59
• /
• 1992

In the concentric tube type recuperator, which is the most typical type of radiation recuperator, installed on the exhaust-secion of industrial furnace, air flows between the adiabatic outer tube and the inner tube in which exhausted gas flows with high temperature. The waste heat of the exhausted gas is transferred to the inner tube, and transferred from the inner tube to the flowing air. The heat transfer by radiation In the concentric tube type recuperator is modeled using spherical harmonics approximation, namely, P-N method and numerically analyzed considering the effect of dynamic flow field. The results are compared with the existing empirical data. In addition, a theoretical method is presented for the analysis of the heat transfer characteristics of a recuperator with a reradiator installed in the inner tube, which causes re-radiant in the inner tube, and the characteristics of the recuperator is analyszed and defined.

• 한국HCI학회:학술대회논문집
• /
• 2008.02a
• /
• pp.813-816
• /
• 2008

In this research we present a method to reconstruct the casting flow simulation result as a 3D model and visualize it on a VR display. First, numerical analysis of heat flow is performed using an existing commercial CAE simulation software. In this process the shape of the original design model is approximated to a regular rectangular grid. The filling ratio and temperature of each voxel are recorded iteratively by predefined number of steps starting from pouring the melted metal into a mold until it is entirely filled. Next we reconstruct the casting by voxels using the simulation result as an input. The color of voxel is determined by mapping the colors to temperature and filling ratio at each step as the flow proceeds. The reconstructed model is visualized on the Projection Table which is one of horizontal-type VR display. It provides active stereoscopic images.

• Fire Science and Engineering
• /
• v.15 no.1
• /
• pp.7-15
• /
• 2001

This paper describes the smoke filling process of a fire field model based on a self-developed SMEP (Smoke Movement Estimating Program) codo to the simulation of fire induced flows in the atrium space (SIVANS atrium at Japan) containing smoke radiation effect. The SMEP using PISO algorithm solves conservation equations for mass, momentum, energy and species, together with those for the modified k-$\varepsilon$ turbulence model with buoyancy term. Also it solves the radiation equation using the discrete ordinates method. The result of the calculated smoke temperature containing radiation effect has shown a better prediction than the result calculated by only convection effect in comparison with the experimental data. This seems to come from the radiation effect of $H_2$O and $CO_2$ gas under smoke productions. Thus, the consideration of the radiation effect under smoke in fire should be necessary in order to get more realistic result. Also the numerical results indicated that the smoke layer is developing at a rate of about 0.1 m/s. It would take about 450 seconds after starting the ultra fast fire of 560 kW that the smoke layer move down to 1.5m above the escape level.

• Journal of Korea Water Resources Association
• /
• v.32 no.3
• /
• pp.301-310
• /
• 1999

In this paper a numerical model and two analytical models in the hydraulically connected stream-aquifer system were analyzed to compare the induced infiltration rate curves derived from each model. And we also examined the effects of anisotropy of hydraulic conductivity and the direction of the ambient ground water flow on the quantification of the induced infiltration rate. The induced infiltration rate curve determined by models is very simple and useful for estimating the induced infiltration rate since it contains only four physical variables such as the induced infiltration rate, the pumping rate, the distance between the pumping well and the stream, and the ambient ground water flow rate. Under the conditions tested in this paper the induced infiltration rate curves resulted from the Wilson's analytical model and FEWA numerical model were in good agreement, and the anisotropic ratio of hydraulic conductivity was evaluated as a physical factor which influences the behaviour of the induced infiltration rate curve. The methods and results of the paper might Icad to improve the understanding of the induced infiltration phenomenon and can be applied to the planning and disign of pumping well and the optimal determination of the induced infiltration rate and pumping rate for water quality management of the water supply wells.

• Journal of the Korean Society of Safety
• /
• v.21 no.2 s.74
• /
• pp.143-149
• /
• 2006

This paper deals with the space-time finite element analysis of two-dimensional vibration problem with a single variable. The method of space-time finite elements enables the simpler solution than the usual finite element analysis with discretization in space only. We present a discretization technique in which finite element approximations are used in time and space simultaneously for a relatively large time period. The weighted residual process is used to formulate a finite element method for a space-time domain. A stability problem is described and some investigations for chosen type of rectangular space-time finite elements are carried out. Instability is caused by a too large time step of successive time steps in the traditional time-dependent problems. It has been shown that the numerical stability of time-stepping on the larger time steps is quite good. The unstructured space-time finite element not only overcomes the shortcomings of the stability in the traditional numerical methods, but it is also endowed with the features of an effective computational technique. Some numerical examples have been presented to illustrate the efficiency of the described method.