• 한국연소학회:학술대회논문집
• /
• 2004.06a
• /
• pp.95-101
• /
• 2004

Dynamic behavior of diffusive-thermal instability in diluted $CH_4/O_2$ diffusion flames is numerically investigated by adopting detailed chemistry and transport. Counterflow diffusion flame is adopted as a model flamelet. Particular attention is focused on the pulsating-instability regime, which arises for Lewis numbers greater than unity, and the instability occurs at high strain rate near extinction condition in this flame configuration. Once a steady flame structure is obtained for a prescribed value of initial strain rate. transient solution of the flame is calculated after a finite amount of strain-rate perturbation is imposed Oil the steady flame. Transient evolution of the flame depends on the initial strain rate and the amount of perturbed strain rate. Basically, the dynamic behaviors can be classified into two types, namely non-oscillatory decaying solution and diverging solution leading to extinction. The peculiar oscillatory solution. which has been found in the previous study adopting one-step chemistry and constant Lewis numbers, is not observed in this study, which is attributed to both convective flow and preferential diffusion effects.

• Journal of the Society of Naval Architects of Korea
• /
• v.48 no.6
• /
• pp.581-591
• /
• 2011

This paper describes the air compressibility effect in the CFD simulation of water impact load prediction. In order to consider the air compressibility effect, two sets of governing equations are employed, namely the incompressible Navier-stokes equations and compressible Navier-Stokes equations that describe general compressible gas flow. In order to describe violent motion of free surface, volume-of-fluid method is utilized. The role of air compressibility is presented by the comparative study of water impact load obtained from two different air models, i.e. the compressible and incompressible air. For both cases, water is considered as incompressible media. Compressible air model shows oscillatory behavior of pressure on the solid surface that may attribute to the air-cushion effect. Incompressible air model showed no such oscillatory behavior in the pressure history. This study also showed that the CFD simulation can capture the formation of air pockets enclosed by water and solid surface, which may be the location where the air compressibility effect is dominant.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.3
• /
• pp.243-259
• /
• 2016

The Richards equation for water movement in unsaturated soil is highly nonlinear partial differential equations which are not solvable analytically unless unrealistic and oversimplifying assumptions are made regarding the attributes, dynamics, and properties of the physical systems. Therefore, conventionally, numerical solutions are the only feasible procedures to model flow in partially saturated porous media. The standard Finite element numerical technique is usually coupled with an Euler time discretizations scheme. Except for the fully explicit forward method, any other Euler time-marching algorithm generates nonlinear algebraic equations which should be solved using iterative procedures such as Newton and Picard iterations. In this study, lumped mass and distributed mass in the frame of Picard and Newton iterative techniques were evaluated to determine the most efficient method to solve the Richards equation with finite element model. The accuracy and computational efficiency of the scheme and of the Picard and Newton models are assessed for three test problems simulating one-dimensional flow processes in unsaturated porous media. Results demonstrated that, the conventional mass distributed finite element method suffers from numerical oscillations at the wetting front, especially for very dry initial conditions. Even though small mesh sizes are applied for all the test problems, it is shown that the traditional mass-distributed scheme can still generate an incorrect response due to the highly nonlinear properties of water flow in unsaturated soil and cause numerical oscillation. On the other hand, non oscillatory solutions are obtained and non-physics solutions for these problems are evaded by using the mass-lumped finite element method.

• Journal of computational fluids engineering
• /
• v.9 no.4
• /
• pp.34-40
• /
• 2004

It is known that the non-planar model of bypass is more profitable to suppress the development of intimal hyperplasia that tends to occur preferentially in regions of low time averaged shear stress and rapid temporal changes in wall shear stress. In this study it was numerically simulated the blood flow in an coronary artery grafted by artificial bypass to determine the flow characteristic variations due to the anastomosis angle changing. 5 different non-planar anastomosis angle models such as 45°, 60°, 90°, 120° and 135° were considered. When the anastomosis angle is higher, the backward flow region is spatially extended near the downstream region of the anastomosis because of the development of horseshoes vortex. For the case of the nan-planar 45° and 60° of anastomosis, the area of low-OSI zone was decreased by 26% and 13% respectively and the time averaged wall shear stress was increased by more than 55% as compared with 45° of planar model. However, both of the area of the low-OSI zone and the time averaged wall shear stress of 90°, 120° model were significantly increased.

• Journal of the korean Society of Automotive Engineers
• /
• v.15 no.2
• /
• pp.92-104
• /
• 1993

