• Journal of The Korean Astronomical Society
• /
• v.37 no.1
• /
• pp.55-59
• /
• 2004

We have developed a two fluid solar wind model from the Sun to 1 AU. Its basic equations are mass, momentum and energy conservations. In these equations, we include a wave mechanism of heating the corona and accelerating the wind. The two fluid model takes into account the power spectrum of Alfvenic wave fluctuation. Model computations have been made to fit observational constraints such as electron($T_e$) and proton($T_p$) temperatures and solar wind speed(V) at 1 AU. As a result, we obtained physical quantities of solar wind as follows: $T_e$ is $7.4{\times}10^5$ K and density(n) is $1.7 {\times}10^7\;cm^{-3}$ in the corona. At 1 AU $T_e$ is $2.1 {\times} 10^5$ K and n is $0.3 cm^{-3}$, and V is $511 km\;s^{-1}$. Our model well explains the heating of protons in the corona and the acceleration of the solar wind.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.557-560
• /
• 2002

In this study, simulation of weak shock waves are peformed by a two-dimensional thermal fluid or compressible fluid model of the lattice Boltzmann method. The shock wave represents an abrupt change in fluids properties, in which finite variations in pressure, internal energies, and density occur over the shock thickness. The characteristics of the proposed model with a simple distribution function is verified by calculation of the sound speeds, and the shock tube problem. The reflection of a weak shock wave by wedge propagating in a channel is performed. The results agree well with those by finite difference method or by experiment. In the simulation of unsteady shock wave diffraction around a sharp corner, we show a flow field of vortical structure near the comer.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.2005-2011
• /
• 2003

Introducing the interfacial pressure jump terms based on the surface tension into the momentum equations of two-phase two-fluid model, the system of governing equations is turned mathematically into the hyperbolic system. The eigenvalues of the equation system become always real representing the void wave and the pressure wave propagation speeds as shown in the previous manuscript. To solve the interfacial pressure jump terms with void fraction gradients implicitly, the conventional semi-implicit method should be modified as an intermediate iteration method for void fraction at fractional time step. This advanced semi-implicit method (ASIM) then becomes stable without conventional additive terms. As a consequence, including the interfacial pressure jump terms with the advanced semi-implicit method, the numerical solutions of typical two-phase problems can be more stable and sound than those calculated exclusively by using any other terms like virtual mass, or artificial viscosity.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.22 no.7
• /
• pp.401-409
• /
• 2018

Performance-based seismic evaluation is usually done by considering simplified models for the liquid storage tanks therefore, it is important to validate those simplified models before conducting such evaluation. The purpose of this study is to compare the seismic response results of the FSI (fluid-structure interaction) model and the simplified models for the cylindrical liquid storage tanks and to verify the applicability of the simplified models for estimating failure probability. Seismic analyses were carried out for two types of storage tanks with different aspect ratios (H/D) of 0.45 and 0.86. FSI model represents detailed 3D fluid-structure interaction model and simplified models are modeled as cantilever mass-spring model, frame type mass-spring model and shell type mass-spring model, considering impulsive and convective components. Seismic analyses were performed with modal analysis followed by time history analysis. Analysis results from all the models were verified by comparing with the results calculated by the code and literature. The results from simplified models show good agreement with the ones from detailed FSI model and calculated results from code and literature, confirming that all three types of simplified models are very valid for conducting failure probability analysis of the cylindrical liquid storage tanks.

• Wind and Structures
• /
• v.22 no.3
• /
• pp.349-368
• /
• 2016

A preconditioning technique is presented for a simultaneous solution to wind-membrane interaction. In the simultaneous equations, a linear elastic model was employed to deal with the fluid-structure data transfer at the interface. A Lagrange multiplier was introduced to impose the specified boundary conditions at the interface and strongly coupled simultaneous equations are derived after space and time discretization. An initial linear elastic model preconditioner and modified one were derived by treating the linearized elastic model equation as a saddle point problem, respectively. Accordingly, initial and modified fluid-structure interaction (FSI) preconditioner for the simultaneous equations were derived based on the initial and modified linear elastic model preconditioners, respectively. Wind-membrane interaction analysis by the proposed preconditioners, for two and three dimensional membranous structures respectively, was performed. Comparison was made between the performance of initial and modified preconditioners by comparing parameters such as iteration numbers, relative residuals and convergence in FSI computation. The results show that the proposed preconditioning technique greatly improves calculation accuracy and efficiency. The priority of the modified FSI preconditioner is verified. The proposed preconditioning technique provides an efficient solution procedure and paves the way for practical application of simultaneous solution for wind-structure interaction computation.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2011.11a
• /
• pp.50-54
• /
• 2011

