• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2004.04a
• /
• pp.88-91
• /
• 2004

The multiphase flow simulator, MPS, is developed based on the fractional flow approach considering tile fully three phase flow with general initial and boundary condition. Most existing fractional flow-based models are limited to two-phase flow and specific boundary conditions. Although there appears a number of three-phase flow models, they were mostly developed using pressure based approaches. As a result, these models require cumbersome variable-switch techniques to deal with phase appearance and disappearance. The use of fractional flow based approach in MPS makes it unnecessary to use variable-switch to handle the change of phase configurations. Also most existing fractional flow based models consider only specific boundary conditions. However, the present model considers general boundary conditions of most possible and plausible cases which consists of ten cases.

• Economic and Environmental Geology
• /
• v.41 no.6
• /
• pp.731-746
• /
• 2008

Most multiphase flow simulators following fractional flow approach assume incompressibility of fluid and matrix or consider only two phase flow (water and air, water and NAPL). However, in this study, mathematical governing equations were developed for fully compressible three-phase flow using fractional flow based approach. Also, fully compressible multiphase flow simulator (CMPS) considering compressibilities of matrix and fluid was developed using the mathematical governing equations. In order to verify CMPS, the CMPS were compared with analytical solution and the existing multiphase flow simulator, MPS, which had been developed for simulating incompressible multiphase flow (Suk and Yeh 2007; Suk and Yeh 2008). According to the results, solutions of CMPS and MPS and analytical solutions are well matched each other. Thus, it is found that CMPS has the capability of simulating compressible three phase flow phenomena assuming compressibilities of fluids and matrix.

• Economic and Environmental Geology
• /
• v.44 no.2
• /
• pp.161-170
• /
• 2011

The multiphase flow simulator, CHEMPS, was developed based on the fractional flow approach reported in the petroleum engineering literature considering fully three phase flow in physically and chemically heterogeneous media. It is a extension of MPS developed by Suk and Yeh (2008) to include the effect of wettability on the migration of NAPL. The fractional flow approach employs water, total liquid saturation and total pressure as the primary variables. Most existing models are limited to two-phase flow and specific boundary conditions when considering physically heterogeneous media. In addition, these models focused mainly on the water-wet media. However, in a real system, variations in wettability between water-wet and oil-wet media often occur. Furthermore, the wetting of porous media by oil can be heterogeneous, or fractional, rather than uniform due to the heterogeneous nature of the subsurface media and the factors that affect the wettability. Therefore, in this study, the chemically heterogeneous media considering fractional wettability as well as physically heterogeneous media were simulated using CHEMPS. In addition, the general boundary conditions were considered to be a combination of two types of boundaries of individual phases, flux-type and Dirichlet type boundaries.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.33-38
• /
• 2001

Bifurcation of unstable symmetric flow patterns to stable asymmetric ones in laminar sudden-expansion flow has been numerically investigated. Computations were carried out for an expansion ratio of 3 and over a range of the flow Reynolds numbers by using numerical methods of second-order time accuracy and a fractional-step method that guarantees divergence-free flowfields at all times. The critical Reynolds number above which bifurcation of pitchfork type to asymmetric flow pattern takes place is lower in a flow with a higher expansion ratio, in agreement with the previously reported results. The bifurcation diagrams show that the bifurcation takes place at a Reynolds number, $Re_c = 86.3$, higher than the value that has been reported. The lower critical Reynolds number may be due to deficiencies in their computations which employed SIMPLE-type relaxation methods rather than the initial-value approach of the present study. Characteristics of the flow development during the transition to asymmetric stable flow have been investigated by using spectral analysis of the velocity signals obtained by the simulations.

• Journal of Mechanical Science and Technology
• /
• v.16 no.5
• /
• pp.683-695
• /
• 2002

The objective of the present study is to analyze the fluid flow with moving boundary using a finite element method. The algorithm uses a fractional step approach that can be used to solve low-speed flow with large density changes due to intense temperature gradients. The explicit Lax-Wendroff scheme is applied to nonlinear convective terms in the momentum equations to prevent checkerboard pressure oscillations. The ALE (Arbitrary Lagrangian Eulerian) method is adopted for moving grids. The numerical algorithm in the present study is validated for two-dimensional unsteady flow in a driven cavity and a natural convection problem. To extend the present numerical method to engine simulations, a piston-driven intake flow with moving boundary is also simulated. The density, temperature and axial velocity profiles are calculated for the three-dimensional unsteady piston-driven intake flow with density changes due to high inlet fluid temperatures using the present algorithm. The calculated results are in good agreement with other numerical and experimental ones.

• Wind and Structures
• /
• v.15 no.3
• /
• pp.189-208
• /
• 2012

