• 대한기계학회논문집
• /
• 제19권8호
• /
• pp.1950-1963
• /
• 1995

A modification of the approximate-factorization method is made to accelerate the convergency rate and to take sufficiently large Courant number without loss of accuracy. And a stable implicit finite-difference scheme for solving the incompressible Navier-Stokes equations employed above modified method is developed. In the present implicit scheme, the volume fluxes with contravariant velocity components and the pressure formulation in curvilinear coordinates is adopted. In order to satisfy the continuity condition completely and to remove spurious errors for the pressure, the Navier-Stokes equations are solved by a modified SMAC scheme using a staggered gird. The upstream-difference scheme such as the QUICK scheme is also employed to the right hand side. The implicit scheme is unconditionally stable and satisfies a diagonally dominant condition for scalar diagonal linear systems of implicit operator on the left hand side. Numerical results for some test calculations of the two-dimensional flow in a square cavity and over a backward-facing step are obtained using both usual approximate-factorization method and the modified one, and compared with each other. It is shown that the present scheme allows a sufficiently large Courant number of O(10$^{2}$) and reduces the computing time.

• Coupled systems mechanics
• /
• 제7권4호
• /
• pp.407-420
• /
• 2018

In solid-liquid two phase flow, the knowledge of how descending solid particles affected by the presence of downstream wall is important. This work studies at what interstitial distance the velocity of a vertically descending sphere is affected by a downstream wall as a consequence of wall-modified hydrodynamic forces through a validated dynamic model. This interstitial distance-the hydrodynamic coupling distance ${\delta}_c-is$ found to decay monotonically with the approach Stokes number St which compares the particle inertia to viscous drag characterized by the quasi-steady Stokes' drag. The scaling relation ${\delta}_c-St-1$ decays monotonically as literature below the value of St equal to 10. However, the faster diminishing rate is found above the threshold value from St=10-40. Furthermore, an empirical relation of ${\delta}_c-St$ shows dependence on the drop height which clearly indicates the non-negligible effect of unsteady hydrodynamic force components, namely the added mass force and the history force. Finally, we attempt a fitting relation which embedded the particle acceleration effect in the dependence of fitting constants on the diameter-scaled drop height.

• 한국전산유체공학회지
• /
• 제10권1호
• /
• pp.67-72
• /
• 2005

This research is based on the iterative space-marching method for incompressible and compressible Navier-Stokes equations[1-4]. A principle of parallel computational schemes construction for steady and unsteady problems is suggested. It is analytically proven that convergence of these schemes is unconditional for incompressible case. When the parallel scheme is used the total volume of computations is the sum of a large number of independent and equal parts. Estimation of the speed-up K shows that K > 1000 in ideal case. First results of using the parallel schemes are presented.

• 대한기계학회논문집B
• /
• 제21권1호
• /
• pp.49-56
• /
• 1997

Heat-transfer enhancement is seeked through modifications of fin surface. Real life plate-fin heat exchangers have complex three-dimensional geometries. Fins can have arrays of dimples and are attached to rows of penetrating tubes. To isolate the effect of surface modification, we model the real flow by a two-dimensional channel flow with a dimple on one side. The flow is analysed by solving the incompressible Navier-Stokes equation by a finite volume method on a generalized boundary-fitted coordinate. Results show a trapped vortex inside the dimple for all cases computed. Local maximum of Nusselt number occurs near the downstream end of the dimple, due to such a vortex. Location of the vortex does not change with respect to the wall temperature change, but moved downstream when Reynolds number increases. This, together with the results that in all cases vortex core is somewhat downstream of the dimple center, suggests that the mean flow above continuously feeds the kinetic energy to the recirculating flow. Heat transfer enhancement and pressure losses are studied through analysing the relevant dimensionless parameters like, Nusselt number and friction factor. In all cases computed, dimpled channel flow experiences less pressure loss than two-dimensional Poiseuille flow.

• 대한기계학회논문집B
• /
• 제26권8호
• /
• pp.1088-1094
• /
• 2002

