• Journal of Ocean Engineering and Technology
• /
• v.13 no.3B
• /
• pp.3-13
• /
• 1999

The diffusion characteristics of the high temperature and salinity water discharged in Chinhae Bay under BMP(Barge-Mounted Plants) desalination processes were simulated to access environmental impact. The 3-D baroclinic flow model is formulated by integrating the basic equations with respect to each control volume and by transforming them into a finite difference form using the space-staggered grid system. With a 3-D baroclinic flow model, the tide-induced and density-induced current was computed and confirmed by comparing with observed data. From the results of numerical experiment, it is expected that the maximum diffusion lengths of the high temperature and salinity which increase $0.6^{circ}C$ and 0.2 after discharging are 1 km and 3.5km, respectively. It may be expected that the discharge has an effect on surrounding area of discharge, but not an effect on whole area of Chinhae Bay.

• Structural Engineering and Mechanics
• /
• v.71 no.5
• /
• pp.469-483
• /
• 2019

In this paper, an improved mathematical model is presented for the bending analysis of doubly curved functionally graded material (FGM) sandwich rhombic conoids. The mathematical model includes expansion of Taylor's series up to the third degree in thickness coordinate and normal curvatures in in-plane displacement fields. The condition of zero-transverse shear strain at upper and lower surface of rhombic conoids is implemented in the present model. The newly introduced feature in the present mathematical model is the simultaneous inclusion of normal curvatures in deformation field and twist curvature in strain-displacement equations. This unique introduction permits the new 2D mathematical model to solve problems of moderately thick and deep doubly curved FGM sandwich rhombic conoids. The distinguishing feature of present shell from the other shells is that maximum transverse deflection does not occur at its center. The proposed new mathematical model is implemented in finite element code written in FORTRAN. The obtained numerical results are compared with the results available in the literature. Once validated, the current model was employed to solve numerous bending problems by varying different parameters like volume fraction indices, skew angles, boundary conditions, thickness scheme, and several geometric parameters.

• Journal of computational fluids engineering
• /
• v.17 no.1
• /
• pp.94-100
• /
• 2012

In this paper, the cavitating flows around a hydrofoil have been numerically investigated by using a 2-d multi-phase RANS flow solver based on pseudo-compressibility and a homogeneous mixture model on unstructured meshes. For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras one equation model was employed for the closure of turbulence. A dual-time stepping method and the Gauss-Seidel iteration were used for unsteady time integration. The phase change rate between the liquid and vapor phases was determined by Merkle's cavitation model based on the difference between local and vapor pressure. Steady state calculations were made for the modified NACA66 hydrofoil at several flow conditions. Good agreements were obtained between the present results and the experiment for the pressure coefficient on a hydrofoil surface. Additional calculation was made for cloud cavitation around the hydrofoil. The observation of the vapor structure, such as cavity size and shape, was made, and the flow characteristics around the cavity were analyzed. Good agreements were obtained between the present results and the experiment for the frequency and the Strouhal number of cavity oscillation.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.6 no.2
• /
• pp.12-20
• /
• 1998

Exter ballistics of a typical high-speed projectile is studied through a flow-visualization experiment and an unstructured grid Navier-Srokes computation. Experiment produced a schlieren photograph that adequately shows the characteristic features of this complex flow, namely two kinds of oblique cone shocks and turbulent wake developing into the downstream. A hybrid scheme of finite volume-element method is used to simulate the compressible Reynolds-Averaged Navier-Stok- es solution on unstructured grids. Osher's approximate Riemann solver is used to discretize the cinvection term. Higher-order spatial accuracy is obtained by MUSCL extension and van Albada ty- pe flux limiter is used to stabilize the numerical oscillation near the solution discontinuity. Accurate Gakerkin method is used to discretize the viscous term. Explict fourth-order Runge-Kutta method is used for the time-stepping, which simplifies the application of MUSCL extension. A two-layer k-$\varepsilon$ turbulence model is used to simulate the turbulent wakes accurately. Axisymmetric folw and two-dimensional flow with an angle of attack have been computed. Grid-dependency is also checked by carrying out the computation with doubled meshes. 2-D calculation shows that effect of angle of attack on the flow field is negligible. Axi-symmetric results of the computation agrees well with the flow visualization. Primary oblique shock is represented within 2-3 meshes in numerical results, and the varicose mode of the vortex shedding is clearly captured in the turbulent wake region.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.7 no.1
• /
• pp.7-11
• /
• 2006

