• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1572-1577
• /
• 2003

High-power pulsed laser ablation under atmospheric pressure is studied utilizing numerical and experimental methods with emphasis on recondensation ratio, and the dynamics of the laser induced vapor flow. In the numerical calculation, the temperature pressure, density and vaporization flux on a solid substrate are first obtained by a heat-transfer computation code based on the enthalpy method, and then the plume dynamics is calculated by using a commercial CFD package. To confirm the computation results, the probe beam deflection technique was utilized for measuring the propagation of a laser induced shock wave. Discontinuities of properties and velocity over the Knudsen layer were investigated. Related with the analysis of the jump condition, the effect of the recondesation ratio on the plume dynamics was examined by comparing the pressure, density, and mass fraction of ablated aluminum vapor. To consider the effect of mass transfer between the ablation plume and air, unlike the most previous investigations, the equation of species conservation is simultaneously solved with the Euler equations. Therefore the numerical model computes not only the propagation of the shock front but also the distribution of the aluminum vapor. To our knowledge, this is the first work that employed a commercial CFD code in the calculation of pulsed ablation phenomena.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.362-367
• /
• 2011

The aerodynamic performance of the pantograph of the next generation high sped train is analyzed. The calculation of the flow around pantograph is carried cut by FLUENT; by the steady state flow calculation with ${\kappa}-{\omega}$ SST turbulence model, the lift force of the pantograph is computed. For the verification of the numerical schemes am grid systems, flow calculations are performed with the pantograph shape which was used at the experiments performed at Railway Technical Research Institute (RTRI) in Japan. Then, the difference of lift force between numerical am experimental results is about 10%. Therefore, selected numerical schemes and the current grid system is adequate for the analysis am prediction of the aerodynamic performance of panthograph system. Based on these numerical schemes am grid system, the flow around pantograph of the next generation high sped train is calculated and the lift force of the pantograph is predicted; the lift force of the pantograph is about 146N.

• Steel and Composite Structures
• /
• v.38 no.1
• /
• pp.79-86
• /
• 2021

Due to the complex geological conditions, a large number of high quality coal seams was buried in the western of China which cannot be mining in open-pit methods. The dynamic properties of that coal cannot be studied easily in real site for the complex working condition. The compound coal blocks made on the basis of the real situation were studied in the laboratory. The physical and mechanical properties of the compound coal blocks and the raw coal were contrasted by using the UCS tests. The results show that the compound coal blocks made by mixing coal powder, cement and water in proportion of 2.5:2:1 are the closest to that of standard raw coal. Then the propagation of strain waves and crushing effects on the coal were studied in the compound coal blocks by using the super dynamic strain test system and the numerical calculated method of ANSYS/LS-DYNA. The results show that the diameter of the crushing zone in the compound coal blocks was similar to that in the numerical results. The fractures distribution in laboratory tests also has a similar trend to the calculation results. The measured strain waves at the distance of 50 cm, 100 cm, and 150 cm from the center of the charge are mainly concerned at -1.0×104 με and have a similar trend as that in the numerical simulation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.4
• /
• pp.433-440
• /
• 2004

Since the numerical analysis was adopted in the mold design, lots of computational methods have been proposed for the simulations of casting processes for the various shaped molds. Today, it is possible to simulate the filling and solidification processes of most casts using the VOF technique. Though the three-dimensional numerical model based on the Cartesian coordinate system can be applied to any shape of cast, it becomes very inefficient when the three-dimensional model is applied to the cast of axi-symmetrical shape since the control volume includes at least 11 of the physical model. In addition, the more meshes should be distributed along the circumferential boundaries of curved shape in the Cartesian coordinate system fur the better results, while such curved circumferential boundary does not need to be considered in the two-dimensional cylindrical coordinate system. This motivates the present study i.e. developing a two-dimensional numerical model for the axi-symmetrically shaped casts. The SIMPLER algorithm, the VOF method, and the equivalent specific heat method have been adopted in the combined algorithm for the flow calculation, the free surface tracking, and the phase change heat transfer, respectively. The numerical model has been applied to the casting process of a pulley, and it was proven that the mesh and time effective calculation was accomplished comparing to the calculation using three-dimensional model.

• Water for future
• /
• v.28 no.4
• /
• pp.125-136
• /
• 1995

