• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.450-455
• /
• 2001

The time-dependent behavior of nonequilibrium condensation of moist air through the Ludwieg tube is investigated with a computational fluid dynamics(CFD) method. The two-dimensional, compressible, Navier-Stokes equations, fully coupled with the condensate droplet growth equations, are numerically solved by a third-order MUSCL type TVD finite-difference scheme, with a second-order fractional time step. Baldwin-Lomax turbulence model is employed to close the governing equations. The computational results are compared with the previous experiments using the Ludwieg tube with a downstream diaphragm. The results clearly show that for an initial relative humidity below 30% there is no periodic oscillation of the condensation shock wave, but for an initial relative humidity over 40% the periodic excursions of the condensation shock occurs in the Ludwieg tube, and the frequency increases with the initial relative humidity. It is also found that total pressure loss due to nonequilibrium condensation in the Ludwieg tube should not be ignored even for a very low initial relative humidity, and the periodic excursions of the condensation shock wave are responsible for the total pressure loss.

• Journal of Mechanical Science and Technology
• /
• v.17 no.12
• /
• pp.2066-2077
• /
• 2003

The time-dependent behavior of unsteady condensation of moist air through the Ludwieg tube is investigated by using a computational fluid dynamics (CFD) work. The two-dimensional, compressible, Navier-Stokes equations, fully coupled with the condensate droplet growth equations, are numerically solved by a third-order MUSCL type TVD finite-difference scheme, with a second-order fractional time step. Baldwin-Lomax turbulence model is employed to close the governing equations. The predicted results are compared with the previous experiments using the Ludwieg tube with a diaphragm downstream. The present computations represent the experimental flows well. The time-dependent unsteady condensation characteristics are discussed based upon the present predicted results. The results obtained clearly show that for an initial relative humidity below 30% there is no periodic oscillation of the condensation shock wave, but for an initial relative humidity over 40% the periodic excursions of the condensation shock occurs in the Ludwieg tube, and the frequency increases with the initial relative humidity. It is also found that total pressure loss due to unsteady condensation in the Ludwieg tube should not be ignored even for a very low initial relative humidity and it results from the periodic excursions of the condensation shock wave.

• Journal of KSNVE
• /
• v.4 no.2
• /
• pp.163-175
• /
• 1994

An analytical method for nautral frequencies of a partially liquid- filled circular cylindrical shell with various boundary conditions is developed by means of the Stokes's transformation and Fourier series expansion on the basis of Sanders' shell equation. The liquid-shell coupled system is divided into two regions for convenient formulation. One is the empty shell region in which the Sanders' shell equations are formulated without the lipuid effect, the other is wetted shell region in which the shell equations are formulated with consideration of the liquid dynamic effect. The shell equations for each regions are combined by the geometry and the force continuities at the junction of the two regions. For the vibration relevant to the liquid motion, the velocity potential of liquid is assumed as a sum of linear combination of suitable harmonic functions in axial direction. The unknown parameters are selected to satisfy the boundary condition along the wetted shell surface. The natural frequencies of the liquid filled cylindraical shells with the clamped- free and the clamped-clamped boundary conditions examined in the previous works, are obtained by this analytical method. The results are compared with the previous works, and excllent agreement is found for the natural frequencies of the shells.

• Coupled systems mechanics
• /
• v.10 no.2
• /
• pp.123-142
• /
• 2021

The paper deals with the study of the dispersion of quasi-Lamb waves in a hydro-elastic system consisting of an elastic plate, barotropic compressible inviscid fluid, and rigid wall. The motion of the plate is described using the exact equations of elastodynamics, however, the flow of the fluid using the linearized equations and relations of the Navier-Stokes equations. The corresponding dispersion equation is obtained and this equation is solved numerically, as a result of which the corresponding dispersion curves are constructed. The main attention is focused on the effect of the presence of the fluid and the effect of the fluid layer thickness (i.e., the fluid depth) on the dispersion curves. The influence of the problem parameters on the dispersion curves related to the quasi-Scholte wave is also considered. As a result of the analyses of the numerical results, concrete conclusions are made about the influence of the fluid depth, the rigid wall restriction on the fluid motion, and the material properties of the constituents on the dispersion curves. During the analyses, the zeroth and the first four modes of the propagating waves are considered.

• Journal of Advanced Marine Engineering and Technology
• /
• v.37 no.6
• /
• pp.617-622
• /
• 2013

The sloshing of liquid inside an ocean floater is caused by disturbances due to waves. For the analysis of sloshing impact within the floater and that of waves on the floater, the coupled analysis method is used. The Stokes $5^{th}$ order non-linear wave theory equations were adapted for wave making. Furthermore, Navier-Stokes equation and Shear-Stress Transport (SST) turbulent model were used to Computational Fluid dynamics, where the ocean floater motions are considered the heave and the pitch motion. The results obtained confirms the mutual relationship between the rigid body motions and that of sloshing, where the sloshing behaviour within the floater is characterized by the wave effects on the floater.

• 한국연소학회:학술대회논문집
• /
• 2002.06a
• /
• pp.123-132
• /
• 2002

