• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.97-103
• /
• 2008

Unsteady three-dimensional flowfields generated by transverse fuel injection into a supersonic mainstream are simulated with a DES turbulence model. Comparisons are made with experimental results in term of the temporal eddy position and eddy formation frequency. The vorticity field around the jet exit is also analyzed to understand the formation mechanism of the jet vortical structures. Results indicate that the DES model correctly predicts the convection characteristics of the large scale eddies. However, it is also observed that the numerical results slightly overpredict the eddy formation frequency. The jet vortical structures are developed from the competing vortices in the recirculation region of upstream boundary.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.97-103
• /
• 2008

Unsteady three-dimensional flowfields generated by transverse fuel injection into a supersonic mainstream are simulated with a DES turbulence model. Comparisons are made with experimental results in term of the temporal eddy position and eddy formation frequency. The vorticity field around the jet exit is also analyzed to understand the formation mechanism of the jet vortical structures. Results indicate that the DES model correctly predicts the convection characteristics of the large scale eddies. However, it is also observed that the numerical results slightly overpredict the eddy formation frequency. The jet vortical structures are developed from the competing vortices in the recirculation region of upstream boundary.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.296-303
• /
• 2000

Study on natural convection in a rectangular enclosure with a heating point has been studied by numerical and experimental methods. The governing equations were solved by a finite volume method, and for pressure term was used a SIMPLE algorithm. The parameters considered for numerical study are positions and surface temperature of a heating point i.e. Y/H = 0.25. 0.5. 0.75 and $11^{\circ}C{\leq}{\Delta}T{\leq}59^{\circ}C$. The results of isotherms and velocity vectors have been represented for various parameters. Based on the numerical data. the mean Nusselt number in the space can be expressed as a function of Grashof number

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.704-708
• /
• 2000

A new efficient numerical method for computing unsteady, incompressible flows, DRANS (Decoupled Reynolds-Averaged Navier-Stokes), is presented. To eliminate the restriction of CFL condition, a fully-implicit time advancement in which the Crank-Nicolson method is used fer both the diffusion and convection terms. is adopted. Based on decomposition method, the velocity-turbulent quantity decoupling is achieved. The additional decoupling of the intermediate velocity components in the convection term is made for the fully-implicit time advancement scheme. Since the iterative procedures for the momentum, ${\kappa}\;and\;{\varepsilon}$ equations are not required, the components decouplings bring fourth the reduction of computational cost. The second-order accuracy in time of the present numerical algorithm is ascertained by computing decaying vortices. The present decoupling method is applied to turbulent boundary layer with local forcing.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.19 no.4
• /
• pp.19-29
• /
• 1982

Galerkin's finite element method is applied to a two-dimensional heat convection-diffusion problem arising in the hydrodynamic lubrication of thrust bearings used in naval vessels. A parabolized thermal energy equation for the lubricant, and thermal diffusion equations for both bearing pad and the collar are treated together, with proper juncture conditions on the interface boundaries. it has been known that a numerical instability arises when the classical Galerkin's method, which is equivalent to a centered difference approximation, is applied to a parabolic-type partial differential equation. Probably the simplest remedy for this instability is to use a one-sided finite difference formula for the first derivative term in the finite difference method. However, in the present coupled heat convection-diffusion problem in which the governing equation is parabolized in a subdomain(Lubricant), uniformly stable numerical solutions for a wide range of the Peclet number are obtained in the numerical test based on Galerkin's classical finite element method. In the present numerical convergence errors in several error norms are presented in the first model problem. Additional numerical results for a more realistic bearing lubrication problem are presented for a second numerical model.

• Journal of Advanced Research in Ocean Engineering
• /
• v.2 no.2
• /
• pp.48-60
• /
• 2016

