• 대한기계학회논문집B
• /
• 제25권7호
• /
• pp.885-894
• /
• 2001

A large eddy simulation(LES) is performed for turbulent flow in a combustion device. The combustion device is simplified as a cylinder with sudden expansion. To promote turbulent mixing and to accommodate flame stability, a flame holder is attached inside the combustion chamber. Emphasis is placed on the flow details with different geometries of the flame holder. The subgrid scale models are applied and validated. The simulation code is constructed by using a general coordinate system based on the physical contravariant velocity components. The calculated Reynolds numbers are 5000 and 50000 based on the bulk velocity and the diameter of inlet pipe. The predicted turbulent statistics are evaluated by comparing with the LDV measurement data. The agreement of LES with the experimental data is shown to be satisfactory.

• Journal of Advanced Marine Engineering and Technology
• /
• 제28권1호
• /
• pp.98-108
• /
• 2004

This study concerns a large eddy simulation (LES) of turbulent flow around a bluff body inside a sudden expansion cylindrical chamber, a configuration which resembles a premixed gas turbine combustor The simulation code is constructed by using the general coordinate system based on the physical contravariant velocity components. The Smagorinsky model is employed and the calculated Reynolds number is 5,000 based on the bulk velocity and the diameter of the inlet pipe. The combined grid technique and cylindrical grid are tested in the numerical simulation with complex geometry. The predicted turbulent statistics are evaluated by comparing with LDV measurement data. The numerical flow visualizations depict the behavior of turbulent mixing process behind the flame holder.

• Journal of Advanced Marine Engineering and Technology
• /
• 제28권3호
• /
• pp.423-433
• /
• 2004

A large eddy simulation (LES) is performed for turbulent flow in a combustion device. The combustion device is simplified as a cylindrical chamber with sudden expansion. A flame holder is attached inside a cylindrical chamber in order to promote turbulent mixing and to accommodate flame stability. The turbulent sub-grid scale models are applied and validated. Emphasis is placed on the evaluation of turbulent model for the LES of complex geometry. The simulation code is constructed by using a general coordinate system based on the physical contravariant velocity components. The calculated Reynolds number is 5000 based on the bulk velocity and the diameter of inlet pipe. The predicted turbulent statistics are evaluated by comparing with the LDV measurement data. The Smagorinsky model coefficients are estimated and the utility of dynamic SGS models are confirmed in the LES of complex geometry.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2008년도 12th Asian Congress of Fluid Mechanics
• /
• pp.60.3-60.3
• /
• 2008

• 대한기계학회논문집B
• /
• 제37권6호
• /
• pp.591-597
• /
• 2013

대부분의 비뉴턴 유체 모델의 경우 격자볼쯔만 법을 이용한 3차원 유동해석에 효율적으로 적용하기 어려운 문제점을 가지고 있다. 이 연구에서는 이를 해결하기 위해 개발된 격자볼쯔만법 전용 비뉴턴 유체 모델인 Hydro-Kinetic 모델을 3차원 격자볼쯔만법 해석에 적용할 수 있도록 수정하고 이를 3차원 급축소 및 급확대 유동에 적용하여 얻은 결과를 통해서 강한 전단유동장에서 HK모델 유체의 거동을 분석하였다. HK모델은 변형률과 완화시간과의 관계를 나타낸 모델로 강한 전단유동장에서 국소적으로 큰 변형률이 발생하는 경우 완화시간이 급격하게 감소하여 해석이 불안정해지는 경향이 있어 격자볼쯔만법 해석안정성을 확보하기 위해서는 HK모델의 인수 ${\Gamma}$와 완화시간을 일정한 구간에서 변화하도록 제한 할 필요가 있다.

• Interaction and multiscale mechanics
• /
• 제6권2호
• /
• pp.157-172
• /
• 2013

This paper presented an investigation of macromolecular suspension in a grooved channel by using the dissipative particle dynamics (DPD) with finitely extensible non-linear elastic (FENE) bead spring chains model. Before studying the movement and evolution of macromolecules, the DPD method was first validated by modeling the simple fluid flow in the grooved channel. For both simple fluid flow and macromolecular suspension, the flow fields were analyzed in detail. It is found that the structure of the grooved channel with sudden contraction and expansion strongly affects the velocity distribution. As the width of the channel reduces, the horizontal velocity increases simultaneously. Vortices can also be found at the top and bottom corners behind the contraction section. For macromolecular suspension, the macromolecular chains influence velocity and density distribution rather than the temperature and pressure. Macromolecules tend to drag simple fluid particles, reducing the velocity with density and velocity fluctuations. Particle trajectories and evolution of macromolecular conformation were investigated. The structure of the grooved channel with sudden contraction and expansion significantly influence the evolution of macromolecular conformation, while macromolecules display adaptivity to adjust their own conformation and angle to suit the structure so as to pass the channel smoothly.