In the present study, the velocity profiles and entrance length of developing transitional pulsating flows are investigated both analytically and experimentally in the entrance region of a square duct. The systems of conservation equations for transitional pulsating flows in a square duct are solved analytically by linearizing the non-linear convective terms. Analytical solutions are obtained in the form of infinite series for velocity pofiles. The experimental study for the air flow in a square duct(40mm*40mm*4000mm) is carried out to measure velocity profiles and other parameters by using a hot-wire anemometer with a data acquisition and processing system. The distribution of velocity profiles( $u_{ps}$ / $u_{m,ta}$) in the decelerating period is higher than in the accelerating period. The distribution of the axial component of the axial component of velocity in the transitional flow is nearly uniform in the central region of the duct, and decrease rapidly near the wall. The entrance length correlation of the transitional pulsating flows in a square duct is obtained to be $L_{e}$/ $D_{h}$=0.83 $A_{1}$R $e_{ta}$ /(.omega. sup+1)$^{2}$TEX>

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.1 s.232
• /
• pp.103-109
• /
• 2005

Airflow in the nasal cavity of a normal Korean adult is investigated experimentally by tomographic PIV measurement. Knowledge of airflow characteristics in nasal cavities is essential to understand the physiology and pathology aspects of nasal breathing. Several studies have utilized physical models of the healthy nasal cavity to investigate the relationship between nasal anatomy and airflow. All of these researches on nasal airflow are under the condition of constant flow-rate. In this study, nasal cavity flow with the physiological period is investigated by tomographic PIV, for the first time. A pumping system that can produce the periodic flow is created. Thanks to a new method for the model casting by a combination of the rapid prototyping and curing of clear silicone, a transparent rectangular box containing the complex nasal cavity can be made for PIV, The CBC PIV algorithm is used for analysis. Phase-averaged mean and RMS velocity distributions are obtained for inspirational and expiration nasal airflows. The comparison with the constant flow case is appreciated. There exist many flow patterns depending on each phase.

• Journal of Biomedical Engineering Research
• /
• v.36 no.6
• /
• pp.264-270
• /
• 2015

A mock circulatory loop system has been developed to construct a simulator for trainees in cardiopulmonary bypass systems or to simulate a test environment for cardiac-assist devices. This paper proposes a computerized mock circulatory loop system whose node is modularized by using a servo control flow regulator to simulate dynamic change of the hemodynamic status. To observe the effect of time-varying resistance, one with hemodynamic properties, the proposed system replicates the planned cross-sectional areas of the outlet of a ventricular assist device in terms of voltage input of a servo valve. The experiment is performed (1) for steady-input commands of selected area sizes and (2) for dynamic commands such as monotonous increase and decrease, and oscillatory functions of the voltage input, and a computer program based on LabVIEW (National Instruments, Austin, USA) processes every measured data and control command to the servo valve. The results show that the pressure and flow at the target points with respect to time-varying resistance match intuitive estimation: the pressure at the outlet and the pressure drop between both sides of the valve increased and the flow at the outlet decreased for increased resistance.

• Journal of the Korean Society of Visualization
• /
• v.20 no.1
• /
• pp.52-63
• /
• 2022

Direct numerical simulation of blood flow in a stenosed, patient-specific carotid artery was conducted to explore the transient behavior of blood flow with special emphasis on the wall-shear stress distribution over the transition region. We assumed the blood as an incompressible Newtonian fluid, and the vessel was treated as a solid wall. The pulsatile boundary condition was applied at the inlet of the carotid. The Reynolds number is 884 based on the inlet diameter, and the maximum flow rate and the corresponding Womersley number is approximately 5.9. We found the transitional behavior during the acceleration and deceleration phases. In order to quantitatively examine the wall-shear stress distribution over the transition region, the probability density function of the wall-shear stress was computed. It showed that the negative wall-shear stress events frequently occur near peak systole. In addition, the oscillatory shear stress index was used to further analyze the relationship with the negative wall-shear stress appearing in the systolic phase.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.18 no.10
• /
• pp.1057-1064
• /
• 2008

In this paper, static and oscillatory instability of nanopipes conveying fluid and modelled as a thin-walled beam is investigated. Effects of boundary conditions and non-classical transverse shear and rotary inertia are incorporated in this study. The governing equations and the three different boundary conditions are derived through Hamilton's principle. Numerical analysis is performed by using extend Galerkin method which enables us to obtain more exact solutions compared with conventional Galerkin method. Variations of critical flow velocity for different boundary conditions of carbon nanopipes are investigated and pertinent conclusion is outlined.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.9
• /
• pp.1210-1218
• /
• 2000

Thermal convection in a horizontal annulus is considered, and the bifurcation phenomena of flows from time-periodic to chaotic convection are numerically investigated. The unsteady two-dimensional streamfunction-vorticity equation is solved with finite difference method. As Rayleigh number is increased, the steady flow bifurcates to a time-periodic flow with a fundamental frequency, and afterwards a period-tripling bifurcation occurs with further increase of the Rayleigh number. Chaotic convection is established after a period-doubling bifurcation. A periodic convection with period 4 appears after the first chaotic convection. At still higher Rayleigh numbers, chaotic flows reappear.