• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.3
• /
• pp.820-833
• /
• 1995

An analysis is performed to examine the equilibrium state and the stability of modeled Reynolds stress equations for homogeneous turbulent shear flows. The system of the governing equations consists of four coupled ordinary differential equations. The equilibrium states are found by the steady state solution of the governing equations. In order to investigate the stability of the system about its state in equilibrium, and eigenvalue problem is constructed. As a result, constraints for the coeffieients in the model equations are obtained by the stability condition of the equilibrium state as well as by their physically realizable bounds. It is observed that the models with pressure-strain rate correlation that are linear in the anisotropy tensor are stable and produce reasonable equilibrium tensor do not behave properly. Stability considerations about three most commonly used models are given in detail in the final section.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.6
• /
• pp.792-801
• /
• 2000

The problem of phase change from liquid to solid in the inviscid plane-stagnation flow is theoretically investigated. The solution at the initial stage of freezing is obtained by expanding it in powers of time, and the final equilibrium state is determined from the steady-state governing equations. The transient solution is dependent on the three dimensionless parameters, but the equilibrium state is determined by one parameter of (temperature ratio/conductivity ratio). The effect of the fluid flow on the growth rate of the solid in the pure conduction problem can be clearly seen from the solution of the initial stage and the final equilibrium state. The characteristics of the transient heat transfer at the surface of the solid and the liquid side of the solid-liquid interface for all the dimensionless parameters are elucidated.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.12 no.1
• /
• pp.1-11
• /
• 2000

This study investigates the problem of phase change from liquid to solid in the inviscid stagnation flow. The solution of dimensionless governing equations is determined by the three dimensionless parameters of (temperature ratio/conductivity ratio), Stefan number, and diffusi-vity ratio. The solution at the initial stage of freezing is obtained by expanding it in powers of time, and the final equilibrium state is determined from the steady-state governing equations. The equilibrium state is dependent on (temperature ratio/conductivity ratio), but is independent of Stefan number and diffusivity ratio. The effect of fluid flow on the pure conduction problem can be clearly seen from the solution of the initial stage and the final equilibrium state, and the characteristics of the solidification process for all the dimensionless parameters are elucidated.

• The Transactions of the Korean Institute of Electrical Engineers
• /
• v.32 no.10
• /
• pp.335-340
• /
• 1983

The analysis of static power converter circuit using state-space method is presented. Semiconductors are modeled in two-state resistors depending on their ON or OFF states. Then the modes of circuit are determined according to the conducting states of semiconductors and different describing matrices are given automatically for each mode. Newton-Raphson algorithm is used as an iterative method for obtaining steady-state solution and an adjoint network is introduced for the efficient and accurate evaluation of the Jacobi matrix in the algorithm. Using the porogram exploited from the above algorithm, it is shown through examples that the results are in good agreement with the analytic solutions and computation time is considerably reduced for obtaining the steady-state solutions.

• Tunnel and Underground Space
• /
• v.13 no.5
• /
• pp.389-396
• /
• 2003

It is shown that the cubic law can be modified regarding the steady-state Navier-Stokes equations by using perturbation approximation method for a sinusoidal aperture variation. In order to adopt the perturbation theory, the sinusoidal function needs to be non-dimensionalized for the amplitude and wavelength. Then, the steady-state Navier-Stokes equations can be solved by expanding the non-dimensionalized stream function with respect to the small value of the parameter (the ratio of the mean aperture to the wavelength), together with the continuity equation. From the approximate solution of the Navier-Stokes equations, the basic cubic law is successfully modified for the steady-state condition and a sinusoidal aperture variation. A finite difference method is adopted to calculate the pressure within a fracture model, and the results of numerical experiments show the accuracy and applicability of the modified cubic law. As a result, it is noted that the modified cubic law, suggested in this study, will be used for the analysis of fluid flow through aperture geometry of sinusoidal distributions.

• International Journal of Aeronautical and Space Sciences
• /
• v.3 no.1
• /
• pp.95-104
• /
• 2002

