• Structural Engineering and Mechanics
• /
• v.66 no.4
• /
• pp.505-513
• /
• 2018

Steam flooding is widely used in heavy oil reservoir with coupling effects among the formation temperature change, fluid flow and solid deformation. The effective stress, porosity and permeability in this process can be affected by the multi-physical coupling of thermal, hydraulic and mechanical processes (THM), resulting in a complex interaction of geomechanical effects and multiphase flow in the porous media. Quantification of the state of deformation and stress in the reservoir is therefore essential for the correct prediction of reservoir efficiency and productivity. This paper presents a coupled fluid flow, thermal and geomechanical model employing a program (MATLAB interface code), which was developed to couple conventional reservoir (ECLIPSE) and geomechanical (ABAQUS) simulators for coupled THM processes in multiphase reservoir modeling. In each simulation cycle, time dependent reservoir pressure and temperature fields obtained from three dimensional compositional reservoir models were transferred into finite element reservoir geomechanical models in ABAQUS as multi-phase flow in deforming reservoirs cannot be performed within ABAQUS and new porosity and permeability are obtained using volumetric strains for the next analysis step. Finally, the proposed approach is illustrated on a complex coupled problem related to steam flooding in an oil reservoir. The reservoir coupled study showed that permeability and porosity increase during the injection scenario and increasing rate around injection wells exceed those of other similar comparable cases. Also, during injection, the uplift occurred very fast just above the injection wells resulting in plastic deformation.

• Wind and Structures
• /
• v.20 no.3
• /
• pp.423-448
• /
• 2015

In this paper the unsteady fluid-structure interaction (FSI) problems with large structural displacement are solved by partitioned solution approaches in the arbitrary Lagrangian-Eulerian finite element framework. The incompressible Navier-Stokes equations are solved by the characteristic-based split (CBS) scheme. Both a rigid body and a geometrically nonlinear solid are considered as the structural models. The latter is solved by Newton-Raphson procedure. The equation governing the structural motion is advanced by Newmark-${\beta}$ method in time. The dynamic mesh is updated by using moving submesh approach that cooperates with the ortho-semi-torsional spring analogy method. A mass source term (MST) is introduced into the CBS scheme to satisfy geometric conservation law. Three partitioned coupling strategies are developed to take FSI into account, involving the explicit, implicit and semi-implicit schemes. The semi-implicit scheme is a mixture of the explicit and implicit coupling schemes due to the fluid projection splitting. In this scheme MST is renewed for interfacial elements. Fixed-point algorithm with Aitken's ${\Delta}^2$ method is carried out to couple different solvers within the implicit and semi-implicit schemes. Flow-induced vibrations of a bridge deck and a flexible cantilever behind an obstacle are analyzed to test the performance of the proposed methods. The overall numerical results agree well with the existing data, demonstrating the validity and applicability of the present approaches.

• Journal of Ocean Engineering and Technology
• /
• v.31 no.1
• /
• pp.14-21
• /
• 2017

The erosion of solid particles in a pipe elbow was numerically investigated. A numerical procedure to estimate the sand erosion rate, as well as the particle motion, in the pipe elbow flow was introduced. This procedure was performed based on the combined empirical erosion model and computational fluid dynamics (CFD) analysis to consider the interaction between the particle motion and the eroded surface. The underlying turbulent flow on an Eulerian frame is described by the Reynolds averaged Navier-Stokes (RANS) equations with a $k-{\epsilon}$ turbulent model. The one-way coupled Eulerian-Lagrangian motion of the air flow and sand particles is employed to simulate the particle trajectories and particle-wall interactions on the pipe surfaces. The predicted CFD erosion magnitudes are compared with experimental data from pipe elbows. The erosion rate results do not reveal a good accordance between the simulation and experimental results. It seems that the CFD shows a slightly over-predicted erosion ratio.

• Journal of the Korean Society of Safety
• /
• v.10 no.3
• /
• pp.68-80
• /
• 1995

A series of parametric calculations have been performed in order to investigate the short-term and short-range plume and puff behavior of toxic gaseous and solid pollutant dispersion in an open atmosphere. The simulation is made by the use of the computer program developed by this laboratory, in which a control-volume based finite-difference method is used together with the SIMPLEC algorithm for the resolution of the pressure-velocity coupling appeared In Wavier-Stokes equation. The Reynolds stresses are solved by the standard two-equation k-$\varepsilon$ model modified for buoyancy together with the RNG(Renormalization Group) k-$\varepsilon$ model. The major parameters considered in this calculation are pollutant gas density and temperature, the relative velocity of pollutants to that of the surrounding atmospheric air, and particulate size and density together with the height released. The flow field is typically characterized by the formation of a strong recirculation region for the case of the low density gases such as $CH_4$ and air due to the strong buoyancy, while the flow is simply declining pattern toward the downstream ground for the case of heavy molecule like the $CH_2C1_2$and $CCl_4$, even for the high temperature, $200^{\circ}C$. The effect of gas temperature and velocity on the flow field together with the particle trajectory are presented and discussed in detail. In general, the results are physically acceptable and consistent.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.4
• /
• pp.485-495
• /
• 2001

