• 유체기계공업학회:학술대회논문집
• /
• 유체기계공업학회 1999년도 유체기계 연구개발 발표회 논문집
• /
• pp.179-186
• /
• 1999

Three-dimensional flow analysis is implemented to investigate the flow through transonic axial-flow compressor rotor(NASA R67), and to evaluate the performances of k-$\epsilon$ and Baldwin-Lomax turbulence models. A finite volume method is used for spatial discretization. And, the equations are solved implicitly in time with the use of approximate factorization. Upwind difference scheme is used for inviscid terms, but viscous terms are centrally differenced. The flux-difference-splitting of Roe is used to obtain fluxes at the cell faces. Numerical analysis is performed near peak efficiency and near stall. And, the results are compared with the experimental data for NASA R67 rotor. Blade-to-Blade Mach number distributions are compared to confirm the accuracy of the code. From the results, we conclude that k-$\epsilon$ model is better for the calculation of flow rate and efficiency than Baldwin-Lomax model. But, the predictions for Mach number and shock structure are almost same.

• 한국해양공학회지
• /
• 제19권5호
• /
• pp.78-82
• /
• 2005

The motion response results of a submerged submarine section in waves are presented. The numerical method is based on Cauchy's integral and 3 degrees-of-freedom motions of submarine sections are calculated in two dimensions, in regular waves. The fully nonlinear free surface and body boundary conditions are applied to the present problem, and the viscous effects on the submarine are modeled by Morison's formulas. The motions of submarine sections in beam sea are directly simulated and the effects of wave frequency, snorkel depth, and bridge are discussed.

• Journal of Mechanical Science and Technology
• /
• 제20권8호
• /
• pp.1284-1291
• /
• 2006

Single phase pressure drops in micro tubes were investigated through an experimental measurement and a numerical simulation. Experimental Po was obtained in circular micro tubes with 87 and $118{\mu}m$ diameter with distilled water. Experiments were carried out in laminar flow region with varying the Re 15-450 for the $87{\mu}m$ diameter tubes and 60-1300 for the $118{\mu}m$ diameter tube. No early transition from laminar to turbulent flow was detected for the experimental range. The computational estimation of pressure drop in the $87{\mu}m$ diameter tube was performed with the aid of CFD software. Boundary conditions from experiments were used for the numerical simulation. The results of experimental and numerical studies showed a good agreement with the conventional macro theory.

• Tribology and Lubricants
• /
• 제13권2호
• /
• pp.89-95
• /
• 1997

This paper considers two dimensional analysis for lubrication between the single piston ring and the cylinder liner. The piston ring is treated as a reciprocating, hydrodynamic bearing with combined sliding and squeeze motion. Reynolds' equation is used, to model lubrication with Reynolds' cavitation boundary condition. This analysis is developed to get the cyclic variation of minimum film thickness and viscous frictional force. Two types of piston ring face shape are considered. This result can be used to study the influence of ring shape design parameter to improve the characteristics of sealing and lubrication.

• Advances in aircraft and spacecraft science
• /
• 제9권5호
• /
• pp.415-431
• /
• 2022

Secondary flows have a huge impact on losses generation in modern low pressure gas turbines (LPTs). At design point, the interaction of the blade profile with the end-wall boundary layer is responsible for up to 40% of total losses. Therefore, predicting accurately the end-wall flow field in a LPT is extremely important in the industrial design phase. Since the inlet boundary layer profile is one of the factors which most affects the evolution of secondary flows, the first main objective of the present work is to investigate the impact of two different inlet conditions on the end-wall flow field of the T106A, a well known LPT cascade. The first condition, labeled in the paper as C1, is represented by uniform conditions at the inlet plane and the second, C2, by a flow characterized by a defined inlet boundary layer profile. The code used for the simulations is based on the Discontinuous Galerkin (DG) formulation and solves the Reynolds-averaged Navier-Stokes (RANS) equations coupled with the Spalart Allmaras turbulence model. Secondly, this work aims at estimating the influence of viscosity and turbulence on the T106A end-wall flow field. In order to do so, RANS results are compared with those obtained from an inviscid simulation with a prescribed inlet total pressure profile, which mimics a boundary layer. A comparison between C1 and C2 results highlights an influence of secondary flows on the flow field up to a significant distance from the end-wall. In particular, the C2 end-wall flow field appears to be characterized by greater over turning and under turning angles and higher total pressure losses. Furthermore, the C2 simulated flow field shows good agreement with experimental and numerical data available in literature. The C2 and inviscid Euler computed flow fields, although globally comparable, present evident differences. The cascade passage simulated with inviscid flow is mainly dominated by a single large and homogeneous vortex structure, less stretched in the spanwise direction and closer to the end-wall than vortical structures computed by compressible flow simulation. It is reasonable, then, asserting that for the chosen test case a great part of the secondary flows details is strongly dependent on viscous phenomena and turbulence.

