• Explosives and Blasting
• /
• v.23 no.3
• /
• pp.27-42
• /
• 2005

An explosion modeling technique was developed by using the spherical discrete element code, $PFC^{3D}$, which can be used to model the dynamic stress wave propagation phenomenon. The modeling technique is simply based on an idea that the explosion pressure should be applied to a $PFC^{3D}$ particle assembly not in the form of an external force (body force), but in the form of a contact force (surface force). The stress wave propagation modeling was conducted by simulating the experimental approach based on the Hopkinson's effect combined with the spatting phenomenon that had previously been developed to determine the dynamic tensile strength of Inada granite. As a result, the stress wave velocity obtained by the proposed modeling technique was 4167 m/s, which is merely $3\%$ lower than the actual wave velocity of 4300 m/s for an Inada granite.

• Bulletin of the Korean Space Science Society
• /
• 2009.10a
• /
• pp.26.2-26.2
• /
• 2009

The Solar and Space Weather Research Group (SOS) in Korea Astronomy and Space Science Institute (KASI) is constructing the Space Weather Prediction Center since 2007. As a part of the project, we are developing automatic real-time system of the global 3-D magnetohydrodynamics (MHD) simulation. The MHD simulation model of earth's magnetosphere is designed as modified leap-frog scheme by T. Ogino, and it was parallelized by using message passing interface (MPI). Our work focuses on the automatic processing about simulation of 3-D MHD model and visualization of the simulation results. We used PC cluster to compute, and virtual reality modeling language (VRML) file format to visualize the MHD simulation. The system can show the variation of earth's magnetosphere by the solar wind in quasi real time. For data assimilation we used four parameters from ACE data; density, pressure, velocity of solar wind, and z component of interplanetary magnetic field (IMF). In this paper, we performed some initial tests and made a animation. The automatic real-time system will be valuable tool to understand the configuration of the solar-terrestrial environment for space weather research.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.697-698
• /
• 2002

The dispersion of an aerosol bolus in acinus is analyzed numerically. Model geometry is a straight duct surrounded by an axisymmetric semicircular annulus which is expanding or contracting with breathing. Unsteady Wavier-Stokes equation is solved by CFX-F3D, an FVM commercial code and the trajectory of massless particle Is computed by Lagrangian method. For steady flow with no wall motion, mean velocity of aerosol bolus in alveolated duct is a little smaller than that in straight duct and dispersion in alveolated duct is comparable with the dispersion in straight tube. For expanding duct mean velocity of aerosol bolus approaches half of that in straight tube and effective diffusivity is smaller than that of straight tube. For contracting duct mean velocity of aerosol bolus becomes slightly larger than that in straight tube and effective diffusivity is comparable with the case of straight tube.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.18 no.11
• /
• pp.690-696
• /
• 2017

Various types of alternative energy sources to petroleum are being developed both domestically and internationally as clean energy that does not emit greenhouse gases. In particular, offshore wind power has been studied because the wind resources are relatively limitless and the wind power is relatively smaller than onshore. In this study, to analyze the scour phenomenon around offshore wind foundations, mono pile and tripod pile foundations were simulated using a FLOW-3D model. The scour phenomenon was evaluated for mono piles: one is a pile with a 5 m diameter and d=1.69 m and the other is a pile with a 5 m diameter. Numerical analysis showed that in the latter, the falling-flow increased and the maximum scour depth occurred more than 1.7 times. For a tripod pile foundation, the measured velocity and the maximum wave condition were applied to the upstream boundary condition, respectively, and the scour phenomenon was evaluated. When the maximum wave condition was applied, the maximum scour depth occurred more than about 1.3 times. When the LES model was applied, the scour depth reached equilibrium, whereas the numerical results of the RNG model show that the scour phenomenon occurred in the entire boundary area and the scour depth did not reach equilibrium. To evaluate the scour phenomenon around offshore wind foundations, it is reasonable to apply the wave condition and the LES turbulence model to numerical model applications.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.17 no.1
• /
• pp.137-145
• /
• 2009

The steady-state kinetics of the selective catalytic reduction (SCR) of $NO_X$ with $NH_3$ has been investigated over a commercial ${V_2}{O_5}/TiO_2$ catalyst. In order to account for the influence of transport effects the kinetics are coupled with a fully transient two-phase 1D+1D monolith channel model. The Langmuir-Hinshelwood (L-H) mechanism is adopted to describe the steady-state kinetic behavior of the ${V_2}{O_5}/TiO_2$ catalyst. The reaction rate expressions are based on previously reported papers and are modified to fit the experimental data. The steady-state chemical reaction scheme used in the present mathematical model has been validated extensively with experimental data of selective $NO_X$ reduction efficiency for a wide range of inlet conditions such as space velocity, oxygen concentrations, water concentration, and $NO_2/NO$ ratio. The parametric investigations are performed to examine how the $NH_3$ slip from a SCR $DeNO_X$ catalyst and the conversion of $NO_X$ are affected by the reaction temperature, $NH_3/NO_X$ feed ratio, and space velocity for feed gas compositions with $NO_2/NO_X$ ratios of 0 and 0.5.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.24 no.3
• /
• pp.189-202
• /
• 2012

