• Journal of computational fluids engineering
• /
• v.21 no.4
• /
• pp.1-10
• /
• 2016

The analysis of characteristics of turbulent flow and thermal boundary layer for natural convection caused by fire along vertical wall is performed. The 4m-high vertical copper plate is heated and kept at a uniform surface temperature of $60^{\circ}C$ and the surrounding fluid (air) is kept at $16.5^{\circ}C$. The flow and temperature is solved by large eddy simulation(LES) of FDS code(Ver.6), in which the viscous-sublayer flow is calculated by Werner-Wengle wall function. The whole analyzed domain is assumed as turbulent region to apply wall function even through the laminar flow is transient to the turbulent flow between $10^9$<$Gr_z$<$10^{10}$ in experiments. The various grids from $7{\times}7{\times}128$ to $18{\times}18{\times}128$ are applied to investigate the sensitivity of wall function to $x^+$ value in LES simulation. The mean velocity and temperature profiles in the turbulent boundary layer are compared with experimental data by Tsuji & Nagano and the results from other LES simulation in which the viscous-sublayer flow is directly solved with many grids. The relationship between heat transfer rate($Nu_z$) and $Gr_zPr$ is investigated and calculated heat transfer rates are compared with theoretical equation and experimental data.

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

The shell and tube heat exchangers installed in the propulsion system test complex (PSTC) at the Naro Space Center heats cryogenic helium to 500 K with a heat transfer oil. As the experimental helium outlet temperature was lower than expected (less than 100 K), the boundary layer effect of the heat transfer oil is predicted to be the cause of the performance deterioration. A computational fluid dynamics (CFD) analysis was performed to verify where the boundary layer effect exists; however, the boundary layer effect has no significant impact on the performance of the heat exchanger. An alternative method to improve the performance of the heat exchanger by changing the heat transfer oil has been discussed in this paper. The low viscosity and high thermal conductivity at high temperature (~500 K) of heat transfer oil at the shell-side are required to improve the thermal performance of the heat exchanger. The experimental performance of the heat exchanger, used to exchange heat between the cryogenic helium and hot heat transfer oil at the PSTC are summarized in this paper.

• Journal of the Korean Society of Visualization
• /
• v.18 no.2
• /
• pp.10-17
• /
• 2020

The curvature of wound boundaries has been identified as a key modulator that determines a type of force responsible for cell migration. While several studies report how certain curvatures of the boundary correlate with the rate at which the wound closes, it remains unclear how these curvatures are spatiotemporally formed to regulate the healing process. We investigated the dynamic changes in the boundary curvatures by visualizing cell migration patterns. Locally, cells at the convex boundary continuously move forward with transmitting kinetic responses behind to the cells away from the boundary, and cells at the concave boundary exhibit dramatic contracting motion, like a purse-string, when they accumulate enough negative curvatures to gain the thrust toward the void. Globally, the dynamics of boundary geometries are controlled by the diffusive flow of cells driven by the density gradient between the wound area and the cell layer.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1310-1317
• /
• 2017

GASFLOW-MPI is a widely used scalable computational fluid dynamics numerical tool to simulate the fluid turbulence behavior, combustion dynamics, and other related thermal-hydraulic phenomena in nuclear power plant containment. An efficient scalable linear solver for the large-scale pressure equation is one of the key issues to ensure the computational efficiency of GASFLOW-MPI. Several advanced Krylov subspace methods and scalable preconditioning methods are compared and analyzed to improve the computational performance. With the help of the powerful computational capability, the large eddy simulation turbulent model is used to resolve more detailed turbulent behaviors. A backward-facing step flow is performed to study the free shear layer, the recirculation region, and the boundary layer, which is widespread in many scientific and engineering applications. Numerical results are compared with the experimental data in the literature and the direct numerical simulation results by GASFLOW-MPI. Both time-averaged velocity profile and turbulent intensity are well consistent with the experimental data and direct numerical simulation result. Furthermore, the frequency spectrum is presented and a -5/3 energy decay is observed for a wide range of frequencies, satisfying the turbulent energy spectrum theory. Parallel scaling tests are also implemented on the KIT/IKET cluster and a linear scaling is realized for GASFLOW-MPI.

