• Journal of Mechanical Science and Technology
• /
• 제16권2호
• /
• pp.262-269
• /
• 2002

The effects of velocity boundary conditions on the structure of methane-air nonpremixed counterflow flames were investigated by two-dimensional numerical simulation. Two low global strain rates, 12 s$\^$-1/ and 20 s$\^$-1/, were considered for comparison with measurements. Buoyancy was conformed to have strong effects on the flame structure at a low global strain rate. It was shown that the location where a top hat velocity profile was imposed is sensitive to the flame structure, and that the computed temperature along the centerline agrees well with the measurements when plug flow was imposed at the inner surface of the screen nearest the duct exit.

• 한국철도학회:학술대회논문집
• /
• 한국철도학회 2008년도 춘계학술대회 논문집
• /
• pp.1740-1746
• /
• 2008

The present study performed numerical investigation to analyze the smoke behavior in the rescue station by using the commercial CFD code (FLUENT Ver 6.3). The present study adopted a 10MW ultrafast mode for simulation, and it also used the MVHS(Modify Volumetric Heat Source) model modified from the original VHS(Volumetric Heat Source) model in order to treat the product generation and the oxygen consumption under the stoichiometric state. In addition, the present simulation includes the species conservation equation for the materialization of heat source and the estimation of smoke movement. From the results, the smoke flows are moving along the ceiling because of thermal buoyancy force and as time goes, the smoke gradually moves downward at the vicinity of the entrance. Moreover, without using ventilation, it is found that the smoke flows no longer spread across the cross-passages because the pressure in the non-accident tunnel is higher than that in the accident tunnel.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2004년도 춘계학술대회
• /
• pp.1909-1914
• /
• 2004

An immersed boundary method for simulation of density-stratified flows is developed and applied to computation of viscous flows over two-dimensional obstacles in a bounded domain under stable density stratification. Density sources/sinks are introduced on the body surface. Two obstacle shapes are used, a vertical barrier and a smooth cosine-shaped hill; weak stratification, defined by $K=ND/{\pi}U{\leq}1$, where U, N, and D are the upstream velocity, buoyancy frequency, and domain height, respectively, is considered. The results are consistent with other authors' calculations, and shed light on computation of density-stratified flows in complex geometries.

• Nuclear Engineering and Technology
• /
• 제41권9호
• /
• pp.1191-1204
• /
• 2009

In this study, a CFD methodology is employed to address the problem of the prediction of the flow in a T-junction. An Unsteady Reynolds Averaged Navier-Stokes (URANS) approach has been selected for its low computational cost. Moreover, Unsteady Reynolds Navier-Stokes methodologies do not need complex boundary formulations for the inlet and the outlet such as those required when using Large Eddy Simulation (LES) or Direct Numerical Simulation (DNS). The results are compared with experimental data and an LES calculation. In the past, URANS has been tried on T-junctions with mixed results. The biggest limit observed was the underestimation of the oscillatory behavior of the temperature. In the present work, we propose a comprehensive approach able to correctly reproduce the root mean square (RMS) of the temperature directly downstream of the T-junction for cases where buoyancy is not present.

• Nuclear Engineering and Technology
• /
• 제54권9호
• /
• pp.3567-3579
• /
• 2022

This paper provides an assessment of fluid transport and mixing processes inside the primary circuit of the test facility ROCOM through the numerical simulation of Test 2.1 with the system code ATHLET. The experiment represents an asymmetric injection of cold and non-borated water into the reactor coolant system (RCS) of a pressurized water reactor (PWR) to restore core cooling, an emergency procedure which may subsequently trigger a core re-criticality. The injection takes place at low velocity under single-phase subcooled conditions and presents a major challenge for the simulation in lumped parameter codes, due to multidimensional effects in horizontal piping and vessel arising from density gradients and gravity forces. Aiming at further validating ATHLET 3-D capabilities against horizontal geometries, the experiment conditions are applied to a ROCOM model, which includes a newly developed horizontal pipe object to enhance code prediction inside coolant loops. The obtained results show code strong simulation capabilities to represent multidimensional flows. Enhanced prediction is observed at the vessel inlet compared to traditional 1-D approach, whereas mixing overprediction from the descending denser plume is observed at the upper-half downcomer region, which leads to eventual deviations at the core inlet.

• 한국해양공학회지
• /
• 제31권6호
• /
• pp.388-396
• /
• 2017

A Lagrangian approach based computational fluid dynamics (CFD) was used to simulate large and/or sharp deformations and fragmentations of interfaces, including free surfaces, through tracing each particle with physical quantities. According to the concept of the particle-based CFD method, it is possible to apply it to both fluid particles and solid particles such as sand, gravel, and rock. However, the presence of more than two different phases in the same domain can make it complicated to calculate the interaction between different phases. In order to solve multiphase problems, particle interaction models for multiphase problems, including surface tension, buoyancy-correction, and interface boundary condition models, were newly adopted into the moving particle semi-implicit (MPS) method. The newly developed MPS method was used to simulate a typical validation problem involving dam breaking. Because the soil and other particles, excluding the water, may have different viscosities, various viscosity coefficients were applied in the simulations for validation. The newly developed and validated MPS method was used to simulate the mobile beds induced by broken dam flows. The effects of the viscosity on soil particles were also investigated.

• 한국해양공학회지
• /
• 제26권1호
• /
• pp.17-26
• /
• 2012

