• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.155-160
• /
• 2004

This paper describes a supersonic separation of air-launching rocket from the mother plane. Three dimensional Euler equations were numerically solved to analyze steady/unsteady state fluid flows. To solve the Euler equations, named CFD-FASTRAN that is commercial computation code. The results of simulation clearly demonstrate effect of shock-expansion wave interaction between the rocket and the mother plane. Moreover, important influential factors at separating stage of the rocket were extracted with a comprehensive analysis. Finally, from the consideration of supersonic-separation, a guideline to safety-separation is given to the design of supersonic air-launching rocket.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.110-119
• /
• 2006

An accurate and cost efficient method PSE is used for the stability analysis of 2D or 3D compressible boundary layers. A highly accurate finite difference PSE code has been developed at a general curvilinear coordinate system using an implicit marching procedure to deal with a broad range of transition predictions problems. Evolution of disturbances in compressible flat plate boundary layers are studied for free-stream Mach numbers ranging from 0 to 1.5. The effect of mean-flow nonparallelism is found to be weak on two dimensional waves and strong on three dimensional waves. The maximum amplification rate increases monotonically with Mach number. The present PSE solutions are compared with previous numerical investigations and experimental results and are found to be in good agreement.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.671-674
• /
• 2002

Pressure compensating temperature control valve(TCV) is one of the important control devices, which is used to maintain the constant temperature of working fluid in power and chemical plants. The ratio of cylinder hole diameters of inlet and outlet is the main design parameters of TCV. So this needs to be investigated to improve the function of control of temperature and void fraction. In this study, numerical analysis is carried out with various ratios of cylinder hole diameters of the inlet and outlet in the TCV. Especial1y, the distribution of the static pressure Is investigated to calculate the new coefficient($C_{\upsilon}$) and resistance coefficient(K). The governing equations are derived from making using of three-dimensional Naver-Stokes equations with standard $k-{\varepsilon}$ turbulence model and SIMPLE algorithm. Using a commercial code, PHOENICS, pressure and flow fields in TCV are calculated with different inlet and outlet diameters of the cylinder hole for cold and hot water passages.

• Coupled systems mechanics
• /
• v.7 no.1
• /
• pp.5-25
• /
• 2018

A fully-coupled thermodynamic-based transient finite element formulation is proposed in this article for electric, magnetic, thermal and mechanic fields interactions limited to the linear case. The governing equations are obtained from conservation principles for both electric and magnetic flux, momentum and energy. A full-interaction among different fields is defined through Helmholtz free-energy potential, which provides that the constitutive equations for corresponding dual variables can be derived consistently. Although the behavior of the material is linear, the coupled interactions with the other fields are not considered limited to the linear case. The implementation is carried out in a research version of the research computer code FEAP by using 8-node isoparametric 3D solid elements. A range of numerical examples are run with the proposed element, from the relatively simple cases of piezoelectric, piezomagnetic, thermoelastic to more complicated combined coupled cases such as piezo-pyro-electric, or piezo-electro-magnetic. In this paper, some of those interactions are illustrated and discussed for a simple geometry.

• Journal of the Korean Society of Tobacco Science
• /
• v.20 no.1
• /
• pp.87-94
• /
• 1998

After we made the computer source code with mathematical model of Muramatsu et al. that was expressed by the set of simultaneous first-order ordinary differential equations in evaporation-pyrolysis zone of cigarette, we simulated the distribution profiles of temperature and density of flue-cured tobacco. Those equations were solved numerically with the Runge-Kutta-Gill algorithm assuming step size of 0.025mm by Muramatsu at at,, but in this study the advanced algorithm of Runge-Kutta 4th Order assuming step size of 0.0005mm. The initial conditions and physical parameters of Muramatsu et at. were used for solving them. The calculated values corresponded well with results of Muramatsu et al., especially the gradient of the temperature profile increased with smoldering speed and the thickness of the evaporation-pyrolysis zone decreased with increasing of smoldering speed. On the other hand, the temperature gradient decreased with increasing of the effective thermal-conductivity value and the thickness of the evaporation-pyrolysis zone increased with the effective thermal-conductivity value.

