• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.10
• /
• pp.1401-1411
• /
• 2003

A two-dimensional direct numerical simulation was performed to investigate the flame structure of CH$_4$$N_2$-air counterflow nonpremixed flame interacting with a single vortex. The detailed transport properties and a modified 16-step augmented reduced mechanism based on Miller and Bowman's detailed chemistry were adopted in this simulation. The characteristic vortex and chemical time scales were introduced to quantify and investigate the extinction phenomenon during a flame-vortex interaction. The results showed that fuel- and air-side vortex cause an unsteady extinction. In this case, the flame interacting with a vortex was extinguished at much larger scalar dissipation rate than steady flame. It was also found that the air-side vortex extinguished a flame more rapidly than the fuel-side vortex. Furthermore, it was noted that the degree of unsteady effect experienced by a flame can be investigated by comparing the above two characteristic time scales, and this analysis could give an appropriate reason for the results of the previously reported experiment.

• Journal of computational fluids engineering
• /
• v.3 no.1
• /
• pp.11-21
• /
• 1998

A vorticity-based method for the numerical solution of the two-dimensional incompressible Navier-Stokes equations is presented. The governing equations for vorticity, velocity and pressure variables are expressed in an integro-differential form. The global coupling between the vorticity and the pressure boundary conditions is fully considered in an iterative procedure when numerical schemes are employed. The finite volume method of the second order TVD scheme is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition. The velocity field is obtained by using the Biot-Savart integral. The Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well established for potential flow analysis. The present formulation is validated by comparison with data from the literature for the two-dimensional cavity flow driven by shear in a square cavity. We take two types of the cavity now: (ⅰ) driven by non-uniform shear on top lid and body forces for which the exact solution exists, and (ⅱ) driven only by uniform shear (of the classical type).

• 한국연소학회:학술대회논문집
• /
• 2002.06a
• /
• pp.52-60
• /
• 2002

A two-dimensional direct numerical simulation is performed to investigate the flame structure of $CH_4/N_2-Air$ counterflow nonpremixed flame interacting with a single vortex. The detailed transport properties and a modified 16-step augmented reduced mechanism based on Miller and Bowman's detailed reaction mechanism are adopted in this calculation. To quantify the strain on flame induced by a vortex, a scalar dissipation rate (SDR) is introduced. Results show that the fuel and air-side vortex cause an unsteady extinction. In this case, the flame interacting with a vortex is extinguished at much larger SDR than steady flame. It is also found that air-side vortex extinguishes a flame more rapidly than fuel-side vortex. The unsteady effect induced by flame-vortex interaction does not lead to a transient OH overshoot of the maximum steady concentration observed in experiment, while $HO_2$ radical increases more than the maximum steady concentration with increasing SDR. In addition, it is seen that NO and $NO_2$ are not sensitive to the unsteady variation of SDR.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.13 no.1
• /
• pp.68-75
• /
• 2005

By utilizing a semi-empirical soot model, the applicability of the laminar flamelet concept fur simulating the formation and oxidation of soot in the laminar diffusion flame has been studied. The source terms for two transport equations of the soot formation and oxidation are calculated in the mixture fraction/scalar dissipation rate space for laminar flamelets and stored in a library. In this study, emphasis is given to the interaction associated with radiation and soot formation. The radiative heat loss is obtained by solving the radiative transfer equation using the unstructured grid finite volume method with the WSGGM. The calculated temperatures and soot volume fractions agree relatively well with the experimental data and the previous numerical results of Kaplan et al. using the detailed chemistry.

• Journal of the Korean Society of Combustion
• /
• v.9 no.2
• /
• pp.10-17
• /
• 2004

A two-dimensional direct numerical simulation is performed to investigate the flame structure of $CH_4/N_2-Air$ counterflow nonpremixed flame interacting with a single vortex. The detailed transport properties and a modified 16-step augmented reduced mechanism based on Miller and Bowman#s detailed reaction mechanism are adopted in this calculation. To quantify the strain on flame induced by a vortex, a scalar dissipation rate (SDR) is introduced. The results show that fuel-side and air-side vortex cause an unsteady extinction. In this case, the flame interacting with a vortex is extinguished at much larger SDR than steady flame. It is also found that air-side vortex extinguishes a flame more rapidly than fuel-side vortex. The unsteady effect induced by flame-vortex interaction does not lead to a transient OH overshoot of the maximum steady concentration observed in experiment, while $HO_2$ radical increases more than the maximum steady concentration with increasing SDR. In addition, it is seen that NO and $NO_2$ are not sensitive to the unsteady variation of SDR.

• Nuclear Engineering and Technology
• /
• v.55 no.4
• /
• pp.1428-1438
• /
• 2023

