• 한국소음진동공학회논문집
• /
• 제16권9호
• /
• pp.937-948
• /
• 2006

In this study, a fluid/structure coupled analysis system for simulating complex flow-induced vibration (FIV) phenomenon of cascades has been developed. The flow is modeled using Euler and Wavier-Stokes equations with different turbulent models. The fluid domains are modeled using the unstructured grid system with dynamic deformations due to the motion of structural boundary. The Spalart-Allmaras (S-A) and the SST ${\kappa}-{\omega}$ turbulent models are used to predict the transonic turbulent flows. A fully implicit time marching scheme based on the Newmark direct integration method is used in order to solve the coupled governing equations for viscous flow-induced vibration phenomena. For the purpose of validation for the developed FIV analysis system, comparison results for computational analyses of steady and unsteady aerodynamics and flutter analyses are presented in the transonic flow region. In addition, flow-induced vibration analyses for the isolated cascade and multi-blades cascade models have been conducted to show the physical fluid-structure interaction effects in the time domain.

• 한국소음진동공학회:학술대회논문집
• /
• 한국소음진동공학회 2006년도 춘계학술대회논문집
• /
• pp.793-802
• /
• 2006

In this study, a fluid/structure coupled analysis system for simulating complex flow-induced vibration (FIV) phenomenon of cascades has been developed. The flow is modeled using Euler and Wavier-Stokes equations with different turbulent models. The fluid domains are modeled using the unstructured grid system with dynamic deformations due to the motion of structural boundary. The Spalart-Allmaras (S-A) and the SST ${\kappa}-{\omega}$ turbulent models are used to predict the transonic turbulent flows. A fully implicit time marching scheme based on the Newmark direct integration method is used in order to solve the coupled governing equations for viscous flow-induced vibration phenomena. For the purpose of validation for the developed FIV analysis system, comparison results for computational analyses of steady and unsteady aerodynamics and flutter analyses are presented in the transonic flow region. In addition, flow-induced vibration analyses for the isolated cascade and multi-blades cascade models have been conducted to show the physical fluid-structure interaction effects in the time domain.

• Steel and Composite Structures
• /
• 제49권5호
• /
• pp.587-602
• /
• 2023

This article presents the numerical modelling of transient heat transfer in highly heterogeneous composite materials where the thermal conductivity, specific heat and density are assumed to be directional-dependent. This article uses a coupled finite element-finite difference scheme to perform the transient heat transfer analysis of unidirectional (1D) and multidirectional (2D/3D) functionally graded composite panels. Here, 1D/2D/3D functionally graded structures are subjected to nonuniform heat source and inhomogeneous boundary conditions. Here, the multidirectional functionally graded materials are modelled by varying material properties in individual or in-combination of spatial directions. Here, fully spatial-dependent material properties are evaluated using Voigt's micromechanics scheme via multivariable power-law functions. The weak form is obtained through the Galerkin method and solved further via the element-space and time-step discretisation through the 2D-isoparametric finite element and the implicit backward finite difference schemes, respectively. The present model is verified by comparing it with the previously reported results and the commercially available finite element tool. The numerous illustrations confirm the significance of boundary conditions and material heterogeneity on the transient temperature responses of 1D/2D/3D functionally graded panels.

• 한국수자원학회논문집
• /
• 제44권6호
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.

• 한국연소학회:학술대회논문집
• /
• 한국연소학회 2002년도 제24회 KOSCO SYMPOSIUM 논문집
• /
• pp.123-132
• /
• 2002

A numerical study was conducted to investigate the combustion phenomena of normal start and unstart processes based on ISL's RAMAC 30 experiments with different diluent amounts and fill pressures in a ram accelerator. The initial projectile launching speed was 1.8 km/s which corresponded to the superdetonative speed of the stoichiometric $H_2/O_2$ mixture diluted with 5 $CO_2$ or 4 $CO_2$. Experiments with same condition except for projectile surface material demonstrated that ignition was successful with an aluminum projectile, but no combustion was observed in case of a steel projectile. In this study, it was found that neither shock nor viscous heating was sufficient to ignite the mixture at a low speed of 1.8 km/s, as was found in the experiments using a steel projectile. However, we could succeed in igniting the mixtures by imposing a minimal amount of additional heat to the combustor section and simulate the normal start and unstart processes found in the experiments with an aluminum projectile. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations coupled with a Baldwin-Lomax turbulence model and detailed chemistry reaction equations of $H_2/O_2/CO_2$ suitable for high-pressure gaseous combustion were considered. The governing equations were discretized by a high order accurate upwind scheme and solved in a fully coupled manner with a fully implicit, time accurate integration method. The numerical results matched almost exactly to the experimental results. As a result, it was found that the normal start and unstart processes depended on the strength of gas mixture, development of shock-induced combustion wave stabilized by the first separation bubble, and its size and location.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 학술대회
• /
• pp.1-13
• /
• 2008