• Journal of Welding and Joining
• /
• v.18 no.4
• /
• pp.76-81
• /
• 2000

In this research, various heat source equations that have been proposed in previous study were calculated and compared with new model in various laser parameters. This is to treat the problem of predicting, by numerical analysis, the thermo-mechanical behaviors of laser spot welding for thin stainless steel plates. A finite element code, ABAQUS is used for the heat transfer analysis with a three-dimensional plane assumption. Experimental studies if the laser spot welding have also bee conducted to validate the numerical models presented. The results suggest that temperature profiles and weld dimensions are varied according to the heat source of the laser beam. For this reason, it is essential to incorporate an accurate description of the heat source.

• Journal of applied mathematics & informatics
• /
• v.7 no.3
• /
• pp.905-915
• /
• 2000

In this article we discuss some aspects of operational Tau Method on delay differential equations and then we apply this method on the differential delay equation defined by $\omega(u)\;=\frac{1}{u}\;for\;1\lequ\leq2$ and $(u\omega(u))'\;=\omega(u-1)\;foru\geq2$, which was introduced by Buchstab. As Khajah et al.[1] applied the Recursive Tau Method on this problem, they had to apply that Method under the Mathematica software to get reasonable accuracy. We present very good results obtained just by applying the Operational Tau Method using a Fortran code. The results show that we can obtain as much accuracy as is allowed by the Fortran compiler and the machine-limitations. The easy applications and reported results concerning the Operational Tau are again confirming the numerical capabilities of this Method to handle problems in different applications.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.138-144
• /
• 1998

The BGK equation, which is the kinetic model equation of Boltzmann equation, is solved using FDDO(finite difference with the discrete-ordinate method) to compute the rarefied flow of monatomic gas. Using reduced velocity distribution and discrete ordinate method, the scalar equation is transformed into a system of hyperbolic equations. High resolution ENO(Essentially Non-Oscillatory) scheme based on Harten-Yee's MFA(Modified Flux Approach) method with Strang-type explicit time integration is applied to solve the system equations. The calculated results are well compared with the experimental density field of NACA0012 airfoil, validating the developed computer code. Next. the computed results of circular cylinder flow for various Knudsen numbers are compared with the DSMC(Direct Simulation Monte Carlo) results by Vogenitz et al. The present scheme is found to be useful and efficient far the analysis of two-dimensional rarefied gas flows, especially in the transitional flow regime, when compared with the DSMC method.

• 한국신재생에너지학회:학술대회논문집
• /
• 2009.11a
• /
• pp.143-146
• /
• 2009

The metal-supported SOFC has beed developed as a new concept of SOFC which has higher mechanical strength. However, the mass transfer rate in this type of SOFC may be decreased due to the contact layer and the support layer and that can cause the low performance. Therefore, numerical analysis of the heat and mass transfer characteristics in a metal-supported solid oxide fuel cell(SOFC) is studied in this paper. Governing equations and electrochemical equations are calculated simultaneously. And the numerical results are compared with the experimental results for the code validation. The current density, temperature, and pressure drop are suggested as numerical results.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.2
• /
• pp.103-121
• /
• 2015

In this paper, we briefly review and describe a projection algorithm for numerically computing the two-dimensional time-dependent incompressible Navier-Stokes equation. The projection method, which was originally introduced by Alexandre Chorin [A.J. Chorin, Numerical solution of the Navier-Stokes equations, Math. Comput., 22 (1968), pp. 745-762], is an effective numerical method for solving time-dependent incompressible fluid flow problems. The key advantage of the projection method is that we do not compute the momentum and the continuity equations at the same time, which is computationally difficult and costly. In the projection method, we compute an intermediate velocity vector field that is then projected onto divergence-free fields to recover the divergence-free velocity. Numerical solutions for flows inside a driven cavity are presented. We also provide the source code for the programs so that interested readers can modify the programs and adapt them for their own purposes.