• Water for future
• /
• v.29 no.2
• /
• pp.199-207
• /
• 1996

To evaluate the usage of $\kappa-\imath$ turbulence closure for the analysis of thermal discharge behavior, a two-dimensional depth-integrated numerical model is developed. The developed model is applied to a steady flow in an open channel with simle geometry and the numerical results agree well with existing experimental data. The adequate simulation of recirculation, reattachment, and excess temperature rise at downstream of the outlet in the channel attributes to the correct calculation of turbulent eddy viscosity and diffusivity by $\kappa-\imath$ turbulence model. For an accurate prediction of thermal discharge behavior, the introduction of buoyancy production term, the modification of source/sink, and the correct input of turbulence constants of the $\kappa-\imath$ turbulence model are required.

• Nuclear Engineering and Technology
• /
• v.41 no.10
• /
• pp.1347-1360
• /
• 2009

A multi-dimensional component for the thermal-hydraulic system analysis code, MARS, was developed for a more realistic three-dimensional analysis of nuclear systems. A three-dimensional and two-fluid model for a two-phase flow in Cartesian and cylindrical coordinates was employed. The governing equations and physical constitutive relationships were extended from those of a one-dimensional version. The numerical solution method adopted a semi-implicit and finite-difference method based on a staggered-grid mesh and a donor-cell scheme. The relevant length scale was very coarse compared to commercial computational fluid dynamics tools. Thus a simple Prandtl's mixing length turbulence model was applied to interpret the turbulent induced momentum and energy diffusivity. Non drag interfacial forces were not considered as in the general nuclear system codes. Several conceptual cases with analytic solutions were chosen and analyzed to assess the fundamental terms. RPI air-water and UPTF 7 tests were simulated and compared to the experimental data. The simulation results for the RPI air-water two-phase flow experiment showed good agreement with the measured void fraction. The simulation results for the UPTF downcomer test 7 were compared to the experiment data and the results from other multi-dimensional system codes for the ECC delivery flow.

• Journal of Climate Change Research
• /
• v.7 no.3
• /
• pp.357-371
• /
• 2016

Buildings in the city acts as a cause of distorted wind direction, wind speed, causing the stagnation of the air flow. In the recent trend of climate change can not but consider the temperature rise of the urbanization. This study was aimed to analyze the thermal comfort of planetary boundary layer in different artificial constructions areas which has a direct impact on urban climate, and estimating the warming phenomena. Envi-met model was used to consider the urban structure associated with urban growth in order to precisely determine the impact of the building on the city weather condition. The analyzed values of thermal comfort index were temperature, wind speed, horizontal and vertical turbulent diffusivity. In particular, analysis of the PPD(Predicted Percentage of Dissatisfied) represents the human thermal comfort. In this study, by adjusting the arrangement and proportion of the top floor building in the urban it was found that the inflow of the fresh air and cooling can be derived low PPD. Vertical heat flux amount of the city caused by climate change was a factor to form a high potential temperature in the city and the accumulation of cold air does not appear near the surface. Based on this, to make the city effectively respond to climate change may require a long-term restructuring of urban spatial structure and density management.

• Journal of Environmental Science International
• /
• v.1 no.2
• /
• pp.105.2-117
• /
• 1992

The formation mechanism of the vortex streets in the lee of the mountain Is Investigated by a three-dimensional numerical model. The model is based upon the hydrostatic Boussinesq equations in which the vertical turbulent momentum flux is estimated by a turbulence parameterization scheme, but the horizontal viscosity is assumed to be constant. The results show that Karman vortex streets can form even without surface friction in a constant ambient flow with uniform stratification. The vortex formation is related to breaking of the mountain wave, which depends on the Froude number (Fr). In the case of a three-dimensional bell-shaped mountain, the wave breaking occurs when Fr is less than about 0.8, while a barman vortex forms when Fr is less than about 0.22. Vortex formation also depends on Reynolds number, which is estimated from the horizontal diffusivity. The vortex formation can be explained by the wave saturation theory given by Lindzen (1981) with some modification. Simulations in this study show that in the case of Karman vortex formation the momentum flux in the lower level is much larger than the saturated momentum flux, whereas it is almost equal to the saturated momentum at the upper levels as expected from the saturation theory. As a result, large flux divergence is produced in the lower layer, the mean flow is decelerated behind the mountain, and the horizontal wind shear forms between unmodified ambient wind. The momentum exchange between the mean flow and the mountain wave is produced by the turbulence within a breaking wave. From the result, well developed vortices like Karman vortex can be formed. . The results of the momentum budget calculated by the hydrostatic model are almost the same as nonhydrostatic results as long as horizontal scale of the mountain is 10 km. A well developed barman vortex similar to the hydrostatic one was simulated in the nonhydrostatic case. Therefore, we conclude that the hydrostatic assumption is adequate to investigate the origin of the Km8n vortex from the viewpoint of wave breaking.

