• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2000년도 추계 학술대회논문집
• /
• pp.129-134
• /
• 2000

An efficient numerical method to solve the unsteady incompressible Navier-Stokes equations is developed. A fully implicit time advancement is employed to avoid the CFL(Courant-Friedrichs-Lewy) restriction, where the Crank-Nicholson discretization is used for both the diffusion and convection terms. Based on a block LU decomposition, velocity-pressure decoupling is achieved in conjunction with the approximate factorization. Main emphasis is placed on the additional decoupling of the intermediate velocity components with only n th time step velocity The temporal second-order accuracy is Preserved with the approximate factorization without any modification of boundary conditions. Since the decoupled momentum equations are solved without iteration, the computational time is reduced significantly. The present decoupling method is validated by solving the turbulent minimal channel flow unit.

• 대한기계학회논문집B
• /
• 제24권10호
• /
• pp.1317-1325
• /
• 2000

A new efficient numerical method for computing three-dimensional, unsteady, incompressible flows is presented. To eliminate the restriction of CFL condition, a fully-implicit time advancement in which the Crank-Nicolson method is used for both the diffusion and convection terms, is adopted. Based on an approximate block LU decomposition method, the velocity -pressure decoupling is achieved. The additional decoupling of the intermediate velocity components in the convection term is made for the fully -implicit time advancement scheme. Since the iterative procedures for the momentum equations are not required, the velocity components decouplings bring forth the reduction of computational cost. The second-order accuracy in time of the present numerical algorithm is ascertained by computing decaying vortices. The present decoupling method is applied to minimal channel flow unit with DNS (Direct Numerical Simulation).

• 산업기술연구
• /
• 제20권A호
• /
• pp.269-278
• /
• 2000

The pressure-time profile of the explosion gases can controlled for the use of cartridge explosive with two techniques known as Decoupling and spacing of the charges. Decoupling consists of a space between the explosive column and wall of the blast hole. Four different decoupling index 1.4, 1.8, 2.34, 3.0 are selected in this field study. The level of ground vibrations with each decoupling index was measured and the empirical particle velocity equation from these data was obtained. The condition of new cracks at blast hole are also examined. As the decoupling index is increased, the level of the blast vibration is decreased. But the cracks in rock masses are efficiently formed to remove the broken rock. The vibration constant associated with test sites is given as $K=1564.5(D.L)^{-1.3233}$ in terms of D.I.(decoupling index).

• 화약ㆍ발파
• /
• 제15권3호
• /
• pp.20-32
• /
• 1997

The pressure-time profile of the explosion gases can be controlled fot the use of cartridge explosives with two techniques Known as Decoupling and Spacing the charges. Decoupling consists in leaving and empty space between the explosive column and wall of the blast hole. Four different decoupling index, 1.4, 1.8, 2.34, 3.0 are selected in this field study. The level of ground vibrations with each decoupling index are measured and the empirical particle vibrations with each decoupling index are measured and the empirical particle velocity equation from these data was obtained. The condition of new cracks at blast hole are also examined. As the decoupling index in increased, the level of the blast vibration is decreased,. But the cracks in rock masses are efficiently formed to remove the broken rock. The vibration constant associated with a given site $K=1564.5(D.I)^{-1.3233}$ in terms of D.I(decopling index).

• 한국전산유체공학회지
• /
• 제7권1호
• /
• pp.10-19
• /
• 2002

The non-staggered(collocated) grid approach in which all the solution variables are located at the centers of control volumes is very popular for incompressible flow analyses because of its numerical efficiency on the curvilinear or unstructured grids. Rhie and Chow's paper is the first in using non-staggered grid method for SIMPLE algorithm, where pressure weighted interpolation was used to prevent decoupling of pressure and velocity. But it has been known that this non-staggered grid method has stability problems when pressure fields are nonlinear like in natural convection flows. Also Rhie-Chow scheme generates large numerical diffusion near curved walls. The cause of these unwanted problems is too large pressure damping term compared to the magnitude of face velocity. In this study the magnitude of pressure damping term of Rhie-Chow's method is limited to 1∼10% of face velocity to prevent physically unreasonable solutions. The wall pressure extrapolation which is necessary for cell-centered FVM is another source of numerical errors. Some methods are applied in a unstructured FV solver and analyzed in view of numerical accuracy. Here, two natural convection problems are solved to check the effect of the Rhie-Chow's method on numerical stability. And numerical diffusion from Rhie-Chow's method is studied by solving the inviscid flow around a circular cylinder.

• 한국수자원학회논문집
• /
• 제44권6호
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.

• 멤브레인
• /
• 제6권2호
• /
• pp.86-95
• /
• 1996

본 연구에서는1축 회전판형 한외여과막 모듈의 순수 투과율 예측모델과 기름 에멀션의 분리 특성 및 투과율 예측모델을 유도하였다. 1축 RDM은 한외여과막(UOP사, 직경 0.22m)으로 제작된 회전판막을 장착하여 $25^{\circ}C$에서 압력, 각속도($\omega$), 농도 변화에 따라 분리특성을 조사하였다. 회전판막의 각속도가 0에서 52.36rad/s로 증가할수록 회전판내 유체가 받는 원심력에 의한 압력 강하와 분리막 표면의 미끄럼 흐름에 의한 압력 강하로 순수 투과율은 최대 각속도에서 3.9% 감소하였다. 원심력과 미끄럼 흐름에 의한 압력 강하는 선속도(${\omega}r$)의 자승에 비례하였다. 회전판의 각속도가 52.36에서 2.62rad/s로 감소할 때 5% 절삭유의 투과율은 30.16% 감소하였고 농도가 낮을수록 막회전과 투과율에 미치는 영향은 적었다. 절삭유의 투과율(J; $kg/m^{2} \cdot s$)은 회전에 의한 압력 강하를 고려한 유효압력차($\Delta P_{T}$ ; Pa), 벌크농도($C_{B}$; %), 선속도($\omega$r; m/s) 등에 영향을 받으며 실험 결과에 저항 모델을 적용하여 다음과 같은 식을 유도하였다.