A numerical study was conducted to investigate the combustion phenomena of normal start and unstart processes based on ISL's RAMAC 30 experiments with different diluent amounts and fill pressures in a ram accelerator. The initial projectile launching speed was 1.8 km/s which corresponded to the superdetonative speed of the stoichiometric $H_2/O_2$ mixture diluted with 5 $CO_2$ or 4 $CO_2$. Experiments with same condition except for projectile surface material demonstrated that ignition was successful with an aluminum projectile, but no combustion was observed in case of a steel projectile. In this study, it was found that neither shock nor viscous heating was sufficient to ignite the mixture at a low speed of 1.8 km/s, as was found in the experiments using a steel projectile. However, we could succeed in igniting the mixtures by imposing a minimal amount of additional heat to the combustor section and simulate the normal start and unstart processes found in the experiments with an aluminum projectile. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations coupled with a Baldwin-Lomax turbulence model and detailed chemistry reaction equations of $H_2/O_2/CO_2$ suitable for high-pressure gaseous combustion were considered. The governing equations were discretized by a high order accurate upwind scheme and solved in a fully coupled manner with a fully implicit, time accurate integration method. The numerical results matched almost exactly to the experimental results. As a result, it was found that the normal start and unstart processes depended on the strength of gas mixture, development of shock-induced combustion wave stabilized by the first separation bubble, and its size and location.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.05a
• /
• pp.247-250
• /
• 2007

Flexible media such as the paper, the film, etc. are thin, light and very flexible. They behave in geometrically nonlinear. Any of small force makes large deformation. So we must including aerodynamic effect when its behavior is predicted. Thus, it becomes fully coupled fluid-structure interaction(FSI) problem. In FSI problems, where the fluid mesh near the structure undergoes large deformations and becomes unacceptably distorted, which drive the time step to a very small value for explicit calculations, the arbitrary Lagrangian-Eulerian(ALE) methods or rezoning are used to create a new undistorted mesh for the fluid domain, which allows the calculations to continue. In this paper, FE sheet model considering geometric nonlinearity is formulated to simulate the behavior of the flexible media. Aerodynamic force to the media by surrounding air is calculated by solving the incompressible Navier-Stokes equations. Q2Q1(Taylor-Hood) element which means biquadratic for velocity and bilinear for pressure is used for fluid domain. Q2Q1 element satisfies LBB condition and any stabilization technique is not needed. In this paper, cantilevered sheet in the viscous incompressible Navier-Stokes flow is simulated to check the mesh motion and numerical integration scheme, and then falling paper in the air is simulated and the effects of some representative parameters are investigated.

• Journal of computational fluids engineering
• /
• v.17 no.3
• /
• pp.25-32
• /
• 2012

This work was devoted to compare two different cavitation models to study the dependency of model constants. The cavitation model of Merkle et al.(2006) and Kunz et al.(2000) were used for the present computational study. The cavitation models were coupled with the incompressible unsteady Reynolds-Averaged Navier-Stokes solver to indicate the vaporization and condensation processes. For this purpose, a preconditioning method was added as the pseudo-time term to solve the unsteady stiffness problems. For the validation of the numerical simulation, the computation was performed for the cavitating flow in a converging-diverging channel. The present results show that Merkle's cavitation model is independent to the model constants, and the higher numerical accuracy over Kunz's cavitation model.

• The KSFM Journal of Fluid Machinery
• /
• v.2 no.1 s.2
• /
• pp.29-34
• /
• 1999

To develop the efficient numerical optimization method for the design of an airfoil, an evaluation of various methods coupled with two-dimensional Naviev-Stokes analysis is presented. Simplex method and Hook-Jeeves method we used as direct search methods, and steepest descent method, conjugate gradient method and DFP method are used as indirect search methods and are tested to determine the search direction. To determine the moving distance, the golden section method and cubic interpolation method are tested. The finite volume method is used to discretize two-dimensional Navier-Stokes equations, and SIMPLEC algorithm is used for a velocity-pressure correction method. For the optimal design of two-dimensional airfoil, maximum thickness, maximum ordinate of camber line and chordwise position of maximum ordinate are chosen as design variables, and the ratio of drag coefficient to lift coefficient is selected as an objective function. From the results, it is found that conjugate gradient method and cubic interpolation method are the most efficient for the determination of search direction and the moving distance, respectively.

• Fire Science and Engineering
• /
• v.28 no.1
• /
• pp.52-57
• /
• 2014

In this paper, the unsteady incompressible Navier-Stokes equations coupled with energy equation were solved to find out the optimal location of electrical heat trace for anti-freeze of water inside the pipe for fire protection. Since the conduction equation of pipe was coupled with the natural convection of water, the analysis of conjugate heat transfer was conducted. A commercial code (ANSYS-FLUENT) based on SIMPLE-type algorithm was used for investigating the unsteady flows and temperature distributions in water region. From the numerical experiments, the isotherms and the vector fields in water region were obtained. Furthermore, it was found that the lowest part of the pipe cross-section was an optimal position of electrical heat trace assuming the constant thermal expansion coefficient of water since the minimum temperature of the water with the position is higher than those with the other positions.