• Environmental Sciences Bulletin of The Korean Environmental Sciences Society
• /
• v.1 no.2
• /
• pp.105-117
• /
• 1997

The formation mechanism of the vortex streets in the lee of the mountain is investigated by a three-dimensional numerical model. The model is based upon the hydrostatic Boussinesq equations in which the vertical turbulent momentum flux is estimated by a turbulence parameterization scheme, but the horizontal viscosity is assumed to be constant. The results show that Karman vortex streets can form even without surface friction in a constant ambient flow with uniform stratification. The vortex formation is related to breaking of the mountain wave, which depends on the Froude number (Fr). In the case of a three-dimensional bell-shaped mountain, the wave breaking occurs when Fr is less than about 0.8, while a Karman vortex forms when Fr is less than about 0.22. Vortex formation also depends on Reynolds number, which is estimated from the horizontal diffusivity. The vortex formation can be explained by the wave saturation theory given by Lindzen (1981) with some modification. Simulations in this study show that in the case of Karman vortex formation the momentum flux in the lower level is much larger than the saturated momentum flux whereas it is almost equal to the saturated momentum at the upper levels as expected from the saturation theory. As a result, large flux divergence is produced in the lower layer, the mean flow is decelerated behind the mountain, and the horizontal wind shear forms between unmodified ambient wind. The momentum exchange between the mean flow and the mountain wave is produced by the turbulence within a breaking wave. From the result, well developed vortices like Karman vortex can be formed. The results of the momentum budget calculated by the hydrostatic model are almost the same as nonhydrostatic results as long as horizontal scale of the mountain is 10 km. A well developed Karman vortex similar to the hydrostatic one was simulated in the nonhydrostatic case. Therefore, we conclude that the hydrostatic assumption is adequate to investigate the origin of the Karman vortex from the viewpoint of wave breaking.

• International Union of Geodesy and Geophysics Korean Journal of Geophysical Research
• /
• v.21 no.1
• /
• pp.17-30
• /
• 1993

The basic processes responsible for the generation of oceanic fronts were reviewed. In particular the process of a shelf sea front produced by tidal stirring was identified from the one dimensional model of the water column in the coastal area, which incorporates the microscale process for the formation of a tidal front. Also a new criterion to predict its location was suggested. The time evolutions of the distributions of density and turbulent kinetic energy calculated from the model show that the criterion for the formation of a thermocline can be predicted as $R{\delta}^4$~ constant for large $\delta$ ($\delta$＞0.5), but the dependence on $\delta$ decreases as $\delta$ goes to O, where $R=H^4Q/{K_b}^3$,{\;}{\delta}=1-Do/H$, Q is the buoyancy flux at the surface, $K_b$ is the eddy diffusivity maintained at the bottom and Do is the depth of a thermocline in the absence of bottom mixing. The depth of a thermocline was found to decrease as the bottom mixing increases for a given value of Do. The results were interpreted in comparison with the previous 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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.5
• /
• pp.1032-1043
• /
• 1989

The heat transfer characteristics of the drag reducing polymer solutions are investigated experimentally in the thermal entrance region of circular tube flows. Fluids used in experiments are the aqueous solutions of high molecular polymer, polyacrylamide Separan AP-273 and the range of polymer concentrations is from 20 to 1000 wppm. Two stainless steel tubes with inside diameter 8.5mm(L/D=712) and 10.3mm(L/D=1160) are used for the heat transfer flow loops. The flow loop is set up to measure friction factors and heat transfer coefficients of test sections in two different modes; the recirculating flow system and once-through flow system. The test tubes are heated directly by electricity to apply the constant heat flux boundary conditions to the wall. Three different types of adaptors are used to observe the effects of the upstream flow conditions of the heat transfer test sections. The viscosity and characteristic relaxation time of the test fluids circulating in the flow system are measured by the capillary tube viscometer and falling ball viscometer at regular time intervals. The installed adaptors exhibit slight effect on the entrance heat transfer of Newtonian fluid. However, no noticeable effects are observed for the entrance heat transfer of the drag reducing fluids. The order of magnitude of the thermal entrance lengths of the drag reducing fluids which follow the minimum friction asymptote is much longer than that of Newtonian fluids in turbulent flows. A new dimensionless parameter, the viscoelastic Graetz number, is defined and all the experimental data are recasted in terms of the viscoelastic Graetz number. The local Nusselt number of the viscoelastic fluids is represented as a function of flow behavior index n and the viscoelastic Graetz number. As degradation continues the viscosity and the characteristic relaxation time of the testing fluids decrease. Weissenberg number defined by the relaxation time and D/V appears to be a proper dimensionless parameter in describing degradation effects on heat transfer of the viscoelastic fluids.