• 천문학회지
• /
• 제34권4호
• /
• pp.237-241
• /
• 2001

Our goal is to present a simple volume-of-fluid type interface-tracking algorithm to compressible two-phase flow in two space dimensions. The algorithm uses a uniform underlying Cartesian grid with some cells cut by the tracked interfaces into two subcells. A volume-moving procedure that consists of two basic steps: (1) the update of volume fractions in each grid cell at the end of the time step, and (2) the reconstruction of interfaces from discrete set of volume fractions, is employed to follow the dynamical behavior of the interface motion. As in the previous work with a surface-tracking procedure for general front tracking (LeVeque & Shyue 1995, 1996), a high resolution finite volume method is then applied on the resulting slightly nonuniform grid to update all the cell values, while the stability of the method is maintained by using a large time step wave propagation approach even in the presence of small cells and the use of a time step with respect to the uniform grid cells. A sample preliminary numerical result for an underwater explosion problem is shown to demonstrate the feasibility of the algorithm for practical problems.

• 에너지공학
• /
• 제12권1호
• /
• pp.29-34
• /
• 2003

압력파의 일종인 열음향파는 압축성유체를 급속히 가열 또는 냉각하는 경계면 근처에서 유체가 순간적으로 압축 및 팽창하는 경우에 발생하는 현상으로 자연대류가 일어나지 않는 우주공간에서는 매우 중요한 열전달 메커니즘이다. 본 연구에서는 공기로 채워진 공간에서 급속한 가열에 의하여 발생한 열음향파의 전달특성을 수치적인 방법에 의하여 평가하고자 유한체적법을 기반으로 비정상 지배방정식을 이산화하였으며, PISO알고리즘과 2계 상향기법을 적용하여 해석을 수행하였다. 안정적인 수치해는 50 $\times$ 800 개의 셀과 1 $\times$ $10^{-9}$ 시간간격을 적용하여 얻을 수 있었으며, 생성된 열음향파는 유체 속을 통과하면서 점성과 열소산에 의하여 점점 감쇠하여 가는 경향을 보였다. 생성된 압력파는 날카로운 전단과 점점 감소하여 길게 늘어지는 후단부를 갖는 형상을 보였다.

• 대한기계학회논문집A
• /
• 제27권4호
• /
• pp.567-576
• /
• 2003

In this study, a robust controller to control pressure in a pneumatic pressure vessel with a long transmission line is proposed. Frequency response of transmission line using compressible fluid is changed by the flowing state of the fluid. So, it a fixed gain controller designed based on a model supposed the flowing state to a specific state, the performance of the control system could be degraded because of the modelling error. The controller designed in this study is composed of two parts. One is a feedback controller to improve a feedback characteristics and to compensate the influence of the variation of transfer characteristics of a transmission line owing to the change of flowing state and the other is a feedforward controller to regulate command fallowing performance. The experimental results with the designed controller show that the robustness of the control system is achieved regardless of the change of the model or the transmission line. Therefore, the designed controller can be utilized for the Performance improvement of a Pressure control system with a long transmission line using compressible fluid.

• 설비공학논문집
• /
• 제14권11호
• /
• pp.914-922
• /
• 2002

Thermo-acoustic waves can be thermally generated in a compressible flow field by rapid heating and cooling, and chemical reaction near the boundary walls. This mechanism is very important in the space environment in which natural convection does not exist. Also this may be a significant factor for heat transfer when the fluids are close to the thermodynamic critical point. In this study, the generation and transmission characteristics of thermo-acoustic waves in an air-filled confined domain with two-step pulsed heating are studied numerically. The governing equations are discretized using control volume method, and are solved using PISO algorithm and second-order upwind scheme. For the purpose of stable solution, time step was set to the order of $1\times10_-9s,\;and\;grids\;are\;50\times2000$. Results show that temperature and pressure distributions of fluid near the boundary wall subjected to a rapid heating are increased abruptly, and the induced thermo-acoustic wave propagates through the fluid until it decays due to viscous and heat dissipation. Pressure waves have sharp front shape and decay with a long tail in the case of step heating, but these waves have sharp pin shape in the case of pulsed heating.

• Structural Engineering and Mechanics
• /
• 제22권6호
• /
• pp.719-739
• /
• 2006

An analytical method for the free vibration of two circular plates coupled with an inviscid and compressible fluid is developed by the Rayleigh-Ritz method. The fluid is bounded by a rigid cylindrical vessel and two circular plates with an unequal thickness and diameter. It was found that the theoretical results could predict well the fluid-coupled natural frequencies with an excellent accuracy when compared with the finite element analysis results. As the fluid thickness increases or the plate thickness difference increases, an abrupt curve veering in the natural frequency loci of the neighboring modes and drastic changes in the corresponding mode shapes are observed. The mode localization frequently appears in the higher modes and in the wide gap between the plates because of a decrease in the fluid coupling owing to the fluid dispersion effect.