Through SCM modeling, we found that the depletion layer in silicon was of a form of a spherical capacitor with the SCM tip biased. Two-dimensional (2D) finite differential method code with a successive over relaxation (SOR) solver has been developed to model the measurements by SCM of a semiconductor wafer that contains an ion-implanted impurity region. Then, we theoretically analyzed the spherical capacitor and determined the total depleted-volume charge Q, capacitance C, and the rate of capacitance change with bias dC/dV. It is very important to observe the depleted carriers' movement in the silicon layer by applying the bias to the tip. So, we calculated the depleted-volume charge, considering different factors such as tip size, oxide thickness, and applied bias (dc+ac), which have an influence on potential and depletion charges.

• Journal of computational fluids engineering
• /
• v.8 no.2
• /
• pp.49-56
• /
• 2003

Three-dimensional endwall flow within a linear cascade passage of high performance turbine blade is simulated with a 3-D Navier-Stokes CFD code (MOSA3D), which is based on body-fitted coordinate system, pressure-correction and finite volume method. The endwall flow characteristics, including the development and generation of horseshoe vortex, passage vortex, etc. are clearly simulated, consistent with the generally known tendency. The effects of both turbulence model and convective differencing scheme on the prediction performance of endwall flow are systematically analyzed in the present paper. The convective scheme is found to have stronger effect than the turbulence model on the prediction performance of endwall flow. The present simulation result also indicates that the suction leg of the horseshoe vortex continues on the suction side until it reaches the trailing edge.

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

A three-dimensional (3D) unstructured hydrodynamic solver for transient two-phase flows has been developed. A two-fluid three-field model was adopted for the two-phase flows. The three fields represent a continuous liquid, an entrained liquid, and a vapour field. The hydrodynamic solver is for the 3D component of a nuclear system code and the component-scale analysis tools for transient two-phase flows. The finite volume method and unstructured grid are adopted, which are useful for the flows in a complicated geometry. The semi-implicit ICE (Implicit Continuous-fluid Eulerian) numerical scheme has been adapted to the unstructured non-staggered grid. This paper presents the numerical method and the preliminary results of the calculations. The results show that the numerical scheme is robust and predicts the phase change and the flow transitions due to boiling and flashing problems well.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.10
• /
• pp.1609-1618
• /
• 1997

Thermo-mechanical behaviors of polymer matrix composite(PMC) with continuous TiNi fiber are studied using theoretical analysis with 1-D analytical model and numerical analysis with 2-D multi-fiber finite element(FE) model. It is found that both compressive stress in matrix and tensile stress in TiNi fiber are the source of strengthening mechanisms and thermo-mechanical coupling. Thermal expansion of continuous TiNi fiber reinforced PMC has been compared with various mechanical behaviors as a function of fiber volume fraction, degree of pre-strain and modulus ratio between TiNi fiber and polymer matrix. Based on the concept of so-called shape memory composite(SMC) with a permanent shape memory effect, the critical modulus ratio is determined to obtain a smart composite with no or minimum thermal deformation. The critical modulus ratio should be a major factor for design and manufacturing of SMC.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.2
• /
• pp.219-235
• /
• 2014

With the development of the technology of underwater moving bodies, the need for developing the knowledge of surface effect interaction of free surface and underwater moving bodies is increased. Hence, the two-phase flow is a subject which is interesting for many researchers all around the world. In this paper, the non-linear free surface deformations which occur during the water-exit of a circular cylinder due to its buoyancy are solved using finite volume discretization based code, and using Volume of Fluid (VOF) scheme for solving two phase flow. Dynamic mesh model is used to simulate dynamic motion of the cylinder. In addition, the effect of cylinder mass in presence of an external force is studied. Moreover, the oblique exit and entry of a circular cylinder with two exit angles is simulated. At last, water-exit of a circular cylinder in six degrees of freedom is simulated in 3D using parallel processing. The simulation errors of present work (using VOF method) for maximum velocity and height of a circular cylinder are less than the corresponding errors of level set method reported by previous researchers. Oblique exit shows interesting results; formation of waves caused by exit of the cylinder, wave motion in horizontal direction and the air trapped between the waves are observable. In 3D simulation the visualization of water motion on the top surface of the cylinder and the free surface breaking on the front and back faces of the 3D cylinder at the exit phase are observed which cannot be seen in 2D simulation. Comparing the results, 3D simulation shows better agreement with experimental data, specially in the maximum height position of the cylinder.

• Proceedings of the KIEE Conference
• /
• 2005.07b
• /
• pp.1306-1308
• /
• 2005

The characteristic of line start permanent magnet(LSPM) motor with overhang structure is analyzed using 3dimensional finite element method(3D FEM). We propose the appropriate overhang length by considering the magnetic flux distribution according to changing overhang length and manufacturing process. 3D FEM analysis is inevitable because the increase of the total flux cannot be appropriately considered in 2D FEM. In equal model specification total flux is recalculated at 3D FEM, because of considering lamination length. This paper compares the characteristic of Back-emf with and without overhang by permanent magnet that has identical volume.