• Korean Journal of Computational Design and Engineering
• /
• v.18 no.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.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2017.05a
• /
• pp.803-810
• /
• 2017

the liquid oxygen transitions to a supercritical state, causing rapid changes in properties and pseudo boiling in supercritical combustion. the combustion reaction operating in a supercritical state depends on the turbulence diffusion caused by difference of density, therefore, a study of the diffusion flow and pseudo boiling is required. Many researchers have studied these phenomena in the supercritical combustion, but A case study by various variables is inadequate. In this study, the flow field and flame structure were investigated numerically by changing the recirculation flow and liquid oxygen core length through oxygen-fuel ratio(O/F), combustor diameter and recess ratio at supercritical pressure condition.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.23 no.3
• /
• pp.18-27
• /
• 2019

The flow mixing characteristics for the heated dual jets were numerically studied by using URANS (unsteady Reynolds-averaged Navier-Stokes). The increased turbulent diffusion was obtained for the compressible flow, and the thermal diffusion of incompressible flow increased more than that of compressible flow. From the results of FFT and phase portraits, periodic and quasi-periodic states were observed as the jet spacing increased. It was observed that linear variations of merging points and combined points were different because unsteady flow determined the flow mixing characteristics for a large jet spacing.

• Journal of Radiation Protection and Research
• /
• v.24 no.1
• /
• pp.31-38
• /
• 1999

A real-time dose assessment system(FADAS : Following Accident Dose Assessment System) has been developed for the real-time accident consequence assessment against a nuclear accident. Field tracer experiment near Younggwang nuclear power plant was performed to improve the accuracy of developed system and to parameterize the site-specific parameters into the FADAS. The mean values and turbulent components of wind profile obtained through field experiment have been reflected to FADAS with site-specific conditions. The simulated results of diffusion model agreed well with the measured data through tracer experiment. The developed system is being used as a basic module of emergency preparedness system in Korea. The diffusion model which can be reflected site-specific parameters will be improved through field experiments continuously.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.34 no.5
• /
• pp.505-513
• /
• 2010

The tribrachial flame has attracted interest as a basic structure of the flame edge. This flame structure helps understand stabilization of laminar flames and re-ignition of turbulent flames. A number of analytical and experimental studies have been carried out on the tribrachial flame. However, the effect of the variation of the flammability limits on the structure of the tribrachial flame has not been studied in detail. In this study, the effect of non-symmetric flammability limits on the flame structure was investigated by adopting a simple numerical scheme based on several laminar flame theories. A fixed velocity field was considered and boundary matching algorithm was used on the premixed branch. The variation of the diffusion branches under the non-symmetric flammability limits and heat loss was investigated. The formation and extinction of the diffusion branch behind the premixed branch were successfully described. This basic study can help understand the fundamental structure of the flame and can form the basis of subsequent detailed studies.

• The Bulletin of The Korean Astronomical Society
• /
• v.46 no.2
• /
• pp.47.3-48
• /
• 2021

We studied the large scale dynamo process in a system forced by helical magnetic field. The dynamo process is basically nonlinear, but can be linearized with 𝛼&𝛽 coefficients and large scale magnetic field $\bar{B}$. This is very useful to the investigation of solar (stellar) dynamo. A coupled semi-analytic equations based on statistical mechanics are used to investigate the exact evolution of 𝛼&𝛽. This equation set needs only magnetic helicity ${\bar{H}}_M({\equiv}{\langle}{\bar{A}}{\cdot}{\bar{B}}{\rangle},\;{\bar{B}}={\nabla}{\times}{\bar{A}})$ and magnetic energy ${\bar{E}}_M({\equiv}{\langle}{\bar{B}}^2{\rangle}/2)$. They are fundamental physics quantities that can be obtained from the dynamo simulation or observation without any artificial modification or assumption. 𝛼 effect is thought to be related to magnetic field amplification. However, in reality the averaged 𝛼 effect decreases very quickly without a significant contribution to ${\bar{B}}$ field amplification. Conversely, 𝛽 effect contributing to the magnetic diffusion maintains a negative value, which plays a key role in the amplification with Laplacian ∇²(= - k²) for the large scale regime. In addition, negative magnetic diffusion accounts for the attenuation of plasma kinetic energy E_V(= 〈 U² 〉/2) (U: plasma velocity) when the system is saturated. The negative magnetic diffusion is from the interaction of advective term - U • ∇ B from magnetic induction equation and the helical velocity field. In more detail, when 'U' is divided into the poloidal component U_pol and toroidal one U_tor in the absence of reflection symmetry, they interact with - B • ∇ U and - U • ∇ B from ∇ × 〈 U × B 〉 leading to 𝛼 effect and (negative) 𝛽 effect, respectively. We discussed this process using the theoretical method and intuitive field structure model supported by the simulation result.

• Korea-Australia Rheology Journal
• /
• v.17 no.3
• /
• pp.131-140
• /
• 2005