• Journal of Advanced Marine Engineering and Technology
• /
• v.20 no.3
• /
• pp.11-18
• /
• 1996

Liquid droplets of water and argon surrounded by their vapor have been simulated by the milecular dynamics method. To explore the surface phenomena of clusters, each molecule is classified into 'liquid', 'surface', or 'vapor' with respect to the number of neighbor molecules. The contribution of a 'surface' molecule of the water cluster to the far infrared spectrum is almist the same as that of the 'liquid' molecule. Hence, the liquid-vapor interface is viewed as geometrically and temporally varying boundary of 'liquid' molecules with only a single layer of 'surface' molecules that might have different characteristics from the 'liquid' molecules. The time scale of the 'phase change' of each molecule is estimated for the argon cluster by observing the instantancous kinetic and potential energies of each molecule. To compare the feature of clusters with macroscopic droplets, the temperature dependence of the surface tension of the argon cluster is estimated.

• Wind and Structures
• /
• v.35 no.1
• /
• pp.1-20
• /
• 2022

This research work presents an experimental study's outcomes to reveal the impact of an O-ring on the flow control over a sphere placed in a turbulent boundary layer. The investigation is performed quantitatively and qualitatively using particle image velocimetry (PIV) and dye visualization. The sphere model having a diamater of 42.5 mm is located in a turbulent boundary layer flow over a smooth plate for gap ratios of 0≤G/D≤1.5 at Reynolds number of 5 × 10³. Flow characteristics, including patterns of instantaneous vorticity, streaklines, time-averaged streamlines, velocity vectors, velocity fluctuations, Reynolds stress correlations, and turbulence kinetic energy (), are compared and discussed for a naked sphere and spheres having O-rings. The boundary layer velocity gradient and proximity of the sphere to the flat plate profoundly influence the flow dynamics. At proximity ratios of G/D=0.1 and 0.25, a wall jet is formed between lower side of the sphere and flat plate, and velocity fluctuations increase in regions close to the wall. At G/D=0.25, the jet flow also induces local flow separations on the flat plate. At higher proximity ratios, the velocity gradient of the boundary layer causes asymmetries in the mean flow characteristics and turbulence values in the wake region. It is observed that the O-ring with various placement angles (𝜃) on the sphere has a considerable alteration in the flow structure and turbulence statistics on the wake. At lower placement angles, where the O-ring is closer to the forward stagnation point of the sphere, the flow control performance of the O-ring is limited; however, its impact on the flow separation becomes pronounced as it is moved away from the forward stagnation point. At G/D=1.50 for O-ring diameters of 4.7 (2 mm) and 7 (3 mm) percent of the sphere diameter, the -ring exhibits remarkable flow control at 𝜃=50° and 𝜃=55° before laminar flow separation occurrence on the sphere surface, respectively. This conclusion is yielded from narrowed wakes and reductions in turbulence statistics compared to the naked sphere model. The O-ring with a diameter of 3 mm and placement angle of 50° exhibits the most effective flow control. It decreases, in sequence, streamwise velocity fluctuations and length of wake recovery region by 45% and 40%, respectively, which can be evaluated as source of decrement in drag force.

• Advances in aircraft and spacecraft science
• /
• v.9 no.5
• /
• pp.415-431
• /
• 2022

