• Journal of computational fluids engineering
• /
• v.16 no.4
• /
• pp.7-13
• /
• 2011

A comparative analysis of thermal models in the lattice Boltzmann method(LBM) for the simulation of laminar natural convection in a square cavity is presented. A HYBRID method, in which the thermal equation is solved by the Navier-Stokes equation method while the mass and momentum conservation are resolved by the lattice Boltzmann method, is introduced and its merits are explained. All the governing equations are discretized on a cell-centered, non-uniform grid using the finite-volume method. The convection terms are treated by a second-order central-difference scheme with a deferred correction method to ensure stability of the solutions. The HYBRID method and the double-population method are applied to the simulation of natural convection in a square cavity and the predicted results are compared with the benchmark solutions given in the literatures. The predicted results are also compared with those by the conventional Navier-Stokes equation method. In general, the present HYBRID method is as accurate as the Navier-Stokes equation method and the double-population method. The HYBRID method shows better convergence and stability than the double-population method. These observations indicate that this HYBRID method is an efficient and economic method for the simulation of incompressible fluid flow and heat transfer problem with the LBM.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.88-97
• /
• 1996

This paper treats an adaptive finite-element method for the viscous compressible flow governed by Navier-Stokes equations in two dimensions. The numerical algorithm is the two-step Taylor-Galerkin mettled using unstructured triangular grids. To increase accuracy and stability, combined moving node method and grid refinement method have been used for grid adaption. Validation of the present algorithm has been made by comparing the present computational results with the existing experimental data and other numerical solutions. Four benchmark problems are solved for demonstration of the present numerical approach. They include a subsonic flow over a flat plate, the Carter flat plate problem, a laminar shock-boundary layer interaction. and finally a laminar flow around NACA0012 airfoil at zero angle of attack and free stream Mach number of 0.85. The results indicates that the present adaptive triangular grid method is accurate and useful for laminar viscous flow calculations.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.6B
• /
• pp.379-385
• /
• 2012

Recent development of computing power and theoretical advances in computational fluid dynamics have made possible numerical simulations of water waves with full Navier-Stokes equations. In this study, an internal wave maker using the mass source function approach was combined with the level set finite element method for generation of waves. The model is first applied to the two-dimensional linear wave generation and propagation. Then, it is applied to the three-dimensional simulation of the same problem. To effectively utilize computational resources and enhance the speed of execution, parallel algorithms are developed and applied for the three-dimensional problem. The results of numerical simulations are compared with theoretical values and good agreements are observed.

• Journal of the Korean Society of Industry Convergence
• /
• v.5 no.1
• /
• pp.45-52
• /
• 2002

In this study, an exact solution of governing differential equation for the bending problem of partially loaded orthotropic rectangular plates is presented and also its computer program is developed. The method requires that two opposite edges be clamped or simply supported, or one edge clamped and the other simply supported. Any combination of boundary conditions could exist along the other edges. The plate could he subjected to uniform, partially uniform, and line loads. The solution for the deflection of rectangular plate is expressed as a Levy type single Fourier series and the loads arc expressed as a corresponding series. The advantage of the solution is that it overcomes the limitations of the previous Navier's and Levy's methods (limitation of boundary condition and loading conditions of plate), it is easy to program on a computer and it becomes fast to solve the bending problem with computer program. Calculations are presented for isotropic and orthotropic plates with different loading and boundary conditions. Comparisons are made for the isotropic plate with various boundary conditions between the result of this paper and the result of Navier, Levy and Szilard. The deflections were in excellent agreement.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.8
• /
• pp.661-669
• /
• 2007

Study on the feasibility assessment for applying the overlap grid method to numerical calculations on the ground effect has been performed. The objective of the present study is to settle the problem in the grid generation process. A low Mach number preconditioned turbulent flow solver using the overlap grid and the multi-block grid methods has been developed and applied to the ground effect simulation around the RAE 101 airfoil. It has been verified that the overlap grid method not only can provide sufficiently accurate solutions but also work out the grid generation problem in the ground effect simulations.

• Structural Engineering and Mechanics
• /
• v.81 no.4
• /
• pp.523-527
• /
• 2022