Lock-in and drag amplification phenomena are studied for a flexible cantilever using a simplified fluid-structure interaction approach. Instead of solving the 3D domain, a simplified setup is devised, in which 2D flow problems are solved on a number of planes parallel to the wind direction and transversal to the structure. On such planes, the incompressible Navier-Stokes equations are solved to estimate the fluid action at different positions of the line-like structure. The fluid flow on each plane is coupled with the structural deformation at the corresponding position, affecting the dynamic behaviour of the system. An Arbitrary Lagrangian-Eulerian (ALE) approach is used to take in account the deformation of the domain, and a fractional-step scheme is used to solve the fluid field. The stabilization of incompressibility and convection is achieved through orthogonal quasi-static subscales, an approach that is believed to provide a first step towards turbulence modelling. In order to model the structural problem, a special one-dimensional element for thin walled cross-section beam is implemented. The standard second-order Bossak method is used for the time integration of the structural problem.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.12
• /
• pp.4003-4012
• /
• 1996

Hydro dynamically developing and simultaneously (hydro dynamically and thermally) developing laminar flows of a Bingham plastic in a circular pipe have been investigated numerically. Solutions have been obtained by using a four-step fractional method combined with an equal order bilinear finite element method. For the hydro dynamically developing flow, shorter entrance length is required to reach fully developed velocity field for larger yield stress and non-monotonic pressure drop along the pipe centerline is observed when the yield stress exceeds a certain critical value. For the simultaneously developing flow, the heat transfer characteristics show the same trends as those predicted for the thermally developing flow (Graetz problem).

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.49-55
• /
• 2007

The quality of chaotic mixing in three-dimensional micro channel flow has been numerically studied using Fractional-step method (FSM) and particle tracking techniques such as $Poincar{\acute{e}}$ section and Lyapunov exponents. The flow was driven by pressure distribution and the chaotic mixing was generated by applying alternating current to electrodes embedded on the bottom wall at a first half period and on the top wall at a second half period. The equations governing the velocity and concentration distributions were solved using FSM based on Finite Volume approach. Results showed that the mixing quality depended significantly on the modulation period. The modulation period for the best mixing performance was determined based on the mixing index for various initial conditions of concentration distribution. The optimal values of modulation period obtained by the particle tracking techniques were compared with those from the solution of concentration distribution equation using FSM and CFX software and the comparison showed their good match.

• Advances in aircraft and spacecraft science
• /
• v.8 no.4
• /
• pp.345-372
• /
• 2021

The analysis of nonlinear vibrations, buckling, post-buckling, flutter boundary determination and post-flutter behavior of a homogeneous curved plate assuming cylindrical bending is conducted in this article. Other assumptions include simply-supported boundary conditions, supersonic aerodynamic flow at the top of the plate, constant pressure conditions below the plate, non-viscous flow model (using first- and third-order piston theory), nonlinear structural model with large deformations, and application of mechanical and thermal loads on the curved plate. The analysis is performed with constant environmental indicators (flow density, heat, Reynolds number and Mach number). The material properties (i.e., coefficient of thermal expansion and modulus of elasticity) are temperature-dependent. The equations are derived using the principle of virtual displacement. Furthermore, based on the definitions of virtual work, the potential and kinetic energy of the final relations in the integral form, and the governing nonlinear differential equations are obtained after fractional integration. This problem is solved using two approaches. The frequency analysis and flutter are studied in the first approach by transferring the handle of ordinary differential equations to the state space, calculating the system Jacobin matrix and analyzing the eigenvalue to determine the instability conditions. The second approach discusses the nonlinear frequency analysis and nonlinear flutter using the semi-analytical solution of governing differential equations based on the weighted residual method. The partial differential equations are converted to ordinary differential equations, after which they are solved based on the Runge-Kutta fourth- and fifth-order methods. The comparison between the results of frequency and flutter analysis of curved plate is linearly and nonlinearly performed for the first time. The results show that the plate curvature has a profound impact on the instability boundary of the plate under supersonic aerodynamic loading. The flutter boundary decreases with growing thermal load and increases with growing curvature.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.8 no.3
• /
• pp.116-121
• /
• 2005

In the present paper, a new calculation algorithm far solving large scale environmental or geophysical flows with free surface is proposed where the non-hydrostatic pressure component is taken into consideration. Predictor-corrector fractional step approach with explicit, forward time marching scheme in the sigma coordinate system is employed. In order to validate the present calculation algorithm and to estimate the effects of non-hydrostatic pressure on resultant flow and free surface movements, example calculations are carried out for typical steady and unsteady flow problems. Present method can be applied to the meso-scale free surface flows with complex bottom topography where MAC-like 3-d hydrodynamic calculations are quite ineffective and uneconomic.