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

A numerical simulation method has been under development for solving turbulent flows around a ship model in maneuvering motion using the Reynolds Averaged Navier-Stokes equations. The method used second-order finite differences, collocated grids, pressure-Poisson equation and four-stage Runge-Kutta scheme as key components of the solution method. A modified Baldwin-Lomax model is used for the turbulence closure. This paper presents a preliminary result of the computational study on turbulent flows past a ship model in drift motion. Calculations are carried out for a Series 60 $C_B=0.6$ ship model, for which detailed experimental data are available. The results of the present calculations are compared with the experimental data for hydrodynamic forces acting on the model as well as velocity distributions at longitudinal sections. Only fair agreements has been achieved. The computational results show the complex asymmetrical shear flow patterns including three-dimensional separations followed by formation of bilge vortices both in bow and stern regions.

• Geomechanics and Engineering
• /
• v.17 no.1
• /
• pp.11-18
• /
• 2019

The formulas offered so far on the settlement of raft footings provide only a rough estimate of the actual settlement. One of the best ways to make an accurate estimation is to conduct 3-dimensional finite element analyses. However, the required procedure for these analyses is comparatively cumbersome and expensive and needs a bit more expertise. In order to address this issue, in this study, a raft footing settlement formula was developed based on ninety finite element model configurations. The formula was derived using multi-parameter exponential regression analyses. The settlement formula incorporates the dimensions and the elastic modulus of a rectangular raft, vertical uniform pressure and soil moduli and Poisson's ratios up to 5 layers. In addition to this, an equation was offered for the estimation of average deflection of the raft. The proposed formula was checked against 3 well-documented case studies. The formula that is derived from 3D finite element analyses is useful in optimising the raft properties.

• Smart Structures and Systems
• /
• v.8 no.5
• /
• pp.433-447
• /
• 2011

The present paper addresses the nonlinear response of a FG square plate with two smart layers as a sensor and actuator under pressure. Geometric nonlinearity was considered in the strain-displacement relation based on the Von-Karman assumption. All the mechanical and electrical properties except Poisson's ratio can vary continuously along the thickness of the plate based on a power function. Electric potential was assumed as a quadratic function along the thickness direction and trigonometric function along the planar coordinate. By evaluating the mechanical and electrical energy, the total energy equation can be minimized with respect to amplitude of displacements and electrical potential. The effect of non homogenous index was investigated on the responses of the system. Obtained results indicate that with increasing the non homogenous index, the displacements and electric potential tend to an asymptotic value. Displacements and electric potential can be presented in terms of planar coordinate system. A linear analysis was employed and then the achieved results are compared with those results that are obtained using the nonlinear analysis. The effect of the geometric nonlinearity is investigated by using the comparison between the linear and nonlinear results. Displacement-load and potential-load curves verified the necessity of a nonlinear analysis rather than a linear analysis. Improvement of the previous results (by the linear analysis) through employing a nonlinear analysis can be presented as novelty of this study.

• Journal of the Korean Vacuum Society
• /
• v.17 no.6
• /
• pp.487-494
• /
• 2008

The numerical simulation of point-to-plane discharge of nitrogen at medium pressure has been achieved by a quasi-2d numerical model. In the model, we calculate the distributions of electric charges which are varying as temporal and space and determine the electric field depending on space charge distribution by solving Poisson's equation. The continuity equations are treated numerically by using FCT (Flux-Corrected Transport) Algorithm and FEM (Finite Element Method). The numerical simulation results make us to understand the physical characteristics of nitrogen discharge at 50 torr. The comparison with experimental results[1] shows a good qualitative agreement.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.2
• /
• pp.70-76
• /
• 2000

A numerical study on the viscous flow around a 2-dimensional circulation control foil is carried out for application on the field of naval architecture and ocean engineering. The governing equations are the RANS and the continuity equations. The equations are discretized by finite difference method and MAC method and the pressure poisson equation is calculate by a SOR method and an O-type non-staggered boundary fitted coordinate system which is overlapped near the slot is used to improve the numerical accuracy. Turbulence is approximated by a modified Baldwin-Lomax turbulence model. In the present paper, the Coanda effect on a 2-dimensional foil of a 20% thickness ellipse with modified rounded trailing edge has been numerically studied. The change in drag and lift of the foil with various jet momentums are calculated and compared to the experimental results to show good agreements.

• Wind and Structures
• /
• v.14 no.5
• /
• pp.465-480
• /
• 2011

