• Journal of computational fluids engineering
• /
• v.10 no.1
• /
• pp.8-15
• /
• 2005

Aerodynamic sounds generated by a uniform flow around a two-dimensional circular cylinder at Re=150 are simulated by applying the finite difference lattice Boltzmann method. Thethird-order-accurate up-wind scheme (UTOPIA) is used for the spatial derivatives, and the second-order-accurate Runge-Kutta scheme is applied for the time marching. We have succeed in capturing very small pressure fluctuations with the same frequency of the Karman vortex street compared with the pressure fluctuation around a circular cylinder. The propagation velocity of the acoustic waves shows that the points of peak pressure are biased upstream due to the Doppler effect in the uniform flow. For the downstream, on the other hand, it is faster. It is also apparent that the amplitude of sound pressure is proportional to r /sup -1/2/,r being the distance from the center of the circular cylinder. To investigate the effect of the lattice dependence, furthermore, 2D computations of the tone noises radiated by a square cylinder and NACA0012 with a blunt trailing edge at high incidence and low Reynolds number are also investigate.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.1
• /
• pp.43-57
• /
• 1999

The fluid flow and heat transfer in a thin liquid film are investigated numerically. The flow Is assumed to be two-dimensional laminar and surface tension is considered. The most important characteristics of this flow is the existence of a hydraulic jump through which the flow undergoes very sharp and discontinuous change. Arbitrary Lagrangian-Eulerian(ALE) method is used to describe moving free boundary and a modified SIMPLE algorithm based on streamline upwind Petrov-Galerkin(SUPG) finite element method is used for time marching iterative solution. The numerical results obtained by solving unsteady full Navier-Stokes equations are presented for planar and radial flows subject to constant wall temperature or constant wall heat flux, and compared with available experimental data. It Is discussed systematically how the inlet Reynolds and Froude numbers and surface tension affect the formation of a hydraulic jump. In particular, the effect of temperature dependent fluid properties is also discussed.

• Journal of Advanced Marine Engineering and Technology
• /
• v.28 no.5
• /
• pp.717-727
• /
• 2004

In this study, we simulate the aerodynamic sounds generated by a two-dimensional circular cylinder in a uniform flow are simulated by applying the finite difference lattice Boltzmann method (FDLBM). The third-order-accurate up-wind scheme (UTOPIA) is used for the spatial derivatives. and the second-order-accurate Runge-Kutta scheme is applied for the time marching. The results show that we successively capture very small acoustic pressure fluctuations with the same frequency of the Karman vortex street compared with the Pressure fluctuation around a circular cylinder The propagation velocity of the acoustic waves shows that the points of peak pressure are biased upstream due to the Doppler effect in the uniform flow For the downstream. on the other hand. it quickly Propagates. It is also apparent that the amplitude of sound Pressure is Proportional to $r^{-1/2}$, r being the distance from the center of the circular cylinder. To investigate the effect of the lattice dependence furthermore a 2D computation of the tone noise radiated by a NACA0012 with a blunt trailing edge at high incidence and low Reynolds number is also investigated.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.804-809
• /
• 2003

In this study, 2D computations of the Aeolian tones for some obstacles (circular cylinder, square cylinder and NACA0012 airfoil) are simulated. First of all, we calculate the flow noise generated by a uniform flow around a two-dimensional circular cylinder at Re=150 are simulated by applying the finite difference lattice Boltzmann method (FDLBM). The third-order-accurate up-wind scheme (UTOPIA) is used for the spatial derivatives, and the second-order-accurate Runge-Kutta scheme is applied for the time marching. The results show that we successively capture very small acoustic pressure fluctuation with the same frequency of the Karman vortex street compared with the pressure fluctuation around a circular cylinder. The propagation velocity of the acoustic waves shows that the points of peak pressure are biased upstream due to the Doppler effect in the uniform flow. For the downstream, on the other hand, it is faster. To investigate the effect of the lattice dependence, furthermore, simulations of the Aeolian tones at the low Reynolds number radiated by a square cylinder and a NACA0012 airfoil with a blunt trailing edge at high incidence are also investigated.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.4
• /
• pp.48-59
• /
• 1996

The numerical step in the unsteady viscous flow analysis can be divided in the space analysis step satisfying continuity equation and the time marching step. In this study the spectral method is applied to solve the pressure Poisson equation in the space analysis step. If the highest order differential term of the pressure Poisson equation is transformed by Fourier series, pressure arid its first derivatives can be expressed by the integral form of Fourier series. So Gibb's phenomena can be eliminated and the spectral method can be applied to non-periodic problems. The numerical analysis of unsteady viscous flow around 2-dimensional circular cylinder and wing is carried out and compared for verification.

