• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.31-35
• /
• 2006

The delta-formulation of the Navier-Stokes equations has been popularly used in the aerodynamics area. Implicit algorithm can be easily implemented in that by using Taylor series expansion. This formulation is extended for an unsteady analysis by using a dual-time integration. In the meanwhile, the incompressible flows with heat transfers which occur in the area of thermo-hydraulics have been solved by a segregated algorithm such as the SIMPLE method, where each equation is discretised by using an under-relaxed deferred correction method and solved sequentially. In this study, the dual-time delta formulation is implemented in the segregated Navier-Stokes solver which is based on the collocated cell-centerd scheme with un unstructured mesh FVM. The pressure correction equation is derived by the SIMPLE method. From this study, it was found that the Euler dual-time method in the delta formulation can be combined with the SIMPLE method.

• Journal of the Society of Naval Architects of Korea
• /
• v.32 no.4
• /
• pp.38-47
• /
• 1995

The three-dimensional multi-layer computer simulation technology for tidal current developed in the previous study is updated to a new version. many improvements are achieved by following changes : (1) No-reflection condition is adopted instead of no-gradient condition as an open boundary condition. (2) Time marching algorithm is changed so that velocity and pressure(surface movement) might be salved in turn at different time step (3) Convection term in equation of motion is estimated by upwind differencing scheme instead of central differencing. The stability is improved considerably and the steady state is achieved within 2 tidal periods which is about 3 times shorter than that of the old version. Moreover, fluctuations in time disappeared by introducing the new time marching technique. An application to the real near shore area(near Inchon harbor) is performed by the new version. Simulated results are compared with those by the simulation total developed in the University of Tokyo. Validity and effectiveness of the two simulation technologies are chocked through the comparative research works.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.137-143
• /
• 1999

The unsteady transonic viscous flow has been analyzed over an oscillatory airfoil. The CSCM(Conservative Supra Characteristic Method) upwind flux difference splitting method and the iterative time marching scheme having first order accuracy in time and second to third order accuracy in space was applied on dynamic meshes. A steady flow field of Mach number 0.7 has been calculated for the verification of unsteady algorithm. The time-accurate unsteady calculations have been done for NACA 0012 airfoil oscillating around quarter chord about freestream Mach number 0.6 on dynamic meshes. The results have been compared with the AGARD Case 3 experimental data. The periodic characteristics have been compared with the experimental results.

• Fire Science and Engineering
• /
• v.11 no.4
• /
• pp.3-14
• /
• 1997

The smoke filling process for the atrium space containing a fire source is simulated using two types of deterministic fire model : Zone model and Field model. The zone model used is the CFAST(version 1.6) model developed at the Building and Fire Research Laboratories, NIST in the USA. The field model is a self-developed frie field model based on Computational Fluid Dynamic (CFD) theories. This article is focused on finding out the smoke movement and temperature distribution in atrium space which is cubic in shape. For solving the liked set of velocity and pressure equation, the PISO algorithm, which strengthened the velocity-pressure coupling, was used. Since PISO algorithm is a time-marching procedure, computing time si very fast. A computational procedure for predicting velocity and temperature distribution in fire-induced flow is based on the solution, in finite volume method and non-staggered grid system, of 3-dimensional equations for the conservation of mass, momentum, energy, species and so forth. The fire model i.e Zone model and Field model predicted similar results for clear heights and the smoke layer temperature.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.60-68
• /
• 1991

The updated Lagrangian Finite Element Method is introduced to analyse rigid body-fluid impact problem which is characterized by incompressible Navier-Stokes equations and impact-contact conditions between free surface and rigid body. For the convenience of numerical computation, velocity fields are splinted into vicous and pressure parts, and then the governing equations and boundary conditions are decomposed in accordance with the decomposition. However, Viscous stresses acting an the solid boundaries are neglected on the assumption that very small velocity gradients may occur during extremely small time interval of the impact. Four coded quadrilateral elements are used to discretize the space domain and the fully explicit time-marching algorithm is employed with a reasonably small time step. At the beginning of each time step, contact velocity of the rigid body is computed from the momentum balance between the body and the fluid. The velocity field is then computed to satisfy the discretized equations of motions and incompressibility and contact constraints as well as an exact free surface boundary condition. At the end of each time step, the fluid domain is updated from the velocity field. In the present time stepping numerical analysis, behaviour of the free surface near the body can be observed without any difficulty which is very important in the water impact problem. The applicability of the algorithm is illustrated by a wedge type falling body problem. The numerical solutions for time-varying pressure distributions and impact loadings acting ion the surface are obtained.

• Journal of Ocean Engineering and Technology
• /
• v.17 no.6
• /
• pp.7-15
• /
• 2003

A Digital Wave Tank simulation technique, based on a finite-difference method and a modified marker-and-cell (MAC) algorithm, is applied in order to investigate the characteristics of nonlinear Tsunami propagations and their interactions with a 2D sloping beach, Ohkushiri Island, and to predict maximum wove run-up around the island. The Navier-Stokes (NS) and continuity equation are governed in the computational domain, and the boundary values are updated at each time step, by a finite-difference time-marching scheme in the frame of the rectangular coordinate system. The fully nonlinear, kinematic, free-surface condition is satisfied by the modified marker-density function technique. The near shore Tsunami is assumed to be a solitary wave, and is generated from the numerical wave-maker in the developed Digital Wave Tank. The simulation results are compared with the experiments and other numerical methods, based on the shallow-water wave theory.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2003.05a
• /
• pp.231-239
• /
• 2003