Accurate simulation of free-surface wave flows around a ship is very important for better hull-form design. In this paper, a computational fluid dynamics (CFD) code which is based on the open source libraries, OpenFOAM, was developed to predict the wave patterns around a ship. Additional anti-diffusion source term for minimizing a numerical diffusion, which was caused by convection differencing scheme, was considered in the volume-fraction transport equation. The influence of the anti-diffusion source term was tested by applying it to free-surface wave flow around the Wigley and KCS model ships. In results, the wave patterns and hull wave profiles of the Wigley and KCS model ships for various anti-diffusion coefficients showed quite close patterns. While, the band width of the water volume-fraction values between 0.1 to 0.9 at the Wigley and KCS model hull surfaces was narrowed by considering the anti-diffusion term. From the results, anti-diffusion source term decreased free-surface smearing.

• Journal of applied mathematics & informatics
• /
• v.28 no.1_2
• /
• pp.31-48
• /
• 2010

In this article, we consider singularly perturbed Boundary Value Problems(BVPs) for second order Ordinary Differential Equations (ODEs) with Neumann boundary condition and discontinuous source term. A parameter-uniform error bound for the solution is established using the Streamline-Diffusion Finite Element Method (SDFEM) on a piecewise uniform meshes. We prove that the method is almost second order of convergence in the maximum norm, independently of the perturbation parameter. Further we derive superconvergence results for scaled derivatives of solution of the same problem. Numerical results are provided to substantiate the theoretical results.

• Journal of The Korean Astronomical Society
• /
• v.19 no.2
• /
• pp.91-107
• /
• 1986

Making use of our extended version of $\ddot{O}pik's$ convection theory, we have calculated magnetic cycle periods of the sun and late type stars by using Parker's dynamo theory, where we have included the non-linear effect. We presented a relationship between the computed cycle period and spectral type to analyze observed magnetic activities of the late type stars and long-term luminosity variations. It is found that (1) the stellar magentic-cycle period increases towards the later spectral type, (2) the rapid rotation facilitates the activity-related luminosity variation of stars later than about K5, (3) differential rotation plays a critical role in determining the magnetic activity-cycle period, and (4) the non-local effect should be taken into account in order to understand the observed long-term luminosity variations.

• Journal of applied mathematics & informatics
• /
• v.23 no.1_2
• /
• pp.141-152
• /
• 2007

In this paper a numerical method is presented to solve singularly perturbed two points boundary value problems for second order ordinary differential equations consisting a discontinuous source term. First, in this method, an asymptotic expansion approximation of the solution of the boundary value problem is constructed using the basic ideas of a well known perturbation method WKB. Then some initial value problems and terminal value problems are constructed such that their solutions are the terms of this asymptotic expansion. These initial value problems are happened to be singularly perturbed problems and therefore fitted mesh method (Shishkin mesh) are used to solve these problems. Necessary error estimates are derived and examples provided to illustrate the method.

• Journal of Nuclear Fuel Cycle and Waste Technology(JNFCWT)
• /
• v.6 no.2
• /
• pp.119-128
• /
• 2008

Thermal assessment of a new CANDU spent fuel disposal system, which improves the retrievability of the spent fuel and enhances the densification factor compared with the Korean Reference disposal System, is carried out in this study. The canisters for CANDU spent fuels are stored for long term and cooled by natural convection in the proposed disposal system for the retrievability. The steady state thermal analyses for proposed CANDU disposal system are carried out with the ANSYS 10.0 CFX code. The thermal analyses are performed through two steps. At the first step, the sensitivity of the disposal tunnel spacing is analysed. The differences of maximum temperatures by several tunnel spacings are calculated at three points in the disposal tunnel. The result shows that the differences of the temperature at the three points are almost negligible because 99% of the decay heat is removed by natural convection. At the second procedure, 60m tunnel spacing with a ventilation system instead of natural convection is considered. The result is applied to the calculation of the canister surface temperature in disposal tunnel as boundary conditions. Consequently, the average and the maximum surface temperature of disposal canisters are $79.9^{\circ}C$ and $119^{\circ}C$, respectively. The inner maximum temperature of a basket in the disposal canister is calculated as $140.9^{\circ}C$. The maximum temperature of the basket meets the thermal requirement for the CANDU spent fuel cladding.