• Journal of Mechanical Science and Technology
• /
• v.20 no.2
• /
• pp.286-294
• /
• 2006

Flows in a ramjet inlet is simulated for the study of the rocket-ramjet transition. The flow is unsteady, two-dimensional axisymmetric, compressible and turbulent. Double time marching method is used for the unsteady calculation and HLLC method is used as a higher order MUSCL method. As for turbulent calculation, $\kappa-\omega$ SST model is used for more accurate viscous calculations. Sinusoidal pressure perturbation is given at the exit and the flow fields at the inlet is studied. The cruise condition as well as the ground test condition are considered. The pressure level for the ground test condition is relatively low and the effect of the pressure perturbation at the combustion chamber is small. The normal shock at the cruise condition is very sensitive to the pressure perturbation and can be easily detached from the cowl when the exit pressure is relatively high. The sudden decrease in the mass flux is observed when the inlet flow becomes subcritical, which can make the inlet incapable. The amplitude of travelling pressure waves becomes larger as the downstream pressure increases, and the wavelength becomes shorter as Mach number increases. The phase difference of the travelling perturbed pressure wave in space is 180 degree.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.6
• /
• pp.1408-1416
• /
• 1990

A powerful and efficient method for finding approximate solutions to initial-boundary-value problems in the transient coupled thermoelasticity is formulated in time domain using the finite element technique with time-marching strategy. The final system equations can be derived by the Guritin's variational principle using the definition of convolution integral. But, the finite element formulation for the equations of motion is modified by differentiating in time. Numerical results to some test problems are compared with analytical and other sophisticated approximate solutions. Stable responces are observed in all the given examples irrespective of incremental time steps and mesh shapes. In addition, it is shown that good numerical results are obtained even in coarser mesh or larger time step comparing to other numerical methods.

• Fire Science and Engineering
• /
• v.11 no.4
• /
• pp.3-14
• /
• 1997

The smoke filling process for the atrium space containing a fire source is simulated using two types of deterministic fire model : Zone model and Field model. The zone model used is the CFAST(version 1.6) model developed at the Building and Fire Research Laboratories, NIST in the USA. The field model is a self-developed frie field model based on Computational Fluid Dynamic (CFD) theories. This article is focused on finding out the smoke movement and temperature distribution in atrium space which is cubic in shape. For solving the liked set of velocity and pressure equation, the PISO algorithm, which strengthened the velocity-pressure coupling, was used. Since PISO algorithm is a time-marching procedure, computing time si very fast. A computational procedure for predicting velocity and temperature distribution in fire-induced flow is based on the solution, in finite volume method and non-staggered grid system, of 3-dimensional equations for the conservation of mass, momentum, energy, species and so forth. The fire model i.e Zone model and Field model predicted similar results for clear heights and the smoke layer temperature.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.21 no.4
• /
• pp.375-382
• /
• 2008

The particular integral formulation for two(2D) and three(3D) dimensional inelastic transient dynamic stress analysis is presented. The elastostatic equation is used for the complementary solution. Using the concept of global shape function, the particular integrals for displacement and traction rates are obtained to approximate acceleration of the inhomogeneous equation. The Houbolt time integration scheme is used for the time-marching process. The Newton-Raphson algorithm for plastic multiplier is used to solve the system equation. Numerical results of four example problems are given to demonstrate the validity and accuracy of the present formulation.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.6
• /
• pp.703-714
• /
• 2010

This paper proposes an approach for staged finite element modeling with coupled seepage and stress analysis. The stage modeling is based on the predefined inter-relationship between the base model and the unit stage models. A unit stage constitutes a complete finite element model, of which the geometries and attributes are subject to changes from stage to stage. The seepage analysis precedes the mechanical stress analysis at every stage. Division of the wet and dry zone and the pore pressures are evaluated from the seepage analysis and used in determining input data for the stress analysis. The results of the stress analysis may also be associated with the pore water pressures. For consolidation analysis, the pore pressure and the displacement variables are mixed in a coupled matrix equation. The time marching solution produces the dissipation of excess pore pressure and variation of stresses with passage of time. For undrained analysis, the excess pore pressures are computed from the stress increment due to loading applied in the unit stage and are used in revising the hydraulic head. The solution results of a unit stage are inherited and accumulated to the subsequent stages through the relationship of the base model and the individual unit stages. Implementation of the proposed approach is outlined on the basis of the core procedures, and numerical examples are presented for demonstration of its application.