• International Journal of Fluid Machinery and Systems
• /
• 제3권4호
• /
• pp.342-351
• /
• 2010

During the operation of a hydro turbine the fluid mechanical pressure loading on the turbine blades provides the driving torque on the turbine shaft. This fluid loading results in a structural load on the component which in turn causes the turbine blade to deflect. Classically, these mechanical stresses and deflections are calculated by means of finite element analysis (FEA) which applies the pressure distribution on the blade surface calculated by computational fluid dynamics (CFD) as a major boundary condition. Such an approach can be seen as a one-way coupled simulation of the fluid structure interaction (FSI) problem. In this analysis the reverse influence of the deformation on the fluid is generally neglected. Especially in axial machines the blade deformation can result in a significant impact on the turbine performance. The present paper analyzes this influence by means of fully two-way coupled FSI simulations of a propeller turbine utilizing two different approaches. The configuration has been simulated by coupling the two commercial solvers ANSYS CFX for the fluid mechanical simulation with ANSYS Classic for the structure mechanical simulation. A detailed comparison of the results for various blade stiffness by means of changing Young's Modulus are presented. The influence of the blade deformation on the runner discharge and performance will be discussed and shows for the configuration investigated no significant influence under normal structural conditions. This study also highlights that a two-way coupled fluid structure interaction simulation of a real engineering configuration is still a challenging task for today's commercially available simulation tools.

• Geomechanics and Engineering
• /
• 제28권5호
• /
• pp.437-454
• /
• 2022

The hydro-mechanical (HM) two-way sequential coupling in the TOUGH2-FLAC3D linking algorithm is validated completely and successfully in both M to H and H to M directions, which are initiated by mechanical surface loading for geomechanical validation and hydrological groundwater pumping for hydrogeological validation, respectively. For such complete and successful validation, a TOUGH2-FLAC3D linked numerical model is developed first by adopting the TOUGH2-FLAC3D linking algorithm, which uses the two-way (fixed-stress split) sequential coupling scheme and the implicit backward time stepping method. Two geomechanical and two hydrogeological validation problems are then simulated using the linked numerical model together with basic validation strategies and prerequisites. The second geomechanical and second hydrogeological validation problems are also associated with the Mandel effect and the Noordbergum and Rhade effects, respectively, which are three phenomenally well-known but numerically challenging HM effects. Finally, sequentially coupled numerical solutions are compared with either analytical solutions (verification) or fully coupled numerical solutions (benchmarking). In all the four validation problems, they show almost perfect to extremely or very good agreement. In addition, the second geomechanical validation problem clearly displays the Mandel effect and suggests a proper or minimum geometrical ratio of the height to the width for the rectangular domain to maximize agreement between the numerical and analytical solutions. In the meantime, the second hydrogeological validation problem clearly displays the Noordbergum and Rhade effects and implies that the HM two-way sequential coupling scheme used in the linked numerical model is as rigorous as the HM two-way full coupling scheme used in a fully coupled numerical model.

• Geomechanics and Engineering
• /
• 제10권4호
• /
• pp.423-436
• /
• 2016

Currently, there is a great interest in the coupling between multiphase fluid flow and geomechanical effects in hydrocarbon reservoirs and surrounding rocks. The ideal solution for this coupled problem is to introduce the geomechanical effects through the stress analysis solution and implement an algorithm, which assures that the equations governing the flow and stress analyses are obeyed in each time step. This paper deals with the implementation of a program (FORTRAN90 interface code), which was developed to couple conventional reservoir (ECLIPSE) and geomechanical (ABAQUS) simulators, using a partial coupling algorithm. The explicit coupled hydro-mechanical behavior of Iranian field during depletion and $CO_2$ injection is studied using the soils consolidation procedure available in ABAQUS. Time dependent reservoir pressure 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. The FEM analysis of the reservoir showed no sign of plastic strain under production and $CO_2$ injection scenarios in any part of the reservoir and the stress paths do not show a critical behavior.

• 한국지반환경공학회 논문집
• /
• 제11권5호
• /
• pp.35-43
• /
• 2010