An efficient segregated algorithm for the coupling of velocity and pressure of incompressible fluid flow, called IDEAL Inner Doubly-Iterative Efficient Algorithm for Linked-Equations), has been proposed by the present authors. In the algorithm there exist double inner iterative processes for pressure equation at each iteration level, which almost completely overcome two approximations in SIMPLE algorithm. Thus the coupling between velocity and pressure is fully guaranteed, greatly enhancing the convergence rate and stability of solution process. The performance of the IDEAL algorithm for three-dimensional incompressible fluid flow and heat transfer problems is analyzed and a systemic comparison is made between the algorithm and three other most widely-used algorithms (SIMPLER, SIMPLEC and PISO). It is found that the IDEAL algorithm is the most robust and the most efficient one among the four algorithms compared. This new algorithm is used for the velocity prediction of a new interface capturing method. VOSET, also proposed by the present author. It is found that the combination of VOSET and IDEAL can appreciably enhance both the interface capture accuracy and convergence rate of computations.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년 추계학술대회논문집
• /
• pp.1-13
• /
• 2008

An efficient segregated algorithm for the coupling of velocity and pressure of incompressible fluid flow, called IDEAL (Inner Doubly-Iterative Efficient Algorithm for Linked-Equations), has been proposed by the present authors. In the algorithm there exist double inner iterative processes for pressure equation at each iteration level, which almost completely overcome two approximations in SIMPLE algorithm. Thus the coupling between velocity and pressure is fully guaranteed, greatly enhancing the convergence rate and stability of solution process. The performance of the IDEAL algorithm for three-dimensional incompressible fluid flow and heat transfer problems is analyzed and a systemic comparison is made between the algorithm and three other most widely-used algorithms (SIMPLER, SIMPLEC and PISO). It is found that the IDEAL algorithm is the most robust and the most efficient one among the four algorithms compared. This new algorithm is used for the velocity prediction of a new interface capturing method -VOSET, also proposed by the present author. It is found that the combination of VOSET and IDEAL can appreciably enhance both the interface capture accuracy and convergence rate of computations.

• 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.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2007년도 춘계학술대회B
• /
• pp.3003-3008
• /
• 2007

Critical nozzle has been frequently employed to measure the flow rate of various gases, but hydrogen gas, especially being at high-pressure condition, was not nearly dealt with the critical nozzle due to treatment danger. According to a few experimental data obtained recently, it was reported that the discharge coefficient of hydrogen gas through the critical nozzle exceeds unity in a specific range of Reynolds number. No detailed explanation on such an unreasonable value was made, but it was vaguely inferred as real gas effects. For the purpose of practical use of high-pressure hydrogen gas, systematic research is required to clarify the critical nozzle flow of high-pressure hydrogen gas. In the present study, a computational fluid dynamics(CFD) method has been applied to predict the critical nozzle flow of high-pressure hydrogen gas. Redlich-Kwong equation of state that take account for the forces and volume of molecules of hydrogen gas were incorporated into the axisymmetric, compressible Navier-Stokes equations. A fully implicit finite volume scheme was used to numerically solve the governing equations. The computational results were validated with some experimental data available. The results show that the coefficient of discharge coefficient is mainly influenced by the compressibility factor and the specific heat ratio, which appear more remarkable as the inlet total pressure of hydrogen gas increases.

• 한국소음진동공학회논문집
• /
• 제19권8호
• /
• pp.764-772
• /
• 2009

In this study, fluid/structure coupled analyses have been conducted for 3-D stator and rotor configuration. Advanced computational analysis system based on computational fluid dynamics(CFD) and computational structural dynamics(CSD) has been developed in order to investigate fluid/structure responses of general stator-rotor configurations. To solve the fluid/structure coupled problems, fluid domains are modeled using the structural grid system with dynamic moving and local deforming techniques. Reynolds-averaged Navier-Stokes equations with Spalart-Allmaras(S-A) and SST ${\kappa}-{\omega}$ turbulence models are solved for unsteady flow problems. A fully implicit time marching scheme based on the Newmark direct integration method is used for computing the coupled aeroelastic governing equations of the 3-D turbine blades for fluid-structure interaction(FSI) problems. Detailed fluid/structure analysis responses for stator-rotor interaction flow conditions are presented to show the physical performance and flow characteristics.