• 대한기계학회논문집B
• /
• 제28권7호
• /
• pp.809-817
• /
• 2004

In this paper, the unsteady behavior of the viscous flow field past an impulsively started elliptic cylinder is studied numerically. In order to analyze flow field, we introduce vortex particle method. The vorticity transport equation is solved by fractional step algorithm which splits into convection term and diffusion term. The convection term is calculated with Biot-Savart law, the no-through boundary condition is employed on solid boundaries. The diffusion term is modified based on the scheme of particle strength exchange. The particle redistributed scheme for general geometry is adapted. The flows around an elliptic cylinder are investigated for various attack angles at Reynolds number 200. The comparison between numerical results of present study and experimental data shows good agreements.

• 대한토목학회논문집
• /
• 제31권3B호
• /
• pp.293-303
• /
• 2011

비압축성 점성 흐름을 수치해석하기 위한 효율적인 대각행렬화된 근사 인수분해(DAF) 알고리즘을 개발하였다. 압력에 근거한 인공압축성(AC) 기법을 이용하여 3차원 정상 비압축성 Navier-Stokes 방정식을 계산한다. AC 형태로 변형된 지배방정식은 2차 정확도의 유한차분법을 이용하여 공간에 대해서 이산화하였다. 이산화된 방정식계를 2차 정확도로 분할하기 위해서 본 연구에서 개발한 DAF 기법을 적용한다. 이 연구의 목적은 이 DAF 기법의 계산상 효율성을 검토하는 것이다. 만곡부를 갖는 사각형 덕트에서 완전히 발달한 층류 흐름과 발달하는 층류흐름 그리고 공동에서의 층류흐름에 대한 DAF 기법의 해석결과를 잘 알려진 4단계 Runge-Kutta(RK4)기법에 의한 해석해와 상대적으로 비교평가 하였다. 공간에 대해서 동일한 이산화기법을 이용하므로 동일한 격자상에서 계산된 DAF기법과 RK4기법의 해는 근본적으로 동일한 반면에, 이들 두기법의 계산상 효율성은 확연히 다른 것으로 나타났다. 본 연구에서 개발된 DAF기법은 적용한 모든 흐름 문제에 대해서 RK4기법에 비해 최소 2배 이상 적은 계산 시간만을 필요로 하는 것으로 나타났다. 이러한 DAF 기법의 계산상 효율성은 계산용량의 추가나 프로그래밍의 추가적인 복잡함이 없이 확보된다.

• Advances in aircraft and spacecraft science
• /
• 제5권3호
• /
• pp.363-383
• /
• 2018

One-dimensional (1D) models of incompressible flows, can be of interest for many applications in which fast resolution times are demanded, such as fluid-structure interaction of flows in compliant pipes and hemodynamics. This work proposes a higher-order 1D theory for the flow-field analysis of incompressible, laminar, and viscous fluids in rigid pipes. This methodology is developed in the domain of the Carrera Unified Formulation (CUF), which was first employed in structural mechanics. In the framework of 1D modelling, CUF allows to express the primary variables (i.e., velocity and pressure fields in the case of incompressible flows) as arbitrary expansions of the generalized unknowns, which are functions of the 1D computational domain coordinate. As a consequence, the governing equations can be expressed in terms of fundamental nuclei, which are invariant of the theory approximation order. Several numerical examples are considered for validating this novel methodology, including simple Poiseuille flows in circular pipes and more complex velocity/pressure profiles of Stokes fluids into non-conventional computational domains. The attention is mainly focused on the use of hierarchical McLaurin polynomials as well as piece-wise nonlocal Lagrange expansions of the generalized unknowns across the pipe section. The preliminary results show the great advantages in terms of computational costs of the proposed method. Furthermore, they provide enough confidence for future extensions to more complex fluid-dynamics problems and fluid-structure interaction analysis.

• Kyungpook Mathematical Journal
• /
• 제57권4호
• /
• pp.651-665
• /
• 2017

Some computational fluid dynamics problems concerning the thin films flow of viscous fluid with a free surface and draining or coating fluid-flow problems can be delineated by third-order ordinary differential equations. In this paper, the aim is to introduce the numerical solutions of the boundary value problems of such equations by variational iteration method. In this paper, it is shown that the third-order boundary value problems can be written as a system of integral equations, which can be solved by using the variational iteration method. These solutions are gleaned in terms of convergent series. Numerical examples are given to depict the method and their convergence.