불포화 투수계수는 모관흡수력(또는 함수비)과 간극률의 함수로 정의되어야 한다. 그러나 기존 상용프로그램이나 문헌에 개발된 모델들은 모관흡수력 만의 함수로써 정의 되어왔다. 사면의 안정성은 수리학적인 측면과 흙의 전단강도 특성들이 모두 고려되어야 한다. 이 두 가지 특성에 대한 해석은 상용 프로그램인 Seep/W, Sigma/W 그리고 Slope/W(Geo-Slope, 2007) 같은 상호 연계가능한 해석 프로그램을 이용한다. 이런 해석 방법으로 강우시 불포화 침투수를 고려하고 흙의 변형을 예측하며 결국에는 사면의 안정성 평가에 예측하고 있다. 불포화 침투수를 해석하는 Seep/W 프로그램은 흙의 변형은 고려하지 못하고 단지 불포화 침투해석만을 수행한다. 그래서 흙의 변형을 고려하기 위해 Sigma/W 프로그램을 연계하여 실제 사면에서 발생하는 침투에 의한 흙의 변형을 모사하고 있다. 이와 같이 동시에 발생하는 침투와 흙의 변형을 구현하기 위해 여러 연구자들이 새로운 모델들을 개발하고 있지만, 현재로서는 각각의 해석을 연계하여 실제로 일어나는 현상에 접근하고 있는 실정이다. 본 연구는 기존 프로그램을 이용하여 연계해석에 의한 결과와 동시해석과 유사한 알고리즘을 이용하여 해석한 결과를 비교하여 동시해석의 타당성을 검증하고자 하였다.

• Geomechanics and Engineering
• /
• 제29권1호
• /
• pp.1-11
• /
• 2022

A semi-analytical solution to responses of overconsolidated (OC) unsaturated soils surrounding an expanding spherical cavity under constant suction condition is presented. To capture the elastoplastic hydro-mechanical property of OC unsaturated soils, the unified hardening (UH) model for OC unsaturated soil is adopted in corporation with a soil-water characteristic curve (SWCC) and two suction yield surfaces. Taking the specific volume, radial stress, tangential stress and degree of saturation as the four basic unknowns, the problem investigated is formulated by solving a set of first-order ordinary differential equations with the help of an auxiliary variable and an iterative algorithm. The present solution is validated by comparing with available solution based on the modified Cam Clay (MCC) model. Parametric studies reveal that the hydraulic and mechanical responses of spherical cavity expanding in unsaturated soils are not only coupled, but also affected by suction and overconsolidation ratio (OCR) significantly. More importantly, whether hydraulic yield will occur or not depends only on the initial relationship between suction yield stress and suction. The presented solution can be used for calibration of some insitu tests in OC unsaturated soil.

• 터널과지하공간
• /
• 제34권2호
• /
• pp.127-142
• /
• 2024

수치해석 기법의 강건성은 다양한 모델링 조건에서 계산 성능이 유지되는 것을 의미하며, 새로운 해석기법 또는 수치코드는 벤치마크 테스트를 통해 강건성이 평가될 필요가 있다. TOUGH-FLAC 모델링 기법은 국내외적으로 이산화탄소 지중저장, 사용후핵연료 지층처분, 지열 개발 등 다양한 분야에 적용되었으며, 실험 계측자료, 다른 수치코드들과의 결과 비교를 통해 모델링 유효성이 분석되었다. 본 연구에서는 해석해를 갖는 열-수리-역학적 복합거동 문제를 토대로 TOUGH-FLAC 기법의 벤치마크 테스트를 수행하였다. 적용된 해석해는 완전히 포화된 지반에 점열원 작용 시 주변매질의 온도, 간극수압, 역학적 거동과 관계되며, 해석해와 수치모사 결과를 비교하여 TOUGH-FLAC 기법의 강건성이 평가되었다. 또한, 열-수리-역학 해석의 연계항, 유체 상변화, 시간증분이 복합거동 계산에 미치는 영향을 조사하였다.

• 터널과지하공간
• /
• 제28권5호
• /
• pp.400-425
• /
• 2018

