• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.1
• /
• pp.252-261
• /
• 1991

A particle trajectory model to simulate two-phase particle-laden crossjets into two-dimensional horizontal free stream has been developed to study the variations of the jet trajectories and velocity variations of the gaseous and the particulate phases. The following conclusions may be drawn from the predicted results, which are in agreement with experimental observations. The penetration of the two-phase jet in a crossflow is greater than that of the single-phase jet. The penetration of particles into the free stream increases with increasing particle size, solids-gas loading ratio and carrier gas to free stream velocity ratio at the jet exit. When the particle size is large, the solid particles separate from the carrier gas , while the particles are completely suspended in the carrier gas for the case of small size particles. As the particle to carrier gas velocity ratio at the jet exit is less than unity, the particles in the vicinity of the jet exit are accelerated by the carrier gas. As the injection angle is increased, the difference of the particle trajectory from that of the pure gas becomes larger. Therefore, it can be concluded that the velocities and trajectories of the particle-laden jets in a crossflow change depending on the solids-gas loading ratio, particle size, carrier gas to free stream velocity ratio and particle to gas velocity ratio at the jet exit.

• Journal of Mechanical Science and Technology
• /
• v.20 no.12
• /
• pp.2218-2229
• /
• 2006

A conservative finite-volume numerical method for unstructured grids with the cell-centered method has been developed for computing flow and heat transfer by combining the attractive features of the existing pressure-based procedures with the advances made in unstructured grid techniques. This method uses an integral form of governing equations for arbitrary convex polyhedra. Care is taken in the discretization and solution procedure to avoid formulations that are cell-shape-specific. A collocated variable arrangement formulation is developed, i.e. all dependent variables such as pressure and velocity are stored at cell centers. For both convective and diffusive fluxes the forms superior to both accuracy and stability are particularly adopted and formulated through a systematic study on the existing approximation ones. Gradients required for the evaluation of diffusion fluxes and for second-order-accurate convective operators are computed by using a linear reconstruction based on the divergence theorem. Momentum interpolation is used to prevent the pressure checkerboarding and a segregated solution strategy is adopted to minimize the storage requirements with the pressure-velocity coupling by the SIMPLE algorithm. An algebraic solver using iterative preconditioned conjugate gradient method is used for the solution of linearized equations. The flow analysis code (PowerCFD) developed by the present method is evaluated for its application to several 2-D structured-mesh benchmark problems using a variety of unstructured quadrilateral and triangular meshes. The present flow analysis code by using unstructured grids with the cell-centered method clearly demonstrate the same accuracy and robustness as that for a typical structured mesh.

• 한국해양학회지
• /
• v.20 no.2
• /
• pp.28-39
• /
• 1985

The driving mechanisms for low level jets(LLJ) and coastal surface maximum winds are studied with observed wind data from June, 1976 through August, 1980 at Port Aransas and Victoria, Texas, in connection with a baroclinic model. This model is developed considering the forcing functions such as the synoptic and meso-scale pressure gradient, the frictional force, and the atmospheric stability. The results show that a LLJis observed on over 95% of the occasions when a nighttime coastal wind maximum occurred. Baroclinicity generated by sloping terrain during the summertime causes the diurnal variation in the thermal field. This thermal wind component would then decrease the prevailing synoptic-scale southerly wind by day and allow it to increase at night. Nighttime atmospheric stability leads to frictional decoupling which enhances the nocturnal LLJ. At the coastal site neutral stability prevails, thus all owing downward transfer of momentum from the nocturnal LLJ and results in the nocturnal coastal surface wind maximum. The height of LLJis not uniquely related to the inversion layer, and the results of the computations using this model show a good agreement with the observations.

• Journal of Astronomy and Space Sciences
• /
• v.24 no.2
• /
• pp.145-154
• /
• 2007

Antenna pointing analysis for a geostationary satellite has been performed for using the NORAD Two-Line-Elements (TLE) converted from osculating Keplerian orbital elements. In order to check the possibility of the reception of the satellite signal, the antenna offset angles have been derived for the Communications, Ocean, and Meteorological Satellite (COMS) which carries out weekly East-West and North-South station-keeping maneuvers and twice a day thruster assisted momentum dumping. Throughout the analysis, it is shown that the use of converted NORAD TLE simplifies the antenna pointing related interfaces in satellite mission control system. For a highly eccentric transfer orbit cases, further analysis presents that the converted NORAD TLE from near apogee gives more favorable results.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.10 no.4
• /
• pp.151-158
• /
• 1990

Various factors may contribute on the mixing processes in the surf zone formed by irregular waves. The turbulence motion driven by wave breaking may be one of the major causes, the effect due to spatial variation on current velocity be a secondary one, and the additional process may result from the irregular superposition of radiation stresses or wave breaking dissipation incurred by random breaking waves in a broadened surf zone. In the present study a numerical model of spectral waves and induced currents was developed using a superposition technique with ${\kappa}-{\varepsilon}$ closure for mixing process and applied to a field situation of longshore current generated by spectral waves on a uniform beach. It was found from the application that the surf-zone mixing processes formed by irregular waves can be well described by using ${\kappa}-{\varepsilon}$ equations if the source of ${\kappa}$ is properly represented. The nonlinear energy transfer was also found to have some influence on the velocity profile of longshore current particularly in very shallow water region near coast.

• Journal of the Korean Society for Precision Engineering
• /
• v.19 no.8
• /
• pp.92-99
• /
• 2002