• International Journal of Air-Conditioning and Refrigeration
• /
• 제7권
• /
• pp.1-10
• /
• 1999

In this study, a numerical analysis is made on the ice-formation for laminar water flow inside a stenotic tube. The study takes into account the interaction between the laminar flow and the stenotic port in the circular tube. The purpose of the present numerical investigation is to assess the effect of a stenotic shape on the instantaneous shape of the flow passage during freezing upstream/downstream of the stenotic channel. In the solution strategy, the present study is substantially distinguished from the existing works in that the complete set of governing equations in both the solid and liquid regions are resolved. In a channel flow between parallel plates, the agreement between the of predictions and the available experimental data is very good. Numerical analyses are performed for parametric variations of the position and heights of stenotic shape and flow rate. The results show that the stenotic shape has the great effect on the thickness of the solidification layer inside the tube. As the height of a stenosis grows and the length of a stenosis decreases, the ice layer thickness near the stenotic port is thinner, due to backward flow caused by the sudden expansion of a water tunnel. It is found that the flow passage has a slight uniform taper up to the stenotic channel, at which a sudden expansion is observed. It is also shown that the ice layer becomes more fat in accordance with its Reynolds number.

• 대한기계학회논문집B
• /
• 제25권7호
• /
• pp.905-914
• /
• 2001

A large eddy simulation(LES) is performed for turbulent flow around a bluff body inside a sudden expansion cylinder chamber, a configuration which resembles a premixed gas turbine combustor. To promote turbulent mixing and to accommodate flame stability, a flame holder is installed inside the combustion chamber. The Smagorinsky model is employed and the calculated Reynolds number is 5,000 based on the bulk velocity and the diameter of the inlet pipe. The simulation code is constructed by using a general coordinate system based on the physical contravariant velocity components. The predicted turbulent statistics are evaluated by comparing them with the laser-doppler velocimetry (LDV) measurement data. The agreement of LES with the experimental data is shown to be satisfactory. Emphasis is placed on the time-dependent evolutions of turbulent vortical structure behind the flame holder. The numerical flow visualizations depict the behavior of large-scale vortices. The turbulent mixing process behind the flame holder is analyzed by visualizing the sectional views of vortical structure.

• 대한설비공학회지:설비저널
• /
• 제14권4호
• /
• pp.293-299
• /
• 1985

This paper describes an economical prediction procedure for heat transfer phenomenon through a channel containing an abrupt asymmetric expansion in flow cross-seetional area. Numerical solutions for the flow field are obtained by the finite difference numerical method applied to the modified boundary layer equations. Modified boundary energy equation is used to analyze heat transfer as modified boundary momentum equation. Predictions of the method compare very favorable with exprimental data. Results of this study by modified boundary layer equation are as follows : 1. The computation time required for the scheme is at least an order of magnitude less than for the numerical solution of the full Navier-stokes and Energy eguations. 2. In laminar flow, the maximum heat transfer occurs downstream of the reattachment point.

• 한국해양공학회지
• /
• 제10권3호
• /
• pp.138-152
• /
• 1996

Many studies of heat transfer on the swirling flow or unswirled flow in a abrupt pipe expansion are widely carried out. The mechanism is not fully found evidently due to the instabilities of flow in a sudden change of the shape and appearance of turbulent shear layers in a recirculation region and secondary vortex near the corner. The purpose of this study is to obtain data through an experimental study of the swirling flow and heat transfer downstream of an abrupt expansion in a circular pipe with uniform heat flux. Experiments were carried out for the turbulent flow nd heat transfer downstream of an abrupt circular pipe expansion. The uniform heat flux condition was imposed to the downstream of the abrupt expansion by using an electrically heated pipe. Experimental data are presented for local heat transfer rates and local axial velocities in the tube downstream of an abrupt 3:1 & 2:1 expansion. Air was used as the working fluid in the upstream tube, the Reynolds number was varied from 60, 00 to 120, 000 and the swirl number range (based on the swirl chamber geometry, i.e. L/d ratio) in which the experiments were conducted were L/d=0, 8 and 16. Axial velocity increased rapidly at r/R=0.35 in the abrupt concentric expansion turbulent flow through the test tube in unswirled flow. It showed that with increasing axial distance the highest axial velocities move toward the tube wall in the case of the swirling flow abrupt expansion. A uniform wall heat flux boundary condition was employed, which resulted in wall-to-bulk temperatures ranging from 24.deg. C to 71.deg. C. In swirling flow, the wall temperature showed a greater increase at L/d=16 than any other L/d. The bulk temperature showed a minimum value at the pipe inlet, it also exhibited a linear increase with axial distance along the pipe. As swirl intensity increased, the location of peak Nu numbers was observed to shift from 4 to 1 step heights downstream of the expansion. This upstream movement of the maximum Nusselt number was accompanied by an increase in its magnitude from 2.2 to 8.8 times larger than fully developed tube flow values.