An Oscillating Water Column (OWC) wave generation system uses the air flow induced by the vertical motion of water column in the air chamber as a driving force of turbine. It is well known that OWC is one of the most efficient devices to harness wave power. This study estimated the air flow velocity from the time variation of the water level fluctuation in the air chamber under regular wave conditions using 3-dimensional numerical irregular wave tank (3D-NIT) model that can simulate the 3-dimensional irregular wave field. The applicability of the 3D-NIT model was validated by comparing numerically predicted air flow velocities with hydraulic experimental results. In addition, the characteristics of air flow frequency spectrum variation due to the incident frequency spectrum change, and the variations of frequency spectrum and wave reflection due to the existence of converter inside the air chamber were discussed. It is found that the phase difference exists in between the air flow velocity and the water level fluctuation inside the air chamber, and the peak frequency of the spectrum in water level fluctuation is amplified by the resonance in the air chamber.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.1
• /
• pp.8-18
• /
• 2015

A numerical simulation for a nonslender BWB UCAV configuration with a rounded leading edge and span of 1.0 m was performed to analyze its aerodynamic characteristics. Numerical results were compared with experimental data obtained at a free stream velocity of 50 m/s and at angles of attack from -4 to $26^{\circ}$. The Reynolds number, based on the mean chord length, is $1.25{\times}106$. 3D multi-block hexahedral grids are used to guarantee good grid quality and to efficiently resolve the boundary layer. Menter's shear stress transport model and two transition models (${\gamma}-Re_{\theta}$ model and ${\gamma}$ model) were used to assess the effect of the laminar/turbulent transition on the flow characteristics. Aerodynamic coefficients, such as drag, lift, and the pitching moment, were compared with experimental data. Drag and lift coefficients of the UCAV were predicted well while the pitching moment coefficient was underpredicted at high angles of attack and influenced strongly by the selected turbulent models. After assessing the pressure distribution, skin friction lines and velocity field around UCAV configuration, it was found that the transition effect should be considered in the prediction of aerodynamic characteristics of vortical flow fields.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2000.10a
• /
• pp.196-203
• /
• 2000

In this study, the evaluation method of diffusion characteristics of Suspended Soli&SS) and the generation limit(source and thick) are investigated, which is significantly affecting on marine examined by construction works such as dredging and reclamation. Dispersion characteristics of SS is examined by hydraulic tests and numerical works in consideration with the Pusan Port. Hydraulic model test was performed in 2-D wave flume to find the limit wave conditon of re-suspension of solid as well as the time dependent characteristics of settlement The results obtainded in the study are as follows; 1) The quantituative evaluation af SS is the basic parameter of marine environmental impact assessment in related with the port development The SS increases as the water content of sea bed solid increases and the density decreases. 2) The sea bed solid in Sinsundai area, Pusan Port has the water content range of 83~157% 3) The ratio of suspension velocity against settlement velocity is about 0.25 and SS concentration converges as the wave heigh. 4) The SS increases 2 time when time step increases 3 time(10 sec to 30 sec) in numerical simulation It means that the effect of the time step should be checked in detail to stable. The diffusion The diffusion coefficient are Affiected senstively in the dispersion process while sea ved friction coefficinet have not strong relation in the simulated area

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.60-60
• /
• 2021

Investigating Backward-Facing Step(BFS) flow is important in that it is a representative case for separation flows in various engineering flow systems. There have been a wide range of experimental, theoretical, and numerical studies to investigate the flow characteristics over BFS, such as flow separation, reattachment length and recirculation zone. However, most of such previous studies were concentrated only on the perpendicular step angle. In this study, several numerical investigations on the flow pattern over BFS with various step angles (10° ~ 90°) and expansion ratios (1.48, 2 and 3.27) under different Reynolds numbers (5000 ~ 64000) were carried out, mainly focused on the reattachment length. The numerical simulations were performed using an open source 3D CFD software, OpenFOAM, in which the velocity profiles and turbulence intensities are calculated by RANS (Reynolds Averaged Navier-Stokes equation) and 3D LES (Large Eddy Simulation) turbulence models. Overall, it shows a good agreement between simulations and the experimental data by Ruck and Makiola (1993). In comparison with the results obtained from RANS and 3D LES, it was shown that 3D LES model can capture much better and more details on the velocity profiles, turbulence intensities, and reattachment length behind the step for relatively low Reynolds number(Re < 11000) cases. However, the simulation results by both of RANS and 3D LES showed very good agreement with the experimental data for the high Reynolds number cases(Re > 11000). For Re > 11000, the reattachment length is no longer dependent on the Reynolds number, and it tends to be nearly constant for the step angles larger than 30°.) Based on the calibrated and validated numerical simulations, several additional numerical simulations were also conducted with higher Reynolds number and another expansion ratio which were not considered in the experiments by Ruck and Makiola (1993).

• Journal of computational fluids engineering
• /
• v.20 no.2
• /
• pp.23-36
• /
• 2015

The $2^{nd}$ generation ice accretion analysis program has been developed and validated for various icing conditions. The essential feature of the $2^{nd}$ generation code lies in its capability of handling general 3-D geometry and improved accuracy. The entire velocity fields are obtained based on Navier-Stokes equations in order to take the massively separated flow field into account. Unlike $1^{st}$ generation code, the droplet trajectories are calculated using Eulerian approach, which is adopted to yield appropriate collection efficiency even in the shadow region. For improved thermodynamic analysis on the surfaces, water film model and modified Messinger model are newly included in the present analysis. The ice shape for a given time step is obtained by considering the exact amount of ice accreted on the surface. Each module of the icing analysis code has been seamlessly integrated on the OpenFOAM platform. The developed code was validated against available experimental data for 2D airfoils and 3D DLR-F4. Due to the lack of experimental data, the computed results of DLR-F4 were compared with those obtained from FENSAP-ICE, which is state-of-the-art 3D icing analysis code. It was clearly shown that the present code produces comparable results to those of FENSAP-ICE, in terms of prediction accuracy and the capability of handling general 3-D geometries.