• 한국전산유체공학회지
• /
• 제21권4호
• /
• pp.11-18
• /
• 2016

A numerical study is conducted on the optical fiber coating flow in a primary coating nozzle consisting of three major parts: a resin chamber, a contraction and a coating die of small diameter. The flow is driven by the optical fiber penetrating the center of the nozzle at a high speed. The axisymmetric two-dimensional flow and heat transfer induced by viscous heating are examined based on the laminar flow assumption. Numerical experiments are performed with varying the convergent angle of nozzle contraction and the optical fiber drawing speed. The numerical results show that for high drawing speed greater than 30 m/s, there is a transition in the essential flow features depending on the convergent angle. For a large convergent angle greater than $30^{\circ}$, unfavorable multicellular flow structures are monitored, which could be associated with wall boundary-layer separation. In the regime of small convergent angle, as the angle increases, the highest resin temperature at the exit of die and the coating thickness decrease but the sensitivity of coating thickness on drawing speed and the maximum shear strain of resin on the optical fiber increase. The effects of the convergent angle are discussed in view of compromise searching for an appropriate angle for high-speed optical fiber coating.

• 한국분말야금학회:학술대회논문집
• /
• 한국분말야금학회 1997년도 추계학술강연 및 발표대회 강연 및 발표논문 초록집
• /
• pp.7-7
• /
• 1997

During sintering of very porous green bodies, as obtained by compaction of hard powders - such as tungsten carbide or ceramics - or by injection moulding, important shrinkage occurs. Due to heterogeneous green density field, gravity effects, friction on the support, thermal gradients, etc., this shrinkage is often non-uniform, which' may induce significant shape changes. As the ratio of compact dimension to powder size is very high, the mechanics of continuum is relevant to model such phenomena. Thus numerical techniques, such as the finite element method can be used to simulate the sintering process and predict the final shape of the sintered part. Such type of simulation has much been developed in the last decade firstly for hot isostatic pressing and next for die compaction. Finite element modelling has been recently applied to free sintering. The simulation of sintering should be based on constitutive equations describing the thermo-mechanical behaviour of the material under any state of stress and any temperature which may arise within the sintering body. These equations can be drawn either from experimental data or from micromechanical models. The experiments usually consist in free sintering and sinter-forging tests. Indeed applying more complex loading conditions at high temperature under controlled atmosphere is delicate. Micromechanical models describe the constitutive behaviour of aggregates of spheres from the deformation of two-sphere contact either by viscous flow or grain boundary diffusion. Such models are not able to describe complex microstructure and mechanisms as observed in real materials but they can give some basic information on the formulation of constitutive equations. Practically both experimental and theoretical approaches can be coupled to identify the constitutive equations. Such procedure has been performed for modelling the sintering of compacts obtained by die pressing of a mixture of tungsten carbide and cobalt powders. The constitutive behaviour of this material during sintering has been described by a linear viscous constitutive model, whose functions have been fitted from results of free sintering and sinter-forging experiments. This model has next been introduced in ABAQUS finite element code to simulate the sintering of heterogeneous green compacts of various geometries at constant temperature. Examples of simulations are shown and compared with experiments.

• 한국음향학회지
• /
• 제33권6호
• /
• pp.358-365
• /
• 2014