• Structural Engineering and Mechanics
• /
• 제20권4호
• /
• pp.389-403
• /
• 2005

In the present paper it is aimed to perform the stochastic dynamic analysis of fluid and fluidstructure systems by using the Lagrangian approach. For that reason, variable-number-nodes twodimensional isoparametric fluid finite elements are programmed in Fortran language by the authors and incorporated into a general-purpose computer program for stochastic dynamic analysis of structure systems, STOCAL. Formulation of the fluid elements includes the effects of compressible wave propagation and surface sloshing motion. For numerical example a rigid fluid tank and a dam-reservoir interaction system are selected and modeled by finite element method. Results obtained from the modal analysis are compared with the results of the analytical and numerical solutions. The Pacoima Dam record S16E component recorded during the San Fernando Earthquake in 1971 is used as a ground motion. The mean of maximum values of displacements and hydrodynamic pressures are compared with the deterministic analysis results.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제17권3호
• /
• pp.209-219
• /
• 2013

The present paper deals with the extension of AUSMPW+ scheme into two-fluid model for multiphase flow. AUSMPW+ scheme is the improvement of a single-phase AUSM+ scheme by designing pressure-based weighting functions to prevent oscillations near a wall and shock instability after a strong shock. Recently, Kitamura and Liou assessed a family of AUSM-type schemes with two-fluid model governing equations [K. Kitamura and M.-S. Liou, Comparative study of AUSM-Family schemes in compressible multi-phase flow simulations, ICCFD7-3702 (2012)]. It was observed that the direct application of the single-phase AUSMPW+ did not provide satisfactory results for most of numerical test cases, which motivates the current study. It turns out that, by designing pressure-based weighting functions, which play a key role in controlling numerical diffusion for two-fluid model, problems reported in can be overcome. Various numerical experiments validate the proposed modification of AUSMPW+ scheme is accurate and robust to solve multiphase flow within the framework of two-fluid model.

• Journal of applied mathematics & informatics
• /
• 제25권1_2호
• /
• pp.51-66
• /
• 2007

A layer of compressible, rotating, elastica-viscous fluid heated & soluted from below is considered in the presence of vertical magnetic field to include the effect of Hall currents. Dispersion relation governing the effect of viscoelasticity, salinity gradient, rotation, magnetic field and Hall currents is derived. For the case of stationary convection, the Rivlin-Erickson fluid behaves like an ordinary Newtonian fluid. The compressibility, stable solute gradient, rotation and magnetic field postpone the onset of thermosolutal instability whereas Hall currents are found to hasten the onset of thermosolutal instability in the absence of rotation. In the presence of rotation, Hall currents postpone/hasten the onset of instability depending upon the value of wavenumbers. Again, the dispersion relation is analyzed numerically & the results depicted graphically. The stable solute gradient and magnetic field (and corresponding Hall currents) introduce oscillatory modes in the system which were non-existent in their absence. The case of overstability is discussed & sufficient conditions for non-existence of overstability are derived.

• 호남수학학술지
• /
• 제36권1호
• /
• pp.11-27
• /
• 2014

In this paper, we propose closures for multi-phase flow models, which satisfy boundary conditions and conservation constraints. The models governing the evolution of the fluid mixing are derived by applying an ensemble averaging procedure to the microphysical equations characterized by distinct phases. We consider compressible multi species multi-phase flow with surface tension and transport.

• International Journal of Naval Architecture and Ocean Engineering
• /
• 제11권1호
• /
• pp.572-583
• /
• 2019

This paper aims to assess the applicability of the Runge Kutta Discontinuous Galerkin-Direct Ghost Fluid Method to the internal explosion inside a water-filled tube, which previously was studied by many researchers in separate works. Once the explosive charge located at the inner center of the water-filled tube explodes, the tube wall is subjected to an extremely high intensity fluid loading and deformed. The deformation causes a modification of the field of fluid flow in the region near the water-structure interface so that has substantial influence on the response of the structure. To connect the structure and the fluid, valid data exchanges along the interface are essential. Classical fluid structure interaction simulations usually employ a matched meshing scheme which discretizes the fluid and structure domains using a single mesh density. The computational cost of fluid structure interaction simulations is usually governed by the structure because the size of time step may be determined by the density of structure mesh. The finer mesh density, the better solution, but more expensive computational cost. To reduce such computational cost, a non-matched meshing scheme which allows for different mesh densities is employed. The coupled numerical approach of this paper has fewer difficulties in the implementation and computation, compared to gas dynamics based approach which requires complicated analytical manipulations. It can also be applied to wider compressible, inviscid fluid flow analyses often found in underwater explosion events.