• 한국전산유체공학회지
• /
• 제6권2호
• /
• pp.1-8
• /
• 2001

Despite the massive effort which has been given to the analysis of the base flows, one commonly occurring case seems to be overlooked. This is for base (rearward facing surface) which is between a subsonic flow and supersonic flow. Potential flows of the air and gas streams are computed for the flow past a separated wake. Then a viscous jet mixing is superimposed on this inviscid solution. Conservation of mass, momentum and energy is achieved by multiple iterations. Despite the iterations, the wake flow field is computed with modest computer requirements.

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

A numerical investigation has been made for flows in an axisymmetric circular cylinder with a rotating cone located at the bottom of the container. The axisymmetric container is completely filled with a viscous fluid. Major parameter for the present research is the vertex angle of the cone, otherwise Reynolds number of fluid and aspect ratio of the vessel is fixed. Main interest is in vortex breakdown of meridional circulation by rotation of the cone with respect to the longitudinal axis of the cylinder. The method to this problem is numerically to integrate momentum and continuity equations on a generalized body fitted grid system. The pattern of vortex breakdown is quite different from that in a right circular cylinder with flat end wall disks. Flow visualization photographs of a preceeding work are compared with the present numerical results.

• 한국해양공학회지
• /
• 제15권4호
• /
• pp.46-52
• /
• 2001

SMAC method, one of the numerical simulation techniques, is modified from the original MAC method for the time-dependent variation of fluid flows. The Navier-Stokes equations for incompressible time-dependent viscous flow are applied, and Also marker particles which present the visualization of fluid flows are used. In this study, two-dimensional numerical simulations of the overtopping are carried out by SMAC method, and the simulation results are visualized, In addition to, motion pictures are made for efficient visualization of the simulation results. This numerical simulation could also be applied to the design of coastal structures as dike and revetment.

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

An efficient Newton-GMRES algorithm is presented for computing two-dimensional steady compressible viscous flows on unstructured hybrid meshes. The scheme is designed on cell-centered finite volume method which accepts general polygonal meshes. Steady-state solution is obtained with pseudo-transient continuation strategy. The preconditioned, restarted general minimum residual(GMRES) method is employed in matrix-free form to solve the linear system arising at each Newton iteration. The incomplete LU fartorization is employed for the preconditioning of linear system. The Spalart-Allmars one equation turbulence model is fully coupled with the flow equations to simulate turbulence effect. The accuracy, efficiency and robustness of the presently developed method are demonstrated on various test problems including laminar and turbulent flows over flat plate and airfoils.

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

Aerodynamic sensitivity analysis is performed for the Navier-Stokes equations coupled with two-equation turbulence models using a discrete adjoint method and a direct differentiation method respectively. Like the mean flow equations, the turbulence model equations are also hand-differentiated to accurately calculate the sensitivity derivatives of flow quantities with respect to design variables in turbulent viscous flows. Both the direct differentiation code and the adjoint variable code adopt the same time integration scheme with the flow solver to efficiently solve the differentiated equations. The sensitivity codes are then compared with the flow solver in terms of solution accuracy, computing time and computer memory requirements. The sensitivity derivatives obtained from the sensitivity codes with different turbulence models are compared with each other. Using two-equation turbulence models, it is observed that a usual assumption of constant turbulent eddy viscosity in adjoint methods may lead to seriously inaccurate results in highly turbulent flows.

• 한국안전학회지
• /
• 제13권4호
• /
• pp.3-8
• /
• 1998

The vortex particle method, which includes viscous effects, consists of diffusion of boundary vorticity and creation of the vortex particles, convection, particle strength exchange, and particle redistribution. Accuracy of the boundary element method is very important since it creates the particles around the body at every time step. A boundary element method based on source panel was investigated as part of computation of unsteady separated flows by rising the vortex particle method. The potential flows were computed around a circular cylinder and a square cylinder. The results around the circular cylinder were compared with the exact solution, and the distribution of vorticity, in particular near the sharp comers of the square cylinder, is scrutinized for different number of panels.