This paper presents a two-phase model for hybrid rocket internal ballistics design using $N_2O$ as oxidizer The two-phase model results are compared with data obtained from static firing test. Two-phase model is suitable for blow-down type with saturated compressible fluid as $N_2O$, presented the result by Part 1. HDPE as Fuel, and $N_2O$ as oxidizer were used during the static firing test. The combustor were designed for an average thrust of 30 kgf where oxidizer tank pressure in set to 50 bar. The numerical results of internal ballistic showed good agreements with static firing test results where thrust, oxidizer tank pressure and chamber pressure are compared.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.13 no.1
• /
• pp.585-598
• /
• 2021

The course-keeping stability of a high speed planing boat should be considered at the design stage for its safe operations. The shape of waterjet intake plane is one of important design parameters of a waterjet propelled planing boat. That has significant influences on the stern flow patterns and pressure distributions. In this study, the effects of the waterjet intake shapes of planing boats on the course-keeping stabilities are investigated. Two kinds of designed planing boats have the same dimensions, but there are differences in waterjet intake plane shapes. Captive and free-running model tests, Computational Fluid Dynamics (CFD) analyses are carried out in order to estimate their hydrodynamic performances including course-keeping stabilities. The results show that the flat and wide waterjet intake plane of the initially designed boat makes the course-keeping stability worse. The waterjet intake shape is redesigned to improve the course-keeping stability. The improved performances are confirmed by free-running model tests and full-scale trials.

• The KSFM Journal of Fluid Machinery
• /
• v.17 no.2
• /
• pp.12-23
• /
• 2014

A parametric study on anti-vortex holes for turbine blade cooling was investigated numerically. Three-dimensional Reynolds-averaged Navier-Stokes equations and shear stress transport turbulence model were used for analysis of anti-vortex film cooling. Validation of numerical results was carried out comparing with experimental data. The cooling performance of anti-vortex holes was assessed by two geometric variables, the ratio of diameters of holes and the lateral distances between the primary hole and anti-vortex hole at blowing ratios of 0.5 and 1.0. The results showed that the spatially-averaged film-cooling effectiveness increases as the ratio of the diameters increases and the distance between the primary hole and anti-vortex hole decreases.

• Plant Journal
• /
• v.14 no.1
• /
• pp.42-46
• /
• 2018

The flow and heat transfer characteristics inside a solar chimney power plant system are analyzed in this article. 3-D model with the $k-{\varepsilon}$ turbulence closure was developed. In this model, to solve the radiative transfer equation the discrete ordinates radiation model was implemented, using a two-band radiation model. To simulate radiation effects from the sun's rays, the solar ray tracing algorithm was coupled to the calculation via a source term in the energy equation. Simulations were carried out for a system with the geometry parameters of the Manzanares power plant. Based on the numerical results, the velocity and temperature distributions were illustrated and the results were validated by comparing with experimental data of the Manzanares prototype power plant. Moreover, temperature profile of the ground surface of the system was illustrated.

• Journal of Ocean Engineering and Technology
• /
• v.29 no.2
• /
• pp.163-168
• /
• 2015

This paper presents the results of experimental work on the wedge slamming impact problem. An experiment was done with a wedge model. The deadrise angle of the wedge was $4^{\circ}$. The model was made in two parts: the outside part was made of a 5-mm-thick steel plate that could be assumed to be a rigid body, and the inside part was made of a thin SUS plate that could be assumed to be an elastic body. Thin SUS plate thicknesses of 2 mm and 3 mm were used to determine the effect of plate rigidity. The drop height was varied from 0.25 m to 1 m to determine the effect of a large deformation.