A Digital Wave Tank simulation technique based on a finite-difference method and a modified marker-and-cell (MAC) algorithm is applied to investigate the characteristics of nonlinear Tsunami propagations and their interactions with a 2D sloping beach and Ohkushiri island, and to predict maximum wave run-up around the island. The Navier-Stokes (NS) and continuity equation are governed in the computational domain and the boundary values updated at each time step by a finite-difference time-marching scheme in the frame of rectangular coordinate system. The fully nonlinear kinematic free-surface condition is satisfied by the modified marker-density function technique. The Nearshore Tsunami is assumed to be a solitary wave and generated from the numerical wavemaker in the developed Digital Wave Tank. The simulation results are compared with the experiments and other numerical methods based on the shallow-water wave theory.

• Journal of Energy Engineering
• /
• v.4 no.1
• /
• pp.126-139
• /
• 1995

저속의 압축성 유동장에서의 원형 실린더 주위의 유동 및 열전달특성을 해석하였다. 비압축성 유동장에서의 실린더 주위의 유동 및 열전달현상에 대하여는 실험과 수치해석을 포함한 광범위한 연구가 진행되어 왔으며 매우 잘 알려져 있다. 실린더 벽면과 주위 유동장의 온도차가 큰 경우, 밀도의 변화가 커지므로 유동장은 압축성 유체가 되나 지배 방정식의 복잡함과 적절한 수치해석 방법의 부족으로 실린더 주위의 유동장을 압축성유체로 해석한 경우는 매우 드물다. 현재 압축성유동 해석에 널리 사용되는 time marching algorithm은 저속의 유동장 해석시 지배방정식에 나타나는 eigenvalue들의 괴리에 의하여 수렴속도가 현저히 떨어지게 된다. 본 연구에서는 이와 같은 난점을 극복할 수 있는 time-derivative preconditioning 방법을 사용하여 온도차가 큰 유동장에서의 강제 및 혼합대류에 대한 계산을 수행하였고 이들의 열전달특성을 비교하였다. 강제대류의 경우 실린더 벽면 온도의 증가에 따른 밀도 감소의 영향은 유동장의 Re수를 감소시켜 확산의 영향을 증가시키면서, 혼합대류의 경우 부력의 영향은 가열되는 유동장의 범위와 재순환 영역을 강제대류에 비해 현저히 감소시킨다.

• The Journal of the Korea institute of electronic communication sciences
• /
• v.8 no.1
• /
• pp.191-197
• /
• 2013

Conventional CT and MRI scans produce cross-section slices of body that are viewed sequentially by radiologists who must imagine or extrapolate from these views what the 3 dimensional anatomy should be. By using sophisticated algorithm and high performance computing, these cross-sections may be rendered as direct 3D representations of human anatomy. The 2D medical image analysis forced to use time-consuming, subjective, error-prone manual techniques, such as slice tracing and region painting, for extracting regions of interest. To overcome the drawbacks of 2D medical image analysis, combining with medical image processing, 3D visualization is essential for extracting anatomical structures and making measurements. We used the gray-level thresholding, region growing, contour following, deformable model to segment human organ and used the feature vectors from texture analysis to detect harmful cancer. We used the perspective projection and marching cube algorithm to render the surface from volumetric MR and CT image data. The 3D visualization of human anatomy and segmented human organ provides valuable benefits for radiation treatment planning, surgical planning, surgery simulation, image guided surgery and interventional imaging applications.

• Coupled systems mechanics
• /
• v.1 no.4
• /
• pp.361-379
• /
• 2012

This paper presents a fully coupled three-dimensional solver for the analysis of interaction between pulsatile flow and large deformation structure. A partitioned time marching algorithm is employed for the solution of the time dependent coupled discretised problem, enabling the use of highly developed, robust and well-tested solvers for each field. Conservative transfer of information at the fluid-structure interface is combined with an effective multi-predict-correct iterative scheme to enable implicit coupling of the interacting fields at each time increment. The three-dimensional unsteady incompressible fluid is solved using a powerful implicit time stepping technique and an ALE formulation for moving boundaries with second-order time accurate is used. A full spectrum of total variational diminishing (TVD) schemes in unstructured grids is allowed implementation for the advection terms and finite element shape functions are used to evaluate the solution and its variation within mesh elements. A finite element dynamic analysis of the highly deformable structure is carried out with a numerical strategy combining the implicit Newmark time integration algorithm with a Newton-Raphson second-order optimisation method. The proposed model is used to predict the wave flow fields of a particular flow-induced vibrational phenomenon, and comparison of the numerical results with available experimental data validates the methodology and assesses its accuracy. Another test case about three-dimensional biomedical model with pulsatile inflow is presented to benchmark the algorithm and to demonstrate the potential applications of this method.