This work investigated the optimal condition for an uniform deposition growth rate in the vertical cylindric CVD chamber. Heat transfer, surface chemical reaction and mass diffusion in the flow field of CVD chamber h,id been computed using Fluent v5.3 code. A SIMPLE based finite Volume Method (FVM) was adopted to solve the fully elliptic equations for momentum, temperature and concentration of a chemical species. The numerical analysis results show good agreements with the measurements obtained by N. Yoshikawa. The results obtained by the numerical analysis showed that the film growth rate in the center of a susceptor is increasing, as the inner flow approaches to the forced convection. To the contrast, as it approaches to the natural convection, that in the outside of a susceptor is increasing. As the Reynolds number increases, the uniformity may not hold due to the larger temperature gradient at a susceptor surface. Therefore, when the temperature gradient on the surface of a susceptor is zero, the film growth rate becomes uniform on most surface.

• Advances in Energy Research
• /
• v.5 no.2
• /
• pp.91-105
• /
• 2017

Thermal hydraulic (TH) analysis of nuclear power reactors is utmost important. In this way, the numerical codes that preparing TH data in reactor core are essential. In this paper, a subchannel analysis of a Russian pressurized water reactor (WWER1000) core with enhanced numerical code is carried out. For this, in fluid domain, the mass, axial and lateral momentum and energy conservation equations for desired control volume are solved, numerically. In the solid domain, the cylindrical heat transfer equation for calculation of radial temperature profile in fuel, gap and clad with finite difference and finite element solvers are considered. The dependence of material properties to fuel burnup with Calza-Bini fuel-gap model is implemented. This model is coupled with Isotope Generation and Depletion Code (ORIGEN2.1). The possibility of central hole consideration in fuel pellet is another advantage of this work. In addition, subchannel to subchannel and subchannel to rod connection data in hexagonal fuel assembly geometry could be prepared, automatically. For a demonstration of code capability, the steady state TH analysis of a the WWER1000 core is compromised with Thermal-hydraulic analysis code (COBRA-EN). By thermal hydraulic parameters averaging Fuel Assembly-to-Fuel Assembly method, the one sixth (symmetry) of the Boushehr Nuclear Power Plant (BNPP) core with regular subchannels are modeled. Comparison between the results of the work and COBRA-EN demonstrates some advantages of the presented code. Using the code the thermal modeling of the fuel rods with considering the fission gas generation would be possible. In addition, this code is compatible with neutronic codes for coupling. This method is faster and more accurate for symmetrical simulation of the core with acceptable results.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.2B
• /
• pp.165-173
• /
• 2011

In this study, FESWMS FST2DH model was used to analyze the change of flow characteristics after making the riparian forest. The additional flow resistance is calculated based on the drag-force concept acting on each tree and the lateral momentum transfer between planted and non-planted zone could be satisfactorily reproduced by parabolic turbulence model in this depth-averaged 2-D numerical model. For model validation, the simulated velocities were compared with the measured data, showing good agreement in both tree density cases of experiments. The previous method using a proper Manning's n coefficient gives reasonable solutions only to evaluate the conveyance, but the calculated approach velocity at each tree was different from realistic value. The proposed procedure could be widely used to evaluate hydraulic effects of riparian trees in practical engineering.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.6 no.1
• /
• pp.25-33
• /
• 1986

Natural river channels usually have a deep section and one or two floodplains, which is called a compound channel. As the general method in the compound channel overestimates the discharge capacity, the momentum transfer due to interaction between the main channel flow and flow over its floodplain must be considered. Scale model experiments are performed for the rectangular main channel with an asymmetrical floodplain. Firstly, velocities are measured at various section grids. Secondary, boundary shear stresses are calculated from velocity distributions. Lastly, in order to determine the apparent shear force, the shear stress distributions are integrated along the wetted perimeter for the full cross-section and equated to the total weight force in the flow direction. The hydraulic characteristics in a compound channel are closely examined with the scales of length, velocity, boundary shear stress, and apparent shear force which are described with the various relationships.

• 한국가시화정보학회:학술대회논문집
• /
• 2007.11a
• /
• pp.134-137
• /
• 2007

The temporal evolution of wake behind a circular cylinder oscillating rotationally with a relatively high forcing frequency has been investigated experimentally using a dynamic PIV technique. Experiments were carried out with varying the frequency ratio $F_R\;(=f_f/f_n)$ in the range from 0.0 (stationary) to 1.6 at oscillation amplitude of ${\theta}_A=30^{\circ}$ and Reynolds number of $Re=4.14{\times}10^3$. Depending on the forcing condition ($F_R$), the flow was divided into three regimes; non-lock-on ($F_R=0.4$), transition ($F_R=0.8$, 1.6) and lock-on regimes ($F_R=1.0$) with markedly different flow structure in the near-wake region behind the cylinder. When the frequency ratio was less than 1.0 ($F_R{\le}1.0$), the rotational oscillatory motion of the cylinder decreased the length of the vortex formation region and enhanced the mutual interaction between large-scale vortices across the wake centerline. The entrainment of ambient fluid seemed to play an important role in controlling the near-wake flow and shear-layer instability. However, the flow characteristics changed markedly beyond the lock-on flow regime ($F_R=1.0$) due to high-frequency forcing. At $F_R=1.6$, the mutual interactions between the vortices shed from both sides of the cylinder were not so strong. Thereby, the flow entrainment and momentum transfer into the wake center region were reduced. In addition, the size of the large-scale vortices decreased since the lateral extent of the wake was suppressed.