• Structural Engineering and Mechanics
• /
• 제66권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
• /
• 제20권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.

• 한국해양공학회지
• /
• 제31권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.

• 한국안전학회지
• /
• 제10권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.

• 대한기계학회논문집B
• /
• 제25권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.

• 한국항공우주학회지
• /
• 제47권10호
• /
• pp.720-727
• /
• 2019

하이브리드 로켓 연소의 저주파수 연소불안정은 고체연료의 열적지연(Thermal Lag)과 경계층 유동 변화에 의한 열전달 진동의 공진에 의해 발생한다. 본 연구는 연료 표면 근처의 경계층 유동의 교란이 어떤 물리적 과정에 의해 발생하여 연소불안정으로 발달하는지를 실험적으로 확인하였다. 특히 산화제의 스월 분사는 연소 안정화에 매우 큰 기여를 하므로 스월 강도를 증가시키며 경계층의 변화와 연소불안정의 발생과정을 연구하였다. 경계층 섭동을 확인하기 위하여 연소 유동장을 가시화하였고 이미지에 대한 POD(Proper Orthogonal Decomposition) 분석을 시도하였다. 스월 강도가 증가할수록 500Hz 대역 고주파수 p', q'의 결합이 약해지며 열적지연과 유사한 주파수 특성을 갖는 Rayleigh Index의 섭동 발생도 약해져 경계층 진동의 발생이 점차 감소하는 것을 관찰하였다. 따라서 고주파수 p', q'의 주기적인 결합에 의한 축 방향 경계층 진동이 나타나면 열적지연 주파수와 공진에 의한 연소불안정이 발생함을 확인하였다.

• 한국안전학회지
• /
• 제9권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
• /
• 제46권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.

• 대한기계학회논문집B
• /
• 제24권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
• /
• 제54권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.