ELM, a characteristic line based method, was applied to advection-dispersion equation, and the results obtained were compared with those of Eulerian schemes(Stone-Brian and QUICKEST). The calculation methods consisted of Lagrangian interpolation scheme and cubic spline interpolation scheme for the advection calculation, and the Crank-Nicholson scheme for the dispersion calculation. The results of numerical methods were as follows: (1) for Gaussian hill: ELM, using Lagrangian interpolation scheme, gave the most accurate computational result, ELM, using cubic spline interpolation scheme, and QUICKEST scheme gave numerical damping for Peclet number 50. Stone-Brian scheme gave phase shift introduced in the numerical solution for Peclet number 10 and 50. (2) for advanced front: All schemes gave accurate computational results for Peclet number 1 and 4. ELM, Lagrangian interpolation scheme, and Stone,Brian scheme gave dissipation error and ELM, using cubic spline interpolation scheme, and QUICKEST scheme gave numerical oscillation for Peclet number 50.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.706-711
• /
• 2003

A numerical method is presented for computing unsteady incompressible two-phase flows with immersed solids. The method is based on a level set technique for capturing the phase interface, which is modified to satisfy a contact angle condition at the solid-fluid interface as well as to achieve mass conservation during the whole calculation procedure. The modified level set method is applied for numerical simulation of bubble deformation in a micro channel with a cylindrical solid block and liquid jet from a micro nozzle.

• Proceedings of the Safety Management and Science Conference
• /
• 2009.11a
• /
• pp.565-572
• /
• 2009

This paper presents the optimization steps with weight and importance of estimated characteristic values in the multiresponse surface analysis(MRA). The research introduces the shape parameter of individual desirability function for relaxation and tighening of specification bounds. The study also proposes the combinded desirability function using arithmetic, geometric and harmonic means considering the measurement unit and numerical pattern.

• Journal of Korea Water Resources Association
• /
• v.34 no.1
• /
• pp.57-65
• /
• 2001

A numerical model for analyzing transcritical flow in open channel is tested to various cases of channel shape. As the numerical models developed for transcritical flow until now mainly focused on the application to only prismatic or hypothetical channels, there are some restrictions to apply the nonprismatic channels. In this study, to verify the accuracy and stability of second-order implicit ENO scheme, the numerical model was applied to the channels which haute the varying channel bed and width. Also the numerical model was applied to unsteady flow as well as steady flow. The study shows that the numerical model provides good accuracy in the calculation of stage and velocity with no numerical oscillation, particularly in the calculation of hydraulic jump and discontinous flow Then the implicit ENO scheme demonstrated good accuracy as a high-resolution scheme and stability as an implicit scheme.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.35 no.2
• /
• pp.41-48
• /
• 2023

In this study, simplified linear numerical method that can simulate wave generation and transformation by a moving bottom is introduced. Numerical analysis is conducted in wave number domain after continuity equation, linear dynamic and kinematic free surface boundary conditions and linear kinematic bottom boundary condition are Fourier transformed, and the results are expressed in space domain by an inverse Fourier transform. In the wavenumber domain, the dynamic free water surface boundary condition and the kinematic free water surface boundary condition are numerically calculated, and the velocity potential in the mean water level (z = 0) satisfies the continuity equation and the kinematic bottom boundary condition. Wave generation and transformation are investigated when the triangular and rectangular shape of bottoms move periodically. The results of the simplified numerical method are compared with the results of previous analytical solutions and agree well with them. Stability of numerical results according to the calculation time interval (Δt) and the calculation wave number interval (Δk) was also investigated. It was found that the numerical results were appropriate when Δt ≤ T(period)/1000 and Δk ≤ π/100.

• Journal of the Korean Society for Railway
• /
• v.13 no.1
• /
• pp.44-50
• /
• 2010

This study is devoted to understand the basic characteristics of the flowfield around a train in evacuated tube and to suggest an efficient numerical method to calculate the flowfield. To get steady-state solution in minimum calculation domain, various boundary condition have been tried for steady calculation and have been compared to the solution of unsteady calculation. At the train velocity of 300km/h, the aerodynamic drag results of both calculation method agreed very well. The drag ratio between on the open filed and in the tube from the calculation result by the suggested numerical method lied in the same fitting curve with that from the filed test of high-speed trains running in the line.