• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.104-111
• /
• 2009

The existing code which solves two-dimensional RANS(Reynolds Averaged Navier-Stokes) equations and 2-equation turbulence model equations was modified to enable numerical simulation of various supersonic flows. For this, various boundary conditions have been implemented to the code. Bleed boundary condition was incorporated into the code for calculating wall mean flow quantities. Furthermore, the boundary conditions for the turbulence quantities along rough surfaces as well as porous walls were applied to the code. The code was verified and validated by comparing the computational results against the experimental data for the supersonic flows over bleed region on a flat plate. Using the newly modified code, numerical simulations were performed and compared with other computational results as well as the experimental data for the supersonic flows over an oblique shock with a bleed region.

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

Autonomous underwater vehicles (AUVs) are unmanned, underwater vessels that are used to investigate sea environments in the study of oceanography. Docking systems are required to increase the capability of the AUVs, to recharge the batteries, and to transmit data in real time for specific underwater works, such as repented jobs at sea bed. This paper presents a visual :em control system used to dock an AUV into an underwater station. A camera mounted at the now center of the AUV is used to guide the AUV into dock. To create the visual servo control system, this paper derives an optical flow model of a camera, where the projected motions of the image plane are described with the rotational and translational velocities of the AUV. This paper combines the optical flow equation of the camera with the AUVs equation of motion, and deriver a state equation for the visual servo AUV. Further, this paper proposes a discrete-time MIMO controller, minimizing a cost function. The control inputs of the AUV are automatically generated with the projected target position on the CCD plane of the camera and with the AUVs motion. To demonstrate the effectiveness of the modeling and the control law of the visual servo AUV simulations on docking the AUV to a target station are performed with the 6-dof nonlinear equations of REMUS AUV and a CCD camera.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03b
• /
• pp.568-571
• /
• 2008

The computation of moving interface by the level set method typically requires reinitializations of level set function. An inaccurate estimation of level set function ${\phi}$ results in incorrect free-surface capturing and thus errors such as mass gain/loss. Therefore, accurate and robust reinitialization process is essential to the free-surface flows. In the present paper, we pursue further development of the reinitialization process, which evaluates directly level set function ${\phi}$ using a normal vector in the interface without solving the re-distancing equation of hyperbolic type. The Taylor-Galerkin approximation and P1P1splitting FEM are adopted to discretize advection equation of the level set function and the Navier-Stokes equation, respectively. Advection equation of free surface and re-initialization process are validated with benchmark problems, i.e., a broken dam flow and time-reversed single vortex flow. The simulation results are in good agreement with the existing results.

• Journal of Mechanical Science and Technology
• /
• v.15 no.3
• /
• pp.392-402
• /
• 2001

The paper reports a multiple source modeling of low-Reynolds-number dissipation rate equation with aids of DNS data. The key features of the model are to satisfy the wall limiting conditions of the individual source terms in the exact dissipation rate equation using the wall damping functions. The wall damping functions are formulated in term of dimensionless dissipation length scale ι(sup)+(sub)D(≡ι(sub)D($\upsilon$$\xi$)(sup)1/4/$\upsilon$) and the invariants of small and large scale turbulence anisotropy tensors. $\alpha$(sub)ij(=$\mu$(sub)i$\mu$(sub)j/$\kappa$-2$\delta$(sub)ij/3) and e(sub)ij(=$\xi$(sub)ij/$\xi$-2$\delta$(sub)ij/3). The model constants are optimized with aids of DNS data in a plane channel flow. Adopting the dissipation length scale as a parameter of damping function, the applicabilities of $\kappa$-$\xi$ model are extended to the turbulent flow calculation of complex flow passages.

• Journal of Advanced Marine Engineering and Technology
• /
• v.34 no.8
• /
• pp.1084-1093
• /
• 2010

Large eddy simulation of turbulent premixed flame stabilized by the bluff body is performed by using sub-grid scale combustion model based on the G-equation describing the flame front propagation. The basic idea of LES modeling is to evaluate the filtered-front speed, which should be enhanced in the grid scale by the scale fluctuations. The dynamic subgrid scale models newly introduced into the G-equation are validated by the premixed combustion flow behind the triangle flame holder. The calculated results can predict the velocity and temperature of the combustion flow in good agreement with the experiment data.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.39 no.1
• /
• pp.64-74
• /
• 1997

