• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2010.11a
• /
• pp.729-732
• /
• 2010

The present study has been motivated by the development of a reliable numerical methodology for simulation of kerosene/LOx coaxial swirl injectors. To deal with thermodynamic non-ideality and anomalies of transport properties pronounced at supercritical pressures, a set of subroutine libraries has been constructed based on the cubic equations of state, and applied to an existing flamelet analysis code. For computational efficiency, two-dimensional axisymmetric RANS formulation with swirl was adopted and validated successfully against an isothermal coaxial swirling jet. For the actual problem with high pressure combustion, however, numerical results show that the RANS models yield excessive production of turbulence probably due to high density gradient magnitude in the vicinity of mixing layer of swirling film flow, and imply strongly further improvement of the turbulence models.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.699-702
• /
• 2004

In this study, the affect of mounting axisymmetrical supersonic inlet to airfoil, which has 65 degree swept angle was numerically investigated. The parameter for this calculation are tree stream Mach number M=2.0 and 2.5, the distance between inlet spike and airfoil lower surface $L_{sw}$/$R_{cowl}$ = 1.21-1.54 and angle of attack to the airfoil 0-4. The mass capture ratio improved 3points in M=2.0 condition and 1points in M=2.5 while the mass capture ratio without airfoil surface was 57% and 71 % for each case. These are the result from increase of density and change of velocity deflection by the shock wave structure formed between inlet and airfoil surface. On the other hand, the distortion of Mach number at cowl lip plane increased by 13% in M=2.0, 3% in M=2.5 condition. The effects of the angle attack on the mass capture ratio is greater than that of the shock wave interaction between inlet and cowl, but the effects to the distortion is smaller in the range of this calculation condition. In the condition of M=2.0 with 4 degrees of angle of attack, inlet distortion of Mach number is mainly caused by the affects of the shock wave interaction between inlet and airfoil surface, while the largest angle of the velocity vector in the radial direction at cowl lip plane is caused by the affect of angle of attack. This large velocity vector made the flow inside the cowl subsonic and caused spillage, which interfere with the boundary layer of airfoil surface.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.11 no.11
• /
• pp.4664-4670
• /
• 2010

Comparison of the classical mobile-immobile zone (MIM) model to the derived model led to several conclusions. If the MIM model is to be applied, the initial concentration in the immobile zone has to be down-scaled by a correction factor that is a function of pore geometry. The MIM model was valid only after sufficiently long time has passed, i.e., only after the diffusion front reaches the deepest pore wall in the immobile zone. The MIM mass-transfer coefficient $\alpha$, was inversely proportional to the square of the pore depth. Also it did not depend on the mobile-zone flow velocity, contrary to the number of laboratory and field observations. The classical MIM model displayed a rapid exponential decay of immobile-zone concentration. Meanwhile at large times, the newly derived model displayed similar exponential decay. This was contrary to the mounting evidence of power-law BTC tails observed in laboratory and field settings.

• Coupled systems mechanics
• /
• v.12 no.1
• /
• pp.41-68
• /
• 2023

The present work is concerned with the study of the influence of inhomogeneous initial stresses in a hollow cylinder containing a compressible inviscid fluid on the propagation of axisymmetric longitudinal waves propagating in this cylinder. The study is carried out using the so-called three-dimensional linearized theory of elastic waves in bodies with initial stresses to describe the motion of the cylinder and using the linearized Euler equations to describe the flow of the compressible inviscid fluid. It is assumed that the inhomogeneous initial stresses in the cylinder are caused by the internal pressure of the fluid. To solve the corresponding eigenvalue problem, the discrete-analytic solution method is applied and the corresponding dispersion equation is obtained, which is solved numerically, after which the corresponding dispersion curves are constructed and analyzed. To obtain these dispersion curves, parameters characterizing the magnitude of the internal pressure, the ratio of the sound velocities in the cylinder material and in the fluid, and the ratio of the material densities of the fluid and the cylinder are introduced. Based on these parameters, the influence of the inhomogeneous initial stresses in the cylinder on the dispersion of the above-mentioned waves in the considered hydro-elastic system is investigated. Moreover, based on these results, appropriate conclusions about this influence are drawn. In particular, it is found that the character of the influence depends on the wavelength. Accordingly, the inhomogeneous initial stresses before (after) a certain value of the wavelength lead to a decrease (increase) of the wave propagation velocity in the zeroth and first modes.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.18 no.4
• /
• pp.274-281
• /
• 2015

Thermal and nuclear power plants on shore commonly use the sea water for cooling facility. Discharged cooling water has the high kinematic energy potential due to amount of water flux. Numerical analysis was made to find the suitable combinations between the arrangement of tidal turbines and the overall dimensions of the discharged channel. Several parameters such as the turbine diameter to inlet size, and the axial distance to turbine size were investigated. Power coefficients for various test conditions were also compared to see the effect of inlet configurations such as single inlet and dual inlet. For the single inlet, the mean power coefficient appeared to be gradually decreased with increasing distance, and the maximum power was obtained when the turbine diameter was same as the inlet diameter. For the dual inlet, the tendency was similar so that the better result when the turbine diameter was same as the inlet diameter. It is expected that the present methodology can be extensively utilized to harness the high kinetic energy flow of the discharge channel of power plant.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.6
• /
• pp.436-449
• /
• 2013

