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

Three dimensional turbulent flow fields around ships are simulated by a numerical method. Reynolds Averaged Navier-Stokes equations are used where Reynolds stresses are approximated by Baldwin-Lomax and Sub-Grid Scale(SGS) turbulence models. Body-fitted coordinate system is introduced to conform three dimensional ship geometries. The governing equations are discretized by a finite volume method. Temporal derivatives are approximated by the forward differencing and the convection terms are approximated by the QUICK or Kawamura scheme. The 2nd-order centered differencing is used for other spatial derivatives. Pressure and velocity fields are simultaneously iterated by the Highly Simplified Marker-And-Cell method. To verify the numerical method and turbulence models, flow fields around ships are simulated and compared to the experiments.

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

A finite element scheme using the concept of PISO method has been developed to solve the Navier-Stokes viscous flows in all speed range. This scheme includes development of new pressure equation that retains both the hyperbolic term related with the density variation and the elliptic term reflecting the incompressibility constraint. The present method is applied to the incompressible two-dimensional driven cavity flow problems(Re=100, 400 and 1,000). For compressible flows, the Carter plate problem(M=3 and Re=1,000) is computed. Finally, we have simulated the shock-boundary layer interaction(M=2 and Re=2.96×10/sup 5/), a more difficult problem, and compared its results with the experiment to demonstrate the shock capturing capability of the present solution algorithm.

• 한국전산유체공학회지
• /
• 제2권2호
• /
• pp.35-43
• /
• 1997

An efficient inverse method for 1.he supersonic/hypersonic axisymmetric body design is developed for the parabolized Navier-Stokes equations. The developed method is examined numerically for three extreme testcases in the supersonic(M/sub ∞/=3.0) and hypersonic(M/sub ∞/=6.28) speeds. The first one is a negative pressure distribution near a vacuum pressure and the second one is a positive pressure distribution over the whole region of the body. The last one is the case of abrupt change of pressure distribution to zero in the forward region of the body. These testcases show the robustness of the method. By introducing a regular-falsi method and by using a not-fully converged inverse solution, the convergence behavior was greatly improved.

• 한국전산유체공학회지
• /
• 제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.

• 한국추진공학회지
• /
• 제14권1호
• /
• pp.11-19
• /
• 2010

온도예조건화 나비어스톡스 방정식의 계산오차를 줄이기 위한 방법을 제시하였다. 이 방법은 또한 기존 예조건화 방법론에도 적용하였다. 제시된 방법의 타당성을 검토하기 위하여 다양한 마하수의 원형 실린더 주위의 단열 층류 점성 유동을 계산하였다. 엔탈피의 재정의를 통해 총엔탈피의 크기 정도를 줄임으로써 계산오차에 의한 온도예조건화의 수렴성 문제가 해결됨을 보였다.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제11권1호
• /
• pp.31-39
• /
• 2007

A pedagogic model for blood flow is introduced to help medicine majors understand a simplified version of Navier-Stokes equations which is known to be a good tool for interpreting the phenomena in blood flow. The pressure gradient consists of a time-independent part known as Hagen-Poiseuille's gradient and a time-dependent part known as Sexl's, and the model formula for the volume rate of blood flow is reduced to a very simple form. For demonstration, the blood rate in human aorta system is analyzed in connection with the time-dependence of pressure gradient. It is shown for Sexl's part that the flow rate lags the pressure gradient by ${\pi}/2$, which is thought to be due to the relaxation process involved.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제15권2호
• /
• pp.143-150
• /
• 2011

We consider the global existence of strong solutions of the 3D incompressible Navier-Stokes equations in a thin periodic domain. We present a simple proof that a strong solution exists globally in time when the initial velocity in $H^1$ and the forcing function in $L^p$(0,${\infty}$;$L^2$), $2{\leq}p{\leq}{\infty}$ satisfy certain condition. This condition is basically similar to that by Iftimie and Raugel[7], which covers larger and larger initial data and forcing functions as the thickness of the domain ${\epsilon}$ tends to zero.

• 한국전산유체공학회지
• /
• 제4권2호
• /
• pp.67-73
• /
• 1999

In the present paper, three types of the flow solvers were used to investigate the influence on the airfoil shape optimization. The adopted equations, i.e., Euler, thin layer Navier-Stokes and full Navier-Stokes ones. are solved using implicit LU-ADI decomposition scheme. The gradient projection method with the sinusoidal function was used as an optimization algorithm. The present numerical method was applied to the drag minimization problems under the initial shape of NACA0012 airfoils.

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

After the collapse of the Tacoma bay bridge at Tacoma Washington, the accurate prediction of aerodynamics became crucial to the sound design of bridges. CFD(Computational Fluid Dynamics) becomes important tool for the prediction on wind effects on the bridge due to the recent development of CFD. The usage of CFD is further prompted by the advantages in using CFD, such as low-cost and fast feed-back of design. In this paper, an unsteady compressible Reynolds averaged Navier-Stokes code is used for the computation of the flow over bridges. Coakley's ��q-${\omega}$ �� two-equation turbulence model is used for the turbulent eddy viscosity. For accurate and stable computations, the local preconditioning method is adapted to the code. Aerodynamic characteristics of a couple bridges are presented to show the validity and the accuracy of the method.