선박용 프로펠러가 수중에서 빠른 속도로 회전할 때, 날개 표면의 국부적인 압력이 낮아짐에 따라 불가피하게 여러 형태의 공동이 발생한다. 이러한 공동현상은 근본적으로 날개 단면의 기하학적 형상 특성과 수동력학적 운용조건에 의해 결정되며, 결과적으로 선박 프로펠러에서 유기되는 수중방사소음은 공동의 발생특성과 직결된다고 할 수 있다. 따라서 저소음 프로펠러 설계를 위해서는 날개 단면의 형상에 따라 발생하는 공동과 그에 따른 소음특성을 이론 및 실험을 통해 정량적으로 평가할 수 있어야 한다. 본 연구에서는 저소음 프로펠러의 설계단계에서부터 적용이 가능한 부분공동성능 해석법 개발 및 날개단면 형상정보 도출을 목표로 선박용 프로펠러 날개 단면에서 발생하는 부분공동 다상 유동의 비점성 수치해석을 수행하였다. 또한 점성해석 상용프로그램인 FLUENT에서 제공하는 난류 및 공동 모델 조합에 따른 결과를 살펴보았으며, 점성 및 비점성 해석 결과를 비교, 평가 하였다.

• 한국해양공학회지
• /
• 제35권4호
• /
• pp.296-304
• /
• 2021

Research on the development of marine renewable energy is actively in progress. Various studies are being conducted on the development of wave energy converters. In this study, a numerical analysis of wave-energy extraction performance was performed according to the body shape and scale of the sloped oscillating water column (OWC) wave energy converter (WEC), which can be connected with the breakwater. The sloped OWC WEC was modeled in the time domain using a two-dimensional fully nonlinear numerical wave tank. The nonlinear free surface condition in the chamber was derived to represent the pneumatic pressure owing to the wave column motion and viscous energy loss at the chamber entrance. The free surface elevations in the sloped chamber were calculated at various incident wave periods. For verification, the results were compared with the 1:20 scaled model test. The maximum wave energy extraction was estimated with a pneumatic damping coefficient. To calculate the energy extraction of the actual size WEC, OWC models approximately 20 times larger than the scale model were calculated, and the viscous damping coefficient according to each size was predicted and applied. It was verified that the energy, owing to the airflow in the chamber, increased as the incident wave period increased, and the maximum efficiency of energy extraction was approximately 40% of the incident wave energy. Under the given incident wave conditions, the maximum extractable wave power at a chamber length of 5 m and a skirt draft of 2 m was approximately 4.59 kW/m.

• 한국지반공학회:학술대회논문집
• /
• 한국지반공학회 2003년도 봄 학술발표회 논문집
• /
• pp.257-264
• /
• 2003

Tunnelling in water bearing soils influences the ground water regime. It has been indicated in the literature that the existence of ground water above a tunnel influences tunnel stability and the settlement profile. Only limited research, however, has been done on ground water movements around tunnels and their influence on tunnel performance. Time dependent soil behaviour can be caused by the changes of pore water pressure and/or the viscous properties of soil(creep) under the stress change resulting from the advance of the tunnel face. De Moor(1989) demonstrated that the time dependent deformations due to tunnelling are mainly the results of pore pressure dissipation and should be interpreted in terms of effective stress changes. Drainage into tunnels is governed by the permeability of the soil, the length of the drainage path and the hydraulic boundary conditions. The potential effect of lime dependent settlement in a shallow tunnel is likely to occur rapidly due to the short drainage path and possibly high coefficient of consolidation. Existing 2D modelling methods are not applicable to these tunnelling problems, as it is difficult to define empirical parameters. In this paper the time-based 2D modelling method is adopted to account for the three dimensional effect and time dependent behaviour during tunnel construction. The effect of coupling between the unloading procedure and consolidation during excavation is profoundly investigated with the method. It is pointed out that realistic modelling can be achieved by defining a proper permeability at the excavation boundary and prescribing appropriate time for excavation Some guidelines for the numerical modelling of drained and undrained excavation has been suggested using characteristic time factor. It is highlighted that certain range of the factor shows combined effect between the unloading procedure due to excavation and consolidation during construction.