• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.15-16
• /
• 2003

An implicit pressure correction method on unstructured Cartesian grid is developed for the incompressible Navier-Stokes equations. An immersed boundary method is also incorporated to treat the body geometry. Tests show that with an appropriate amount of dissipation, the method is second order accurate both in time and space. The driven cavity flows with and without immersed bodies are computed to demonstrate the capability of the present scheme.

• 한국방재학회:학술대회논문집
• /
• 2007.02a
• /
• pp.181-185
• /
• 2007

A free-surface correction(FSC) method is presented to solve the 3-D shallow water equations. Using the mode splitting process, FSC method can simulate shallow water flows under the hydrostatic assumption. For the hydrostatic pressure calculation, the momentum equations are firstly discretized using a semi-implicit scheme over the vertical direction leading to the tri-diagonal matrix systems. A semi-implicit scheme has been adopted to reduce the numerical instability caused by relatively small vertical length scale compare to horizontal one. and, as the free surface correction step the final horizontal velocity fields are corrected after the final surface elevations are obtained. Finally, the vertical final velocity fields can be calculated from the continuity equation. The numerical model is applied to the calculation of the simulation of flow fields in a rectangular open channel with the tidal influence. The comparisons with the analytical solutions show overall good agreements between the numerical results and analytical solutions.

• Journal of computational fluids engineering
• /
• v.13 no.3
• /
• pp.51-61
• /
• 2008

The SMAC (Simplified Marker And Cell) algorithm is extended for an application to thermal non-equilibrium two-phase flows in light water nuclear reactors (LWRs). A two-fluid three-field model is adopted and a multi-dimensional unstructured grid is used for complicated geometries. The phase change and the time derivative terms appearing in the continuity equations are implemented implicitly in a pressure correction equation. The energy equations are decoupled from the momentum equations for faster convergence. The verification of the present numerical method was carried out against a set of test problems which includes the single and the two-phase flows. The results are also compared to those of the semi-implicit ICE method, where the energy equations are coupled with the momentum equation for pressure correction.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.31-35
• /
• 2006

The delta-formulation of the Navier-Stokes equations has been popularly used in the aerodynamics area. Implicit algorithm can be easily implemented in that by using Taylor series expansion. This formulation is extended for an unsteady analysis by using a dual-time integration. In the meanwhile, the incompressible flows with heat transfers which occur in the area of thermo-hydraulics have been solved by a segregated algorithm such as the SIMPLE method, where each equation is discretised by using an under-relaxed deferred correction method and solved sequentially. In this study, the dual-time delta formulation is implemented in the segregated Navier-Stokes solver which is based on the collocated cell-centerd scheme with un unstructured mesh FVM. The pressure correction equation is derived by the SIMPLE method. From this study, it was found that the Euler dual-time method in the delta formulation can be combined with the SIMPLE method.

• International Journal of Aeronautical and Space Sciences
• /
• v.2 no.2
• /
• pp.28-38
• /
• 2001

A numerical method for the assessment and correction of tunnel wall interference effects on forced-oscillation testing is presented. The method is based on the wall pressure signature method using computed wall pressure distributions. The wall pressure field is computed using unsteady three-dimensional full Navier-Stokes solver for a 70-degree pitching delta wing in a wind tunnel. Approximately-factorized alternate direction implicit (AF-ADI) scheme is advanced in time by solving block tri-diagonal matrices. The algebraic Baldwin-Lomax turbulence, model is included to simulate the turbulent flow effect. Also, dual time sub-iteration with, local, time stepping is implemented to improve the convergence. The computed wall pressure field is then imposed as boundary conditions for Euler re-simulation to obtain the interference flow field. The static computation shows good agreement with experiments. The dynamic computation demonstrates reasonable physical phenomena with a good convergence history. The effects of the tunnel wall in upwash and blockage are analyzed using the computed interference flow field for several reduced frequencies and amplitudes. The corrected results by pressure signature method agree well with the results of free air conditions.

• Coupled systems mechanics
• /
• v.7 no.4
• /
• pp.421-440
• /
• 2018