In this paper, a numerical analysis is carried out to study the characteristics of supercavitating flows and the drag of relatively simple two-dimensional and axisymmetric bodies which can be used for supercavity generation device, cavitator, of a high-speed underwater vehicle. In order to investigate the suitability of numerical models, cavity flows around the hemispherical head form and two-dimensional wedge are calculated with combinations of three turbulence models(standard $k-{\epsilon}$, realizable $k-{\epsilon}$, Reynolds stress) and two cavitation models(Schnerr-Sauer, Zwart-Gerber-Belamri). From the results, it is confirmed that the calculated cavity flow is more affected by the turbulence model than the cavitation model. For the calculation of steady state cavity flows, the convergence in case of the realizable $k-{\epsilon}$ model is better than the other turbulence models. The numerical result of the Schnerr-Sauer cavitation model is changed less by turbulence model and more robust than the Zwart-Gerber-Belamri model. Thus the realizable $k-{\epsilon}$ turbulence model and the Schnerr-Sauer cavitation model are applied to calculate supercavitating flows around disks, two dimensional $10^{\circ}$ and $30^{\circ}$ wedges. In case of the disk, the cavitation number dependences of the cavity size and the drag coefficient predicted are similar to either experimental data or Reichardt's semi-empirical equations, but the drag coefficient is overestimated about 3% higher than the Reichardt's equation. In case of the wedges, the cavitation number dependences of the cavity size are similar to experimental data and Newman's linear theory, and the agreement of the cavity length predicted and Newman's linear theory becomes better as decreasing cavitation number. However, the drag coefficients of wedges agree more with experimental data than those of Newman's analytic solution. The cavitation number dependences of the drag coefficients of both the disk and the wedge appear linear and simple formula for estimating the drag of supercavitating disks and wedges are suggested. Consequently, the CFD scheme of this study can be applied for numerical analysis of supercavitating flows of the cavitator and the cavitator design.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.14 no.11
• /
• pp.5906-5914
• /
• 2013

Recently in Korea, the policy is being proceeded to build a intergenerational housing on artificial ground of railroad site for utilizing rental house. Due to narrow space of rail road site, suitable method have to be developed such as micropiles which is known as a method of a fast construction. However, If micropile is used as foundations for the super structure, construction cost is increases compared with other pile. Consequently, new concept micropile proposed to improve both bearing capacity and cost efficiency of general micropile. New concept micropile consists of waveform cement grout surrounding tread bar that formed by grouting the soil layer with jet grouting method as control the grout pressure and flow. The micropile with waveform is expected to decrease the construction cost by cut down pile length of general micropile. This paper examined the behavior of the new concept micropile with waveform subjected to axial load using two-dimensional axisymmetric numerical analyses method. According to the numerical result, there will cost effectiveness as the pile displacement decreased despite the length of waveform micropile is down about 5% from a general micropile under the same loading condition. Also, the effect of skin friction force which mobilized from the waveform of micropile appeared at relatively soft ground.

• Journal of the Society of Naval Architects of Korea
• /
• v.55 no.3
• /
• pp.252-264
• /
• 2018

The calculation of steady-state cavitating flows around Supercavitating Underwater Bodies (SUB's), which consist of a circular disk head (cavitator), a conical fore-body, a cylindrical middle-body and either a boat-tail or a flare-tail, are carried out. To calculate the axisymmetric cavitating flow, used is a commercial computational fluid dynamics code based on the finite volume method, Fluent. From the analysis of numerical results, the cavity and drag, affected by the fore-body and tail of the SUB's, are investigated. Firstly, the effect of the fore-body shape is investigated with the same disk cavitator and a cylindrical rear-body of fixed diameter. Then with the same cavitator and a fixed fore-body, the effect of the rear-body shape is investigated. Before the cavity generated by the cavitator covers the slant of fore-bodies sufficiently, the larger the cone angle of the fore-body(i.e., the shorter the slant length), the larger the drag and the slower the development of cavity. After the cavity covers the fore-body completely so that the pressure drag component of the body is vanished, the characteristics of drag-velocity curves are identical. Also, as the tail angle is bigger, the cavity generated by the cavitator is suppressed further and the drag becomes larger. The peak of the drag appears for the flare-tail, i.e., when the tail angle is positive(+). On the contrary, the trough of the drag appears for the boat-tail, i.e., when the tail angle is negative(-). When the tail angle is 5 degrees, the peak of the drag appears at the body speed of 80m/s and the value of the drag is 43% larger than that at the design speed of 100m/s. When the tail angle is -5 degrees, the trough of the total drag appears at 75m/s and that drag is 30% smaller than that of the cavitator, which means the rest of the body has a negative drag.