Secondary flows have a huge impact on losses generation in modern low pressure gas turbines (LPTs). At design point, the interaction of the blade profile with the end-wall boundary layer is responsible for up to 40% of total losses. Therefore, predicting accurately the end-wall flow field in a LPT is extremely important in the industrial design phase. Since the inlet boundary layer profile is one of the factors which most affects the evolution of secondary flows, the first main objective of the present work is to investigate the impact of two different inlet conditions on the end-wall flow field of the T106A, a well known LPT cascade. The first condition, labeled in the paper as C1, is represented by uniform conditions at the inlet plane and the second, C2, by a flow characterized by a defined inlet boundary layer profile. The code used for the simulations is based on the Discontinuous Galerkin (DG) formulation and solves the Reynolds-averaged Navier-Stokes (RANS) equations coupled with the Spalart Allmaras turbulence model. Secondly, this work aims at estimating the influence of viscosity and turbulence on the T106A end-wall flow field. In order to do so, RANS results are compared with those obtained from an inviscid simulation with a prescribed inlet total pressure profile, which mimics a boundary layer. A comparison between C1 and C2 results highlights an influence of secondary flows on the flow field up to a significant distance from the end-wall. In particular, the C2 end-wall flow field appears to be characterized by greater over turning and under turning angles and higher total pressure losses. Furthermore, the C2 simulated flow field shows good agreement with experimental and numerical data available in literature. The C2 and inviscid Euler computed flow fields, although globally comparable, present evident differences. The cascade passage simulated with inviscid flow is mainly dominated by a single large and homogeneous vortex structure, less stretched in the spanwise direction and closer to the end-wall than vortical structures computed by compressible flow simulation. It is reasonable, then, asserting that for the chosen test case a great part of the secondary flows details is strongly dependent on viscous phenomena and turbulence.

• Journal of the korean society of oceanography
• /
• v.34 no.1
• /
• pp.1-21
• /
• 1999

A two-layer quasi-geostrophic numerical model is used to investigate the temporal variability of the East Korea Warm Current (EKWC), especially the separation from the Korean coast and the generation of warm eddies. An attention is given on the active role of the nonlinear boundary layer process. For this, an idealized flat bottom model of the East Sea is forced with the annual mean wind curl and with the inflow-outflow specified at the Korea (Tsushima) and Tsugaru Straits. Two types of separation mechanisms are identified. The first one is influenced by the westward movement of the recirculating leg of the EKWC (externally driven separation),the second one is solely driven by the boundary layer dynamics (internally driven separation). However, these two processes are not independent, and usually coexist. It is hypothesized that 'internally driven separation' arises as the result of relative vorticity production at the wall, its subsequent advection via the EKWC, and its accumulation up to a critical level characterized by the separation of the boundary flow from the coast. It is found that the sharp southeastern corner of the Korean peninsula provides a favorable condition for the accumulation of relative vorticity. The separation of the EKWC usually accompanies the generation of a warm eddy with a diameter of about 120 km. The warm eddy has a typical layer-averaged velocity of 0.3 m/s and its lifespan is up to a year. In general, the characteristics of the simulated warm eddy are compatible with observations. A conclusion is therefore drawn that the variability of the EKWC is at least partially self-excited, not being influenced by any sources of perturbation in the forcing field, and that the likely source of the variability is the barotropic instability although the extent of contribution from the baroclinic instability remains unknown. The effects of the seasonal wind curl and inflow-outflow strength are also investigated.

• The Journal of the Acoustical Society of Korea
• /
• v.39 no.4
• /
• pp.364-370
• /
• 2020

Rocket noise generated from the exhaust plume produces the enormous acoustic loading, which adversely affects the integrity of the electronic components and payload (satellite) at liftoff. The prediction of rocket noise consists of two steps: the supersonic jet exhaust is simulated by a method of the Computational Fluid Dynamics (CFD), and an acoustic transport method, such as the Helmholtz-Kirchhoff integral, is applied to predict the noise field. One of the difficulties in the CFD step is to remove the boundary reflection artifacts from the finite computation boundary. In general, artificial damping, known as a sponge layer, is added nearby the boundary to attenuate these reflected waves but this layer demands a large computational area and an optimization procedure of related parameters. In this paper, a cost-efficient way to separate the reflected waves based on the two microphone method is firstly introduced and applied to the computation result of a laboratory-scale supersonic jet noise without sponge layers.

• 전기의세계
• /
• v.29 no.4
• /
• pp.256-259
• /
• 1980

A singular pertubation theory is applied to obtain an approximate solution for suboptimal control of nuclear reactors with spatially distributed parameters. The inverse of the neutron velocity is regarded as a small perturbing parameter, and the model, adopted for simplicity, is a cylindrically symmetrical reactor whose dynamics are described by the one group diffusion equation with one delayed neutron group. The Helmholtz mode expansion is used for the application of the optimal theory for lumped parameter systems to the spatially distributed parameter systems. An asymptotic expansion of the feedback gain matrix is obtained with construction of the boundary layer correction up to the first order.