This paper presents a similarity solution for the directional casting of binary peritectic alloys in the presence of shrinkage-induced flow. The present model retains essential ingredients of alloy solidification, such as temperature-solute coupling, macrosegregation, solid-liquid property differences, and finite back diffusion in the primary phase. An algorithm for simultaneously determining the peritectic and liquidus positions is newly developed, which proves to be more efficient and stable than the existing scheme. Sample calculations are performed for both hypo- and hyper-peritectic compositions. The results show that the present analysis is capable of properly resolving the solidification characteristics of peritectic alloys so that it can be used for validating numerical models as a test solution.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.47 no.10
• /
• pp.720-727
• /
• 2019

Resonating thermal lags of solid fuel with heat transfer oscillations generated by boundary layer oscillation is the primary mechanism of the occurrence of the LFI (Low Frequency Combustion Instability) in hybrid rocket combustion. This study was experimentally attempted to confirm that how the boundary layer was perturbed and led to the LFI. Special attention was also made on oxidizer swirl injection to investigate the contribution to combustion stabilization. Also the overall behavior of fluctuating boundary layer flow and the occurrence of the LFI was monitored as swirl intensity increased. Fluctuating boundary layer was successfully monitored by the captured image and POD (Proper Orthogonal Decomposition) analysis. In the results, oscillating boundary layer became stabilized as the swirl intensity increases. And the coupling strength between high frequency p', q' diminished and periodical amplification of RI (Rayleigh Index) with similar frequency band of thermal lag was also decreased. Thus, results confirmed that oscillating axial boundary layer triggered by periodic coupling of high frequency p', q' is the primary mechanism to excite thermal resonance with thermal lag characteristics of solid fuel.

• Journal of the Korean Society of Safety
• /
• v.9 no.1
• /
• pp.146-154
• /
• 1994

A series of numerical calculations are performed in order to investigate the dispersion mechanism of toxic gaseous and solid pollutants in extremely short-term and short range. The calculations are carried out in an open space characterized by turbulent boundary layer. The simulation is made by the use of numerical model, in which a control-volume based finite difference method is used together with the SIMPLEC algorithm for the resolution of the pressure-velocity coupling problem. The Reynolds stresses are solved by two-equation, k-$\varepsilon$ model modified for buoyancy. The major parameters consider-ed in this study are temperature, velocity and Injection height of toxic gases, environmental conditions such as temperature and velocity of free stream air, and topographic factor. The results are presented and discussed in detail. The flow field is commonly characterized by the formation of a strong recirculation zone due to the upward motion of the hot toxic gas and ground shear stress. The driving force of the upward motion is explained by the effect of thermal buoyancy of hot gas and the difference of inlet velocity between toxic gas and free stream.

• Structural Engineering and Mechanics
• /
• v.46 no.2
• /
• pp.197-212
• /
• 2013

With more and more high-rise building being constructed in recent decades, bluff body flow with high Reynolds number and large scale dimensions has become an important topic in theoretical researches and engineering applications. In view of mechanics, the key problems in such flow are high Reynolds number turbulence and fluid-solid interaction. Aiming at such problems, a parallel fluid-structure interaction method based on socket parallel architecture was established and combined with the methods and models of large eddy simulation developed by authors recently. The new method is validated by the full two-way FSI simulations of 1:375 CAARC building model with Re = 70000 and a full scale Taipei101 high-rise building with Re = 1e8, The results obtained show that the proposed method and models is potential to perform high-Reynolds number LES and high-efficiency two-way coupling between detailed fluid dynamics computing and solid structure dynamics computing so that the detailed wind induced responses for high-rise buildings can be resolved practically.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.3
• /
• pp.426-434
• /
• 2000

This paper deals with a generalized similarity solution for the one-dimensional solidification of ternary or higher-order multicomponent alloys. The present approach not only retains the existing features of binary systems such as temperature- solute coupling, shrinkage-induced flow, solid-liquid property differences, and finite back diffusion, but also is capable of handling a multicomponent alloy without restrictions on the partition coefficient and microsegregation parameter. For an alloy of N-solute species, governing equations in the mushy region reduce to (N+2) nonlinear ordinary differential equations via similarity transformation, which are to be solved along with the closed-form solutions for the solid and liquid regions. A linearized correction scheme adopted in the solution procedure facilitates to determine the solidus and liquidus positions stably. The result for a sample ternary alloy agrees excellently with the numerical prediction as well as the reported similarity solution. Additional calculations are also presented to show the utility of this study. Finally, it is concluded that the present analysis includes the previous analytical approaches as subsets.

• Nuclear Engineering and Technology
• /
• v.54 no.3
• /
• pp.849-859
• /
• 2022

Feed-water valves (FWVs) are used to regulate the flow rate of water entering steam generators, which are very important devices in nuclear power plants. Due to the working environment of relatively high pressure and temperature, there is strength failure problem of valve body in some cases. Based on the thermo-fluid-solid coupling model, the valve body stress of the feed-water valve in the opening process is investigated. The flow field characteristics inside the valve and temperature change of the valve body with time are studied. The stress analysis of the valve body is carried out considering mechanical stress and thermal stress comprehensively. The results show that the area with relatively high-velocity area moves gradually from the bottom of the cross section to the top of the cross section with the increase of the opening degree. The whole valve body reaches the same temperature of 250 ℃ at the time of 1894 s. The maximum stress of the valve body meets the design requirements by stress assessment. This work can be referred for the design of FWVs and other similar valves.