In this study, a computational method for wall shear stress combined with an implicit direct-forcing immersed boundary method is presented. Near the immersed boundaries, the sub-grid stress is determined by a wall model in which the wall shear stress is directly calculated from the Lagrangian force on the immersed boundary. A coupling mathematical model of the transition process for a model Francis turbine comprising turbulent flow and rotating rigid guide vanes is established. The spatiotemporal distributions of pressure, velocity, vorticity and turbulent quantity are gained with the transient process; the drag and lift coefficients as well as other forces (moments) are also obtained as functions of the attack angle. At the same time, analysis is conducted of the characteristics of pressure pulsation, velocity stripes and vortex structure at some key parts of flowing passage. The coupling relations among the turbulent flow, the dynamical force (moment) response of blade and the rotating of guide vane are also obtained.

• Journal of Korean Port Research
• /
• v.13 no.1
• /
• pp.167-174
• /
• 1999

Numerical analysis of two-fluid flows for both water and air is carried out. Free-Surface flows with an arbitrary deformation have been simulated around two dimensional submerged hydrofoil. The computation is performed using a finite volume method with unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell-wise local mesh refinement. the integration in space is of second order based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels The linear equation systems are solved by conjugate gradient type solvers and the non-linearity of equations is accounted for through picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations the continuity equation the conservation equation of one species and the equations or two turbulence quantities.

• Journal of Mechanical Science and Technology
• /
• v.14 no.7
• /
• pp.789-795
• /
• 2000

Free-surface flows with an arbitrary deformation, induced by a submerged hydrofoil, are simulated numerically, considering two-fluid flows of both water and air. The computation is performed by a finite volume method using unstructured meshes and an interface capturing scheme to determine the shape of the free surface. The method uses control volumes with an arbitrary number of faces and allows cell wise local mesh refinement. The integration in space is of second order, based on midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation in time through three time levels. The linear equations are solved by conjugate gradient type solvers, and the non-linearity of equations is accounted for through Picard iterations. The solution method is of pressure-correction type and solves sequentially the linearized momentum equations, the continuity equation, the conservation equation of one species, and the equations for two turbulence quantities. Finally, a comparison is quantitatively made at the same speed between the computation and experiment in which the grid sensitivity is numerically checked.

• Journal of Advanced Marine Engineering and Technology
• /
• v.32 no.8
• /
• pp.1215-1220
• /
• 2008

To obtain the shape of the free surface more accurately, computations are carried out by a finite volume method using unstructured meshes and an interface capturing method. Free-surface flow, which is very important in the fields of ship and marine engineering, is numerically simulated for flows of both water and air. Control volumes are used with an arbitrary number of faces and allows a local mesh refinement. The integration is of second order, with a midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation. The solution method of pressure-correction type solves sequentially equations of momentum, continuity, conservation, and two-equations turbulence model. Comparison are quantitatively made between the computation and experiment in order to confirm the solution method.

• The Magazine of the Society of Air-Conditioning and Refrigerating Engineers of Korea
• /
• v.17 no.4
• /
• pp.382-394
• /
• 1988

Dynamic characteristics of thermally-forced stratification process in a square enclosure with a linear temperature profile at the side walls have been investigated through flow visualization experiment and numerical analysis. The experiment was performed on air with the Rayleigh numbers of order $10^5$. A particle tracer method is used for the flow visualization and to obtain a sudden linear temperature profile at the side walls copper blocks which already have a linear temperature profile are come into contact with the thin copper plates of the test section. Immediately a meridional circulation is developed and heat transfer takes place from the wall to the interior region by circulation of fluid and finally a thermal stratification is achieved. In the numerical study, QUICK scheme for convective terms, SIMPLE algorithm for pressure correction, and the implicit method for the time marching are adopted for the integration of conservation equations. Comparison of flow visualization and numerical results shows that the developing flow patterns are very similar in dynamic nature even though there is a time lag due to the inevitable time delay in setting up a linear temperature profile. For high Rayleigh numbers, the oscillatory motion is likely to take place and stratified region is extended. However, initial temperature adjustment process is much slower than that for low Rayleigh numbers.