In this paper, a stabilized large eddy simulation technique is developed to predict turbulent flow with high Reynolds number. Streamline Upwind Petrov-Galerkin (SUPG) stabilized method and three-step technique are both implemented for the finite element formulation of Smagorinsky sub-grid scale (SGS) model. Temporal discretization is performed using three-step technique with viscous term treated implicitly. And the pressure is computed from Poisson equation derived from the incompressible condition. Then two numerical examples of turbulent flow with high Reynolds number are discussed. One is lid driven flow at Re = $10^5$ in a triangular cavity, the other is turbulent flow past a square cylinder at Re = 22000. Results show that the present technique can effectively suppress the instabilities of turbulent flow caused by traditional FEM and well predict the unsteady flow even with coarse mesh.

• Journal of the Korean Magnetics Society
• /
• v.4 no.1
• /
• pp.56-61
• /
• 1994

Magnetic fluid is ferromagnetic material in liquid state, so the surface configuration of magnetic fluid affects the magnetic field, and vice versa. To analyze the devices with magnetic fluid, the magnetic field equations and hydrodynamic equation should be solved simultaneously. This paper presents the numerical algorithm to obtain the surface configuration of fluid under the influence of gravity, pressure and magnetic field without conventional sim¬plified assumption. The algorithm consists of nonlinear finite element method and ferro-hydrodynamics, such as Poisson equations and Bernoulli equations, respectively The simulated configurations of fluid are compared with experimental results, and the influence of the amount of fluid and pole piece shape on the seal capacities are analyzed.

• Journal of Korea Water Resources Association
• /
• v.43 no.9
• /
• pp.801-811
• /
• 2010

A three-dimensional numerical model called NEWTANK is employed to investigate solitary wave run-up with an internal wave-maker on a steep slope. The numerical model solves the spatially averaged Navier-Stokes equations for two-phase flows. The LES (large-eddy-simulation) approach is adopted to model the turbulence effect by using the Smagorinsky SGS (sub-grid scale) closure model. A two-step projection method is adopted in numerical solutions, aided by the Bi-CGSTAB (Bi-Conjugate Gradient Stabilized) method to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate VOF (volume-of-fluid) method is used to track the distorted and broken free surface. A solitary wave is first internally generated and propagated over a constant water depth in the three-dimensional domain. Numerically predicted results are compared with analytical solutions and numerical errors are analyzed in detail. The model is then applied to study solitary wave run-up on a steep slope and the obtained results are compared with available laboratory measurements.

• Journal of Wetlands Research
• /
• v.16 no.3
• /
• pp.327-336
• /
• 2014

Submerged thin walls are extreme case of submerged rectangular blocks, and could be used for many purposes in rivers or coastal zones, e.g. to tsunami. To understand flow characteristics including flow and pressure fields around a specific submerged thin wall a numerical model was applied which includes computation of hydrodynamic pressure on ${\sigma}$-coordinate. ${\sigma}$-coordinate has strong merits for simulation of subcritical flow over mild-sloped beds. On the other hand ${\sigma}$-coordinate is quite poor to treat sharp structures on the bed. There have been a few trials to incorporate dynamic pressure in ${\sigma}$-coordinate by some researchers. One of the previous approaches includes process of sloving the Poisson equation. However, the above method includes many high-order terms, and requires long cpu for simulation. Another method SOLA was developed by Hirt et al. for computation of dynamic pressure, but it was valid for straight grid system only. Previous SOLA was modified for ${\sigma}$-coordinate for the present purpose and was adopted in a model system, CST3D. Computed flow field shows reasonable behaviour including vorticity is much stronger than the upstream and downstream of the structure. The model was verified to laboratory experiments at a 2DV flume. Time-average flow vectors were measured by using one-dimensional electro-magnetic velocimeter. Computed flow field agrees well with the measured flow field within 10 % error from the speed point of view at 5 profiles. It is thought that the modified SOLA scheme is useful for ${\sigma}$-coordinate system.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.5
• /
• pp.1158-1174
• /
• 1988

A numerical technique is developed for the solution of fully developed turbulent recirculating flow in the passage of variable area using the non-orthogonal coordinate transformation. In the numerical analysis, primitive pressure-velocity finite difference equations were solved by SIMPLER algorithm with 2-equation turbulence model and algebraic stress model (ASM). QUICK scheme on the differencing of convective terms which is free from the inaccuracies of numerical diffusion has been applied to the variable grids and the results compared with those from HYBRID scheme. In order to test the effect of streamline curvatures on turbulent diffusion Lee and Choi streamline curvature correction model which has been obtained by modifying the Leschziner and Rodi's model is testes. The ASM was also employed and the results are compared to those from another turbulence model. The results show that difference of convective differencing schemes and turbulence models give significant differences in the prediction of velocity fields in the expansion region and outlet region of the combustor, however show little differences in the parallel flow region.