A CELTHYM(CEll-based Long-term HYdrologic Model), a pre-processor and a post-processor that can be integrated with geographic information system(GIS) were developed to predict the stream flow from the small agricultural watershed on the daily basis. The CELTHYM calculates the direct runoff from a grid using SCS curve number method and then sum up all of cells with respect to a sub-catchment area belonged to a stream grid and integrated to an outlet. Base flow of a watershed outlet was computed by integrating of the base flow of each stream grid that was averaged the sub-catchment deep-percolation and calculated with the release rate. Two kind of water budget equation were used to compute the water balance in a grid that was classified into not paddy field and paddy field. One of the two equation is a soil water balance equation to account the soil moisture of the upland, forest and excluding paddy field grid. The other is a paddy water balance equation for the paddy field, calculating the ponding depth, the effective rainfall, the deep percolation and the evapotranspiration.

• Journal of the Society of Naval Architects of Korea
• /
• v.43 no.2 s.146
• /
• pp.177-185
• /
• 2006

Cavitation is the dominant noise source of the marine vehicle. Of the various types of cavitation , tip vortex cavitation is the first appearance type of marine propeller cavitation and it generates high frequency noise. In this study, tip vortex cavitation behavior and noise are numerically investigated. A numerical scheme using Eulerian flow field computation and Lagrangian particle trace approach is applied to simulate the tip vortex cavitation on the hydrofoil. Vortex flow field is simulated by combined Moore and Saffman's vortex core radius equation and Sculley vortex model. Tip vortex cavitation behavior is analyzed by coupled Rayleigh-Plesset equation and trajectory equation. The cavitation nuclei are distributed and released in the vortex flow result. Vortex cavitation trajectories and radius variations are computed according to nuclei initial size. Noise is analyzed using time dependent cavitation bubble position and radius data. This study may lay the foundation for future work on vortex cavitation study and it will provide a basis for proper underwater propeller noise control strategies.

• Korea-Australia Rheology Journal
• /
• v.15 no.2
• /
• pp.83-90
• /
• 2003

In cases of the microfluidic channel, the electrokinetic influence on the transport behavior can be found. The externally applied body force originated from the electrostatic interaction between the nonlinear Poisson-Boltzmann field and the flow-induced electrical field is applied in the equation of motion. The electrostatic potential profile is computed a priori by applying the finite difference scheme, and an analytical solution to the Navier-Stokes equation of motion for slit-like microchannel is obtained via the Green's function. An explicit analytical expression for the induced electrokinetic potential is derived as functions of relevant physicochemical parameters. The effects of the electric double layer, the zeta potential of the solid surface, and the charge condition of the channel wall on the velocity profile as well as the electroviscous behavior are examined. With increases in either electric double layer or zeta potential, the average fluid velocity in the channel of same charge is entirely reduced, whereas the electroviscous effect becomes stronger. We observed an opposite behavior in the channel of opposite charge, where the attractive electrostatic interactions are presented.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.1 no.1
• /
• pp.15-22
• /
• 2000

There are lots of industrial machines of which functions are achieved by operation of multi-phase fluids. Some of them take advantage of the characteristics of counter-current two-phase flow The maximum flow rates of gas and liquid phases which flow in opposite-directions (counter-current flow) are limited by a phenomenon known as a Counter-Current Flow Limitation (CCFL or Flooding) The mass and momentum conservation equations for each Phase were established to build a first-order hyperbolic partial derivative equations system. A new CCFL model is developed based on the characteristic equation of the hyperbolic PDE system. The present model has its applicationto the case in which a non-uniform flow is developed around a square or sharp-edged entrance of liquid phase. The model is able to he used to Predict the operating-limit of components in which mass and heat transfer are taking place between liquid and gas phases.

• Structural Engineering and Mechanics
• /
• v.37 no.4
• /
• pp.401-413
• /
• 2011

The aerodynamic stability of orthotropic tensioned membrane structures with rectangular plane is theoretically studied under the uniform ideal potential flow. The aerodynamic force acting on the membrane surface is determined by the potential flow theory in fluid mechanics and the thin airfoil theory in aerodynamics. Then, based on the large amplitude theory and the D'Alembert's principle, the interaction governing equation of wind-structure is established. Under the circumstances of single mode response, the Bubnov-Galerkin approximate method is applied to transform the complicated interaction equation into a system of second order nonlinear differential equation with constant coefficients. Through judging the stability of the system characteristic equation, the critical divergence instability wind velocity is determined. Finally, from different parametric analysis, we can conclude that it has positive significance to consider the characteristics of orthotropic and large amplitude for preventing the instability destruction of structures.