In this work we show the results of our most recent Direct Numerical Simulations (DNS) of turbulent viscoelastic channel flow using spectral spatial approximations and a stabilizing artificial diffusion in the viscoelastic constitutive model. The Finite-Elasticity Non-Linear Elastic Dumbbell model with the Peterlin approximation (FENE-P) is used to represent the effect of polymer molecules in solution, The corresponding rheological parameters are chosen so that to get closer to the conditions corresponding to maximum drag reduction: A high extensibility parameter (60) and a moderate solvent viscosity ratio (0.8) are used with two different friction Weissenberg numbers (50 and 100). We then first find that the corresponding achieved drag reduction, in the range of friction Reynolds numbers used in this work (180-590), is insensitive to the Reynolds number (in accordance to previous work). The obtained drag reduction is at the level of $49\%\;and\;63\%$, for the friction Weissenberg numbers 50 and 100, respectively. The largest value is substantially higher than any of our previous simulations, performed at more moderate levels of viscoelasticity (i.e. higher viscosity ratio and smaller extensibility parameter values). Therefore, the maximum extensional viscosity exhibited by the modeled system and the friction Weissenberg number can still be considered as the dominant factors determining the levels of drag reduction. These can reach high values, even for of dilute polymer solution (the system modeled by the FENE-P model), provided the flow viscoelasticity is high, corresponding to a high polymer molecular weight (which translates to a high extensibility parameter) and a high friction Weissenberg number. Based on that and the changes observed in the turbulent structure and in the most prevalent statistics, as presented in this work, we can still rationalize for an increasing extensional resistance-based drag reduction mechanism as the most prevalent mechanism for drag reduction, the same one evidenced in our previous work: As the polymer elasticity increases, so does the resistance offered to extensional deformation. That, in turn, changes the structure of the most energy-containing turbulent eddies (they become wider, more well correlated, and weaker in intensity) so that they become less efficient in transferring momentum, thus leading to drag reduction. Such a continuum, rheology-based, mechanism has first been proposed in the early 70s independently by Metzner and Lamley and is to be contrasted against any molecularly based explanations.

• Nuclear Engineering and Technology
• /
• v.55 no.5
• /
• pp.1802-1813
• /
• 2023

Conjugate heat transfer between liquid metal and solid is a common phenomenon in a liquid-metal-cooled fast reactor's fuel assembly and heat exchanger, dramatically affecting the reactor's safety and economy. Therefore, comprehensively studying the sophisticated conjugate heat transfer in a liquid-metal-cooled fast reactor is profound. However, it has been evidenced that the traditional Simple Gradient Diffusion Hypothesis (SGDH), assuming a constant turbulent Prandtl number (Pr_t,, usually 0.85 - 1.0), is inappropriate in the Computational Fluid Dynamics (CFD) simulations of liquid metal. In recent decades, numerous studies have been performed on the four-equation model, which is expected to improve the precision of liquid metal's CFD simulations but has not been introduced into the conjugate heat transfer calculation between liquid metal and solid. Consequently, a four-equation model, consisting of the Abe k - ε turbulence model and the Manservisi k_𝜃 - ε_𝜃 heat transfer model, is applied to study the conjugate heat transfer concerning liquid metal in the present work. To verify the numerical validity of the four-equation model used in the conjugate heat transfer simulations, we reproduce Johnson's experiments of the liquid lead-bismuth-cooled turbulent pipe flow using the four-equation model and the traditional SGDH model. The simulation results obtained with different models are compared with the available experimental data, revealing that the relative errors of the local Nusselt number and mean heat transfer coefficient obtained with the four-equation model are considerably reduced compared with the SGDH model. Then, the thermal-hydraulic characteristics of liquid metal turbulent pipe flow obtained with the four-equation model are analyzed. Moreover, the impact of the turbulence model used in the four-equation model on overall simulation performance is investigated. At last, the effectiveness of the four-equation model in the CFD simulations of liquid sodium conjugate heat transfer is assessed. This paper mainly proves that it is feasible to use the four-equation model in the study of liquid metal conjugate heat transfer and provides a reference for the research of conjugate heat transfer in a liquid-metal-cooled fast reactor.

• The Magazine of the Society of Air-Conditioning and Refrigerating Engineers of Korea
• /
• v.17 no.3
• /
• pp.238-248
• /
• 1988

The flow pattern inside the power condenser is generally known to be very complicated due to the phase change and turbulence effects as well as the effect of condenser geometry. In the present study, the flow pattern inside the power condenser was numerically simulated with a personal computer. The widely known CHAMPION 2/E/FIX(Concentration, Heat and Momentum Program Instruction Outfit, 2D/Elliptic/Fixed grid) computer code was modified for this purpose. The flow was asssumed to be two-dimensional and steady-state, and the tube bank was considered to be homogeneous porous medium. Simple turbulent diffusion coefficients based on the appropriate experiments were obtained for the computation. Through this analytical approach, the flow pattern could be predicted fairly well. The computational results also show that the location of the air vent plays an important key role in determining the efficiency of the condenser.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.10a
• /
• pp.129-134
• /
• 2000

An efficient numerical method to solve the unsteady incompressible Navier-Stokes equations is developed. A fully implicit time advancement is employed to avoid the CFL(Courant-Friedrichs-Lewy) restriction, where the Crank-Nicholson discretization is used for both the diffusion and convection terms. Based on a block LU decomposition, velocity-pressure decoupling is achieved in conjunction with the approximate factorization. Main emphasis is placed on the additional decoupling of the intermediate velocity components with only n th time step velocity The temporal second-order accuracy is Preserved with the approximate factorization without any modification of boundary conditions. Since the decoupled momentum equations are solved without iteration, the computational time is reduced significantly. The present decoupling method is validated by solving the turbulent minimal channel flow unit.