• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.7
• /
• pp.907-915
• /
• 2000

The preconditioned Krylov subspace methods were applied to the incompressible Navier-Stoke's equations for convergence acceleration. Three of the Krylov subspace methods combined with the five of the preconditioners were tested to solve the lid-driven cavity flow problem. The MILU preconditioned CG method showed very fast and stable convergency. The combination of GMRES/MILU-CG solver for momentum and pressure correction equations was found less dependency on the number of the grid points among them. A guide line for stopping inner iterations for each equation is offered.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.11 no.4
• /
• pp.55-68
• /
• 2007

First-order least-squares method of a distributed optimal control problem for the incompressible Navier-Stokes equations is considered. An optimality system for the optimal solution are reformulated to the equivalent first-order system by introducing velocity-flux variables and then the least-squares functional corresponding to the system is defined in terms of the sum of the squared $L^2$ norm of the residual equations of the system. The optimal error estimates for least-squares finite element approximations are obtained.

• Nuclear Engineering and Technology
• /
• v.41 no.9
• /
• pp.1191-1204
• /
• 2009

In this study, a CFD methodology is employed to address the problem of the prediction of the flow in a T-junction. An Unsteady Reynolds Averaged Navier-Stokes (URANS) approach has been selected for its low computational cost. Moreover, Unsteady Reynolds Navier-Stokes methodologies do not need complex boundary formulations for the inlet and the outlet such as those required when using Large Eddy Simulation (LES) or Direct Numerical Simulation (DNS). The results are compared with experimental data and an LES calculation. In the past, URANS has been tried on T-junctions with mixed results. The biggest limit observed was the underestimation of the oscillatory behavior of the temperature. In the present work, we propose a comprehensive approach able to correctly reproduce the root mean square (RMS) of the temperature directly downstream of the T-junction for cases where buoyancy is not present.

• Journal of applied mathematics & informatics
• /
• v.37 no.5_6
• /
• pp.507-519
• /
• 2019

The ground for this paper is to examine the generalized Stokes' first and second issues for an incompressible couple pressure liquid under isothermal conditions. Exact solutions for each problem are acquired by using the Laplace transform (LT) with respect to the time variable t and the sine Fourier transform (FT) with respect to the y-variable. Further, a comparison is given of the obtained results and the results of Devakar and Lyengar [1] and by using the four inverse Laplace transform algorithms (Stehfest's, Tzou's, Talbot, Fourier series) in the space time domain utilizing a numerical methodology. Moreover, velocity profiles are plotted and considered for various occasions and distinctive estimations of couple stress parameters. At the end, the outcomes are exhibited by graphs and in tabular forms.

• Proceedings of the Korean Institute of Navigation and Port Research Conference
• /
• v.1
• /
• pp.181-186
• /
• 2006

In this paper, generation mechanisms of ocean freak waves are briefly introduced in the context of wave-current interaction phenomena. The present model of the fluid motion is based on the Navier-Stokes equations incorporating velocity-pressure formulation because of need to model the nonlinear wave interaction with spatially non-uniform current field. In order to deal with the free surface motion, an Arbitrary Lagrangian-Eulerian (ALE) description is adopted. As an accurate and efficient numerical tool, the spectral element method is presented with general features and specific treatment for the wave-current interaction problem. As an intermediate stage of development, solution procedure and characteristics aspects of the present modeling and numerical method are addressed in detail, and preliminary numerical results prove its accuracy and convergence.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.6
• /
• pp.751-759
• /
• 2007

We have developed several methods for the optimization problem having large-scale and highly nonlinear system. First, step by step method in optimization process was employed to improve the convergence. In addition, techniques of furnishing good initial guesses for analysis using sensitivity information acquired from optimization iteration, and of manipulating analysis/optimization convergency criterion motivated from simultaneous technique were used. We applied them to flow control problem and verified their efficiency and robustness. However, they are based on quasi-Newton method that approximate the Hessian matrix using exact first derivatives. However solution of the Navier-Stokes equations are very cost, so we want to improve the efficiency of the optimization algorithm as much as possible. Thus we develop a true Newton method that uses exact Hessian matrix. And we apply that to the three-dimensional problem of flow around a sphere. This problem is certainly intractable with existing methods for optimal flow control. However, we can attack such problems with the methods that we developed previously and true Newton method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.3
• /
• pp.170-177
• /
• 2009

In the present study, a finite element analysis of conjugate heat transfer problem inside a cavity with a heat-generating conducting body, where constant heat flux is generated, is conducted. A conduction heat transfer problem inside the solid body is automatically coupled with natural convection inside the cavity by using a finite element formulation. A finite element formulation based on SIMPLE type algorithm is adopted for the solution of the incompressible Navier-Stokes equations coupled with energy equation. The proposed algorithm is verified by solving the benchmark problem of conjugate heat transfer inside a cavity having a centered body. Then a conjugate natural heat transfer problem inside a cavity having a heat-generating conducting body with constant heat flux is solved and the effect of the Rayleigh number on the heat transfer characteristics inside a cavity is investigated.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.575-589
• /
• 2008

In the present work, we investigate the hydrodynamic behavior of a turbulent bore, such as tsunami bore and tidal bore, generated by the removal of a gate with water impounded on one side. The bore generation system is similar to that used in a general dam-break problem. In order to the numerical simulation of the formation and propagation of a bore, we consider the incompressible flows of two immiscible fluids, liquid and gas, governed by the Navier-Stokes equations. The interface tracking between two fluids is achieved by the volume-of-fluid (VOF) technique and the M-type cubic interpolated propagation (MCIP) scheme is used to solve the Navier-Stokes equations. The MCIP method is a low diffusive and stable scheme and is generally extended the original one-dimensional CIP to higher dimensions, using a fractional step technique. Further, large eddy simulation (LES) closure scheme, a cost-effective approach to turbulence simulation, is used to predict the evolution of quantities associated with turbulence. In order to verify the applicability of the developed numerical model to the bore simulation, laboratory experiments are performed in a wave tank. Comparisons are made between the numerical results by the present model and the experimental data and good agreement is achieved.

• Journal of the Society of Naval Architects of Korea
• /
• v.32 no.3
• /
• pp.62-71
• /
• 1995

In the numerical analysis of incompressible unsteady Navier-stokes equation, large time is required for solving the pressure Poisson equation of the elliptic type at each time step. In this paper, a numerical analysis by the direct method is carried out to solve the pressure Poisson equation and the computing time is analyzed as mesh size increases. The pressure Poisson equation can be transformed to the boundary value problem by the Green theorem. The computing time for the convolution type of the domain integral can be reduced by using F.F.T. and the computing time in the direct method depends entirely on obtaining the solution of the boundary value problem. The numerical analysis on the known solutions is carried out and compared for the verification of the direct method. And the numerical analysis on the body boundary and domain decomposition problem are carried out with the computing time less than O($n^{3}$) in the (n.n) mesh.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.98-111
• /
• 1993

Two-dimensional nonlinear free-surface wave problems are analyzed with consideration of viscosity. Navier-Stokes equation and continuity equation are solved by the application of Finite Analytic Method, and MAC scheme is used far the treatment of free surface. Surface tension effect is also considered and laminar flow is assumed. The free-surface waves in shallow water, the flows around a vortex-pair with free surface and the wave ahead of a rectangular body are simulated to test the present numerical scheme. In the shallow water problem, viscous effect due to the friction on the bottom is observed. In the second problem, the approach of a vortex-pair to the free surface is simulated to examine the interaction of vortex-pair with the free surface. In the third problem, the wave ahead of a semi-infinite floating body is simulated.