A novel Fission Diffusion Synthetic Acceleration (FDSA) method is developed and implemented as a part of a hybrid neutronics method for source convergence acceleration and variance reduction in Monte Carlo (MC) criticality calculations. The acceleration of the MC calculation stems from constructing a synthetic operator and solving a low-order problem using information obtained from previous MC calculations. By applying the P₁ approximation, two correction terms, one for the scalar flux and the other for the current, can be solved in the low-order problem and applied to the transport solution. A variety of one-dimensional (1-D) and two-dimensional (2-D) numerical tests are constructed to demonstrate the performance of FDSA in comparison with the standalone MC method and the coupled MC and Coarse Mesh Finite Difference (MC-CMFD) method on both intended purposes. The comparison results show that the acceleration by a factor of 3-10 can be expected for source convergence and the reduction in MC variance is comparable to CMFD in both slab and full core geometries, although the effectiveness of such hybrid methods is limited to systems with small dominance ratios.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.105-105
• /
• 2021

본 연구에서는 2차원 비선형 천수모형과 수심평균된 스칼라 이송모형을 해석하는 수치모형에 대해 기술하였다. 수치모형의 정확성을 보장함과 동시에 안정성을 높이기 위해 유한체적법, 플럭스 재구성 및 minmod 제한자를 사용하였다. 비선형 천수방정식의 이송항과 바닥 경사항은 계산된 수심의 양수 보존과 흐름의 정상 상태를 보장하기 위한 second order well-balanced positivity preserving central-upwind method를 이용하여 수치적으로 이산화되었다. 마찬가지로, 이송-확산 방정식 내 이송항은 동일한 2차 풍상차분법을 통해 수치적으로 풀이하였다. 격자점 경계면에서의 불연속으로 인한 수치진동을 방지하기 위해 이송항의 계산에 포함된 보존항의 차이로 인해 발생하는 스칼라의 수치확산을 최소화하기 위해 무차원의 비소산함수를 도입하였다. 또한, 확산항은 유한차분법을 이용하여 이산화하였다. 제안된 수치모형은 시간미분항의 계산을 위해 오일러 기법을 적용하여 계산된 수심 및 스칼라의 양수 보존여부와 함께 정지된 흐름의 정상 상태의 보존여부를 확인하였다. 제안된 수치모형의 해석 정확성을 평가하기 위해 1, 2차원 공간 내 다양한 흐름 조건에서의 해석해를 이용한 3개의 벤치마크 테스트를 수행하였다. 평균 제곱근 오차(Root Mean Squared Error, RMSE)를 산정하여 수치모형의 성능을 정량적으로 평가하였으며, 비소산함수를 적용함에 따라 스칼라의 수치확산이 감소하게 되었음을 확인하였다. 또한, 세 차례의 벤치마크 테스트 결과는 공통적으로 수치모형에 의해 계산된 결과값이 비소산함수를 고려함에 따라 해석해와 잘 일치함을 확인하였다.

• Asian Journal of Atmospheric Environment
• /
• v.6 no.3
• /
• pp.222-233
• /
• 2012

The nocturnal low-level jet makes a significant impact on carbon and water exchanges and turbulent mixing processes in the atmospheric boundary layer. This study reports a case study of nocturnal surface fluxes such as $CO_2$ and water vapor in the surface layer observed at a flat and homogeneous site in the presence of low-level jets (LLJs). In particular, it documents the temporal evolution of the overlying jets and the coincident response of surface fluxes. The present study highlights several factors linking the evolution of low-level jets to surface fluxes: 1) wavelet analysis shows that turbulent fluxes have similar time scales with temporal scale of LLJ evolution; 2) turbulent mixing is enhanced during the transition period of low-level jets; and 3) $CO_2$, water vapor and heat show dissimilarity from momentum during the period. We also found that LLJ activity is related not only to turbulent motions but also to the divergence of mean flow. An examination of scalar profiles and turbulence data reveal that LLJs transport $CO_2$ and water vapor by advection in the stable boundary layer, suggesting that surface fluxes obtained from the micrometeorological method such as nocturnal boundary layer budget technique should carefully interpreted in the presence of LLJs.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.15-28
• /
• 1998

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

• Environmental Engineering Research
• /
• v.17 no.2
• /
• pp.83-88
• /
• 2012

Currently, the technology of $CO_2$ capture and storage (CCS) has become the main issue for climate change and global warming. Among CCS technologies, the prediction of $CO_2$ behavior underground is very critical for $CO_2$ storage design, especially for its safety. Hence, the purpose of this paper is to model and simulate $CO_2$ flow and its heat transfer characteristics in a storage site, for more accurate evaluation of the safety for $CO_2$ storage process. In the present study, as part of the storage design, a micro pore-scale model was developed to mimic real porous structure, and computational fluid dynamics was applied to calculate the $CO_2$ flow and thermal fields in the micro pore-scale porous structure. Three different configurations of 3-dimensional (3D) micro-pore structures were developed, and compared. In particular, the technique of assigning random pore size in 3D porous media was considered. For the computation, physical conditions such as temperature and pressure were set up, equivalent to the underground condition at which the $CO_2$ fluid was injected. From the results, the characteristics of the flow and thermal fields of $CO_2$ were scrutinized, and the influence of the configuration of the micro-pore structure on the flow and scalar transport was investigated.