The objective of this study is to simulate the flow in pipes with various boundary conditions. Free-pressure fluid model, is used in the pipe based on Navier-Stokes equation. The models are solved by using the numerical method. A problem called "stability of pipes" is used in order to compare frequency and critical fluid velocity. When the initial conditions of problem satisfied the instability conditions, the free-pressure model could accurately predict discontinuities in the solution field. Employing nonlinear strains-displacements, stress-strain energy method the governing equations were derived using Hamilton's principal. Differential quadrature method (DQM) is used for obtaining the frequency and critical fluid velocity. The results of this paper are analyzed by hyperbolic numerical method. Results show that the level of numerical diffusion in the solution field and the range of well-posedness are two important criteria for selecting the two-fluid models. The solutions for predicting the flow variables is approximately equal to the two-pressure model 2. Therefore, the predicted pressure changes profile in the two-pressure model is more consistent with actual physics. Therefore, in numerical modeling of gas-liquid two-phase flows in the vertical pipe, the present model can be applied.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.13 no.4
• /
• pp.19-24
• /
• 1976

Some method of analysis of rectangular plates under distributed load of various intensity with all edges built in are presented in. Analysis of many structures such as bottom, side shell, and deck plate of ship hull, and flat slab, deck systems of bridges is a problem of plate with continuous supports or clamped edges. When the four edges of rectangular plate is simply supported, the double fourier series solution developed by Navier can represent an exact result of this problem. If two opposite edges are simply supported, Levy's method is available to give an "exact" solution. When the loading condition and boundary condition of a plate does not fall into these cases, no simple analytic method seems to be feasible. Analysis of a plate under distributed loads of various intensity with all edges built in is carried out by applying Navier solution and Levy's method as well as "Principle of Superposition" In discussing this problem we start with the solution of the problem for a simply supported rectangular plate and superpose on the deflection of such a plate the deflections of the plate by slopes distributed along the all edges. These slopes we adjust in such a manner as to satisfy the condition of no rotation at the boundary of the clamped plate. This method can be applied for the cases of plates under irregularly distributed loads of various intensity with two opposite edges simply supported and the other two edges clamped and all edges simply supported and this method can also be used to solve the influence values of deflection, moment and etc. at arbitrary position of plates under the live load.

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

We present a numerical simulation technique and some preliminary results of the impact and spreading of a droplet containing particles on the solid substrate in 2D. We used the 2nd-order Adams-Bashforth / Crank-Nicholson method to solve the Navier-Stokes equation and employed the level-set method with the continuous surface stress for description of droplet spreading with interfacial tension. The impact velocity has been generated by the instantaneous gravity. The distributed Lagrangian-multipliers method has been combined for the implicit treatment of rigid particles and the discontinuous Galerkin method has been used for the stabilization of the interface advection equation. We investigated the droplet spreading by the inertial force and discussed effects of the presence of particles on the spreading behavior using an example problem. We observed reduced oscillation and spread for the particulate droplet.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.75-80
• /
• 2005

The present work was undertaken to numerically analyze the defrosting phenomena of windshield glass. In order to analysis the phase change from frost to water on windshield glass by discharging hot air from a defroster nozzle, the flow and the temperature field of the cabin interior, the heat transfer through the windshield glass, and the phase change of frost should be solve simultaneously. In the present work, the flow field was obtained by solving 3-D incompressible Navier-Stokes equations, and the temperature field was computed from the incompressible energy equation. The phase change process was solved by the enthalpy method. For the code validation, the temperature and the phase change of the driven cavity were calculated. The calculation showed a good agreement with other numerical results. Then, the present code was applied to the defrosting problem of a real automobile, and a good agreement with the experimental data was also obtained.

• Journal of Mechanical Science and Technology
• /
• v.18 no.4
• /
• pp.664-670
• /
• 2004

Navier-Stokes computation based on a new simplified model is proposed to investigate the interactions of a moving shock wave with multiple vortices arranged in the serial manner. This model problem simulates shock-vortexlet interactions at the shear layer of a compressible vortex often observed in the experiment. Applying the Foppl's idea, we extended the Rankin's model generally used for the description of a single vortex to the multi-vortex version. The acoustic pulses accelerated and decelerated are successively generated and propagated from each shock-vortex interaction, which simply explains the genesis of eccentrically diverging acoustic waves appearing in the experimental photograph.