A channel flow with a high Reynolds number but coarse grids is numerically studied to investigate the prediction possibility of its turbulence which is three-dimensional and time-dependent. In the present paper, a Reynolds-Averaged Navier-Stokes (RANS) model, a Large Eddy Simulation (LES) and a Navier-Stokes equation with no model are tested with a new approach of hybrid RANS/LES, which reduces to RANS model in the boundary layers and at separation, and to Smagorinsky-like LES downstream of separation, and then compared with each other. It is found that the simulations of hybrid RANS/LES method sustain turbulence like those of LES and with no model, and the results are stable and fairly accurate. This indicates strongly that gradual improvements could lead to a simple, stable, and accurate approach to predict turbulence phenomena of wall-bounded flow.

• 대한기계학회논문집B
• /
• 제21권2호
• /
• pp.202-212
• /
• 1997

In this study, performance of the multi-stage explicit methods for numerical computation of the unsteady Navier-Stokes equations is investigated. Three methods under consideration are 1 st-, 2 nd-, and 4 th-order Runge-Kutta (R-K) methods. Compared in this estimation is stability, accuracy, and CPU time of each method. The computational codes developed are applied to the two-dimensional flow in a square cavity driven by an oscillating lid. It turned out that at Reynolds number 400, the 1 st-order R-K method is the best, while at 3200 the 2 nd-order R-K is recommended. At higher Reynolds numbers, it is conjectured that the 4 th-order R-K method will be the best algorithm among three due to its highest stability.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2001년도 춘계 학술대회논문집
• /
• pp.149-154
• /
• 2001

Parallel performance of a Myrinet based PC-cluster was tested and compared with a conventional Fast-Ethernet system. A preconditioned Navier-Stokes code was parallelized with domain decomposition technique, and used for the parallel performance test. Speed-up ratio was examined as a major performance parameter depending on the number of processor and the network topology. As was expected, Myrinet system shows a superior parallel performance to the Fast-Ethernet system even with a single network adpater for a dual processor SMP machine. A test for the dependency on problem size also shows that network communication speed is a crucial factor for parallelized computational fluid dynamics analysis and the Myrinet system is a plausible candidate for high performance parallel computing system.

• Journal of Astronomy and Space Sciences
• /
• 제29권1호
• /
• pp.85-91
• /
• 2012

Optical interferometry and polarimetry have separately provided new insights into stellar astronomy, especially in the fields of fundamental parameters and atmospheric models. We present: scientific justifications for "full-Stokes" optical interferometric polarimetry (OIP); updated instrument requirements; preliminary beam combiner designs; polarimeter design; end-to-end OIP data reduction; and realistic reimaged full-Stokes models of Be stars with a suitable number of telescopes plus noise sources. All of this work represents preliminary research to construct an OIP beam combiner.

• 한국전산유체공학회지
• /
• 제1권1호
• /
• pp.123-141
• /
• 1996

In this paper the upwind flux difference splitting Navier-Stokes method has been applied to study the perfect gas and the equilibrium chemically reacting hypersonic flow over an axisymmetric sphere-cone(5°) geometry. The effective gamma(γ), enthalpy to internal energy ratio was used to couple chemistry with the fluid mechanics for equilibrium chemically reacting air. The test case condition was at altitude(30km) and Mach number(15). The equilibrium shock thickness over the blunt body region was much thinner than that of perfect gas shock. The pressure difference between perfect gas and equilibrium gas was about 3 ∼ 5 percent. The heat transfer coefficient were also calculated. The results were compared with VSL results in order to validate the current numerical analysis. The results from current method were compared well VSL results ; however, not well at near nose. The proper boundary condition and grid system will be studied to improve the solution quality.

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

In this paper the upwind flux difference splitting Navier-Stokes method has been applied to study the perfect gas and the equilibrium chemically reacting hypersonic flow over an axisymmetric sphere-cone($5^{\circ}$) geometry. The effective gamma($\bar{r}$), enthalpy to internal energy ratio was used to couple chemistry with the fluid mechanics for equilibrium chemically reacting air. The test case condition was at altitude(30Km) and Mach number(15). The equilibrium shock thickness over the blunt body region was much thinner than that of perfect gas shock. The pressure difference between perfect gas and equilibrium gas was about $3\sim5$ percent. The skin friction coefficient and heat transfer coefficient were also calculated.