Numerical solution of the full Navier-Stokes equation as well as the energy equation has been obtained for the unsteady natural convection in a rectangular enclosure. One side wall was maintained at very high temperature simulating fires. Especially the effect of surface radiation was taken into account. While the enclosed air was assumed to be transparent, the internal walls directly interacted one another through the surface radiation. Due to a significant temperature difference in the flow field, the equation of state was used instead of the Boussinesq approximation. It was found that the rapid heating of the adiabatic ceiling and floor by the incoming radiation from the hot wall made the evolution at thermo-fluid field highly unstable in the initial period. Therefore, the secondary cells brought about at the floor region greatly affected the heat transfer mechanism inside the enclosure. The heat transfer rate was augmented by the radiation, resulting in requiring less time for the flow to reach the steady state. At the steady state neglecting radiation two internal hydraulic jumps were clearly observed in upper/left as well as in lower/right comer. However, the hydraulic jump in the lower/right comer could not be observed for the case including radiation due to its high momentum flow over the bottom wall. Radiation resulted in a faster establishment of the steady state phenomena.

• Nuclear Engineering and Technology
• /
• v.45 no.2
• /
• pp.191-202
• /
• 2013

A numerical analysis of thermal stratification in the upper plenum of the MONJU fast breeder reactor was performed. Calculations were performed for a 1/6 simplified model of the MONJU reactor using the commercial code, CFX-13. To better resolve the geometrically complex upper core structure of the MONJU reactor, the porous media approach was adopted for the simulation. First, a steady state solution was obtained and the transient solutions were then obtained for the turbine trip test conducted in December 1995. The time dependent inlet conditions for the mass flow rate and temperature were provided by JAEA. Good agreement with the experimental data was observed for steady state solution. The numerical solution of the transient analysis shows the formation of thermal stratification within the upper plenum of the reactor vessel during the turbine trip test. The temporal variations of temperature were predicted accurately by the present method in the initial rapid coastdown period (~300 seconds). However, transient numerical solutions show a faster thermal mixing than that observed in the experiment after the initial coastdown period. A nearly homogenization of the temperature field in the upper plenum is predicted after about 900 seconds, which is a much shorter-term thermal stratification than the experimental data indicates. This discrepancy may be due to the shortcoming of the turbulence models available in the CFX-13 code for a natural convection flow with thermal stratification.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.16 no.3
• /
• pp.181-192
• /
• 2012

Here in this paper, the steady paint film flow on a vertical wall of a non-Newtonian pseudo plastic fluid for drainage problem has been investigated. The exact solution of the nonlinear problem is obtained for the velocity profile. Also the average velocity, volume flux, shear stress on the wall, force to hold the wall in position and normal stress difference have been derived. We retrieve Newtonian case, when material constant ${\mu}_1$ and relaxation time ${\lambda}_1$ equal zero. The results for co-rotational Maxwell fluid is also obtained by taking material constant ${\mu}_1$ = 0. The effect of the zero shear viscosity ${\eta}_0$, the material constant ${\mu}_1$, the relaxation time ${\lambda}_1$ and gravitational force on the velocity profile for drainage problem are discussed and plotted.

• Journal of the Chungcheong Mathematical Society
• /
• v.31 no.1
• /
• pp.75-87
• /
• 2018

In this paper, we consider ratio-dependent predator-prey models with disease in the prey under Neumann boundary condition. We investigate sufficient conditions for the existence and non-existence of non-constant positive steady-state solutions by the effects of the induced diffusion rates.

• Journal of Power Electronics
• /
• v.18 no.5
• /
• pp.1398-1408
• /
• 2018

In order to further improve the steady-state control performance of model predictive torque control (MPTC), a double-vector-based model predictive torque control without a weighting factor is proposed in this paper. The extended voltage vectors synthesized by two basic voltage vectors are used to increase the number of feasible voltage vectors. Therefore, the control precision of the torque and the stator flux along with the steady-state performance can be improved. To avoid testing all of the feasible voltage vectors, the solution of deadbeat torque control is calculated to predict the reference voltage vector. Thus, the candidate voltage vectors, which need to be evaluated by a cost function, can be reduced based on the sector position of the predicted reference voltage vector. Furthermore, a cost function, which only includes a reference voltage tracking error, is designed to eliminate the weighting factor. Moreover, two voltage vectors are applied during one control period, and their durations are calculated based on the principle of reference voltage tracking error minimization. Finally, the proposed method is tested by simulations and experiments.