본 논문에서는 국제공동연구인 DECOVALEX-2019 프로젝트 Task B의 연구결과와 현황을 소개하였다. Task B의 주제는 'Fault slip modelling'으로 유체의 주입으로 인해 발생하는 단층의 재활성(미끄러짐, 전단파괴)과 수리역학적 거동을 예측할 수 있는 해석기법을 개발하는 데에 그 목적이 있다. 1단계 연구는 참가팀들이 연구주제에 대해 숙지하고, 벤치마크 모델을 대상으로 단층의 투수특성과 역학적 거동의 상호작용을 모사할 수 있는 해석코드를 개발할 수 있도록 하는 준비 단계의 연구이다. 본 연구에서는 TOUGH-FLAC 연동해석 기법을 사용하여 물 주입으로 인한 단층의 수리역학적 연계거동을 모사하였다. TOUGH2 해석에서는 단층을 Darcy의 법칙과 삼승법칙을 따르는 연속체 요소로 모델링하였으며, FLAC3D 해석에서는 미끄러짐과 개폐가 허용되는 불연속 인터페이스 요소를 통해 모사하였다. 두 가지 수리간극모델에 대하여 수리역학적 커플링 관계식을 수치화하였으며, 연속체 요소(수리모델)와 인터페이스 요소(역학모델)의 거동을 연계할 수 있는 해석기법을 제시하였다. 또한, 단층의 역학적 변형(간극의 변화)으로 인한 수리물성 변화와 기하학적 변화(해석 메쉬의 변형)를 수리해석에 반영할 수 있는 해석기법을 개발하였다. 다양한 압력의 물을 단계적으로 주입하고 이로 인해 유도되는 단층의 탄성거동 및 전단파괴(미끄러짐)에 대해 살펴보았으며, 수리간극의 변화 양상과 원인, 압력 분포와 주입율의 관계 등을 면밀히 검토하였다. 해석 결과, 본 연구에서 개발한 해석기법이 물 주입으로 인한 단층의 미끄러짐 거동을 합리적인 수준에서 재현할 수 있는 것으로 판단할 수 있었다. 본 연구의 해석모델은 Task B에 참여하는 국외 연구팀들과의 의견 교류와 워크숍을 통해 지속적으로 개선하는 한편, 향후 연구의 현장시험에 적용하여 타당성을 검증할 예정이다.

• Coupled systems mechanics
• /
• 제6권4호
• /
• pp.439-464
• /
• 2017

The hydro-elastic system consisting of a pre-stretched highly elastic plate, compressible Newtonian viscous fluid, and the rigid wall is considered and it is assumed that on the plate a lineal-located time-harmonic force acts. It is required to investigate the dynamic behavior of this system and determine how the problem parameters and especially the pre-straining of the plate acts on this behavior. The elasticity relations of the plate are described through the harmonic potential and linearized (with respect to perturbations caused by external time-harmonic force) form of these relations is used in the present investigation. The plane-strain state in the plate is considered and the motion of that is described within the scope of the three-dimensional linearized equations of elastic waves in elastic bodies with initial stresses. The motion of the fluid is described by the linearized Navier-Stokes equations and it is considered the plane-parallel flow of this fluid. The Fourier transform with respect to the space coordinate is applied for a solution to the corresponding boundary-value problem. Numerical results on the frequency response of the interface normal stress and normal velocity and the influence of the initial stretching of the plate on this response are presented and discussed. In particular, it is established that the initial stretching of the plate can decrease significantly the absolute values of the aforementioned quantities.

• Coupled systems mechanics
• /
• 제8권3호
• /
• pp.199-218
• /
• 2019

This paper studies the forced vibration of the hydro-elastic system consisting of the anisotropic (orthotropic) plate, compressible viscous fluid and rigid wall within the scope of the exact equations and relations of elastodynamics for anisotropic bodies for describing of the plate motion, and with utilizing the linearized exact Navier-Stokes equations for describing of the fluid flow. For solution of the corresponding boundary value problem it is employed time-harmonic presentation of the sought values with respect to time and the Fourier transform with respect to the space coordinate on the coordinate axis directed along the plate length. Numerical results on the pressure acting on the interface plane between the plate and fluid are presented and discussed. The main aim in this discussion is focused on the study of the influence of the plate material anisotropy on the frequency response of the mentioned pressure. In particular, it is established that under fixed values of the shear modulus of the plate material a decrease in the values of the modulus of elasticity of the plate material in the direction of plate length causes to increase of the absolute values of the interface pressure. The numerical results are presented not only for the viscous fluid case but also for the inviscid fluid case.

• Coupled systems mechanics
• /
• 제10권2호
• /
• pp.123-142
• /
• 2021

The paper deals with the study of the dispersion of quasi-Lamb waves in a hydro-elastic system consisting of an elastic plate, barotropic compressible inviscid fluid, and rigid wall. The motion of the plate is described using the exact equations of elastodynamics, however, the flow of the fluid using the linearized equations and relations of the Navier-Stokes equations. The corresponding dispersion equation is obtained and this equation is solved numerically, as a result of which the corresponding dispersion curves are constructed. The main attention is focused on the effect of the presence of the fluid and the effect of the fluid layer thickness (i.e., the fluid depth) on the dispersion curves. The influence of the problem parameters on the dispersion curves related to the quasi-Scholte wave is also considered. As a result of the analyses of the numerical results, concrete conclusions are made about the influence of the fluid depth, the rigid wall restriction on the fluid motion, and the material properties of the constituents on the dispersion curves. During the analyses, the zeroth and the first four modes of the propagating waves are considered.