A coupled analysis of a floating crane barge with a crane wire and hanging structure is carried out in thetime domain. The motion analysis of the crane barge is based on the floating multi-body dynamics, and thecrane wire is modeled as a simple spring tension. The hanging structure is assumed to be a rigid body with 3 degree-of-freedom translational motion. In this study, numerical simulations were conducted at three different stages. First, the developed code was validated by comparing the time-domain motion response of a crane barge with the frequency-domain results. Then, a coupled analysis of a crane barge and simple structure hanging by the crane wire was performed using the present scheme. The motion response and wire tension from the present calculations are compared with the results of OrcaFlex. The agreement between the two sets of results isfairly good. Last, lowering simulations in regular and irregular waves were conducted considering buoyancy changes in the hanging structure. The effects of the wave conditions, structure's weight, wire length, and lowering speed on the wire tension are considered.

• 대한기계학회논문집
• /
• 제18권10호
• /
• pp.2712-2723
• /
• 1994

Numerical simulation of an axisymmetric ethylene-air jet diffusion flame has been carried out in order to investigate flame dynamics and soot formation. The model solves the time-dependent Navier-Stokes equations and includes models for soot formation, chemical reaction, molecular diffusion, thermal conduction, and radiation. Numerically FCT(Flux Corrected Transport) and DOM(Discrete Ordinate Method) methos are used for convection and radiation trasport respectively. Simulation was conducted for a 5 cm/sec fuel jet flowing into a coflowing air stream. The maximum flame temperature was found to be approximately 2100 K, and was located at an axial position of approximately 5 cm from the base of the flame. The maximum soot volume fraction was about $7{\times}10^{-7}$, and was located within the high temperature region where the fuel mole fraction ranges from 0.01 to 0.1. The buoyancy-driven low-frequency(12~13 Hz) structures convected along the outer region of the flame were captured. In case without radiation trasport, the maximum temperature was higher by 150 K than in case with radiation. Also the maximum soot volume fraction reached about $8{\times}10^{-6}$. As the the hydrocarbon fuel forms many soot particles, the radiation transport becomes to play a more important role.

• 한국CDE학회논문집
• /
• 제18권6호
• /
• pp.399-406
• /
• 2013

This paper describes a new simulation technique for advection-diffusion phenomena over the sea surface using the lattice Boltzmann method (LBM), capable of predicting oil dispersion from tankers. The LBM is used to solve the pollutant transport problem within the framework of the ocean environment. The sea space is represented by the lattices, where each lattice has the information on oil transportation. Since dispersed oils (i.e., oil droplets) at sea are transported by convection due to waves, buoyancy, and turbulent diffusion, the conservation of mass and many physical oil transport rules were used in the prediction model. Since the LBM is modeled using the uniform lattices and simple rules, it can be easily accelerated by the parallel mechanism, for example, GPU-accelerated method. The proposed model using the LBM is used to simulate a simple pollution event with the oil pollutants of 10,000 kL. The simulation results indicate that the LBM method accelerated with the GPU is 6 times faster than that without the GPU.

• 한국수자원학회:학술대회논문집
• /
• 한국수자원학회 2019년도 학술발표회
• /
• pp.98-98
• /
• 2019

중력류 또는 밀도류는 주변 유체에 비해 상대적으로 밀도가 큰 유체가 밀도차에 의한 추진력으로 흐르는 것이다. 중력류의 수치모델링에는 두 가지 어려움이 있다. 즉, 적합한 지배방정식을 구성하여 적용하는 것 그리고 난류의 영향을 합리적으로 반영하는 것이다. 기존 중력류 해석을 위한 지배방정식들은 유체의 연속방정식과 운동량 방정식 그리고 밀도 또는 농도의 이송방정식을 조합하여 구성된다. 이들 지배방정식을 이용한 연구들은 대부분 두 유체 사이의 밀도차가 충분히 작아서 밀도 변동(variations)의 영향은 오로지 부력항에서만 유지된다는 Boussinesq 근사에 근거를 둔다. 그리고 이송방정식에서 밀도 또는 농도의 확산계수을 점성계수의 함수로 표현하기 위해서 Schmidt 수를 이용한다. 수치모델링에서 Schimdt 수는 상수값을 적용하지만, 이 값은 밀도의 연직방향 경사에 근거한 부력빈도(buoyancy frequency)와 난류량의 따라 큰 차이를 보이는 것으로 알려져있다. 한편, 표준 통계학적 난류모델과 벽함수를 적용한 수치모델링은 초기 중력에 의해서 무너지는(slumping) 단계를 넘어 관성력으로 추진되는 단계와 점성 효과가 지배적인 단계에서는 정확도에 현저히 낮아지기 때문에 대부분 큰와모의(large-eddy simulation, LES) 또는 DNS(direct numerical simulation)수준의 고해상도(high-resolution) 해석기법을 적용하여 공학적인 문제에 적용하는 데는 한계가 있다. 이 연구에서는 Boussinesq 근사와 Schmidt 수를 사용하지 않으며, LES 보다 적용이 용이한 DES (detached-eddy simulation)기법을 조합한 다상흐름 수치모델을 적용하여 중력류를 해석을 시도하였다. 수치해석결과를 실험값과 함께 기존 수치모델링 기법으로 구한 수치해와 비교분석하여 이 연구에서 개발 및 적용된 수치모델링 기법의 적용성을 평가한다.