• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.215-220
• /
• 2008

A numerical method for gas-liquid two-phase flow is applied to solve shock-bubble interaction problems. The present method employs a finite-difference Runge-Kutta method and Roe's flux difference splitting approximation with the MUSCL-TVD scheme. A homogeneous equilibrium cavitation model is used. By this method, a Riemann problem for shock tube was computed for validation. Then, shock-bubble interaction problems between cylindrical bubbles located in the liquid and incident liquid shock wave are computed.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.156-157
• /
• 2003

A simple method is proposed to split the flux vector of the Euler equations by introducing two artificial wave speeds. The direction of wave propagation can be adjusted by these two wave speeds. This idea greatly simplifies the upwinding, and leads to a new family of upwind schemes. Numerical flux function for multi-dimensional Euler equations is formulated for any grid system, structured or unstructured. A remarkable simplicity of the scheme is that it successfully achieves one-sided approximation for all waves without recourse to any matrix operation. Moreover, its accuracy is comparable with the exact Riemann solver. For 1-D Euler equations, the scheme actually surpasses the exact solver in avoiding expansion shocks without any additional entropy fix. The scheme can exactly resolve stationary contact discontinuities, and it is also freed of the carbuncle problem in multi­dimensional computations.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2011.04a
• /
• pp.258-259
• /
• 2011

• International Journal of Aeronautical and Space Sciences
• /
• v.3 no.1
• /
• pp.74-85
• /
• 2002

The numerical experiment has been conducted to investigate the unsteady shock wave reflecting phenomena. The cell-vertex finite-volume, Roe's upwind flux difference splitting method with unstructured grid is implemented to solve unsteady Euler equations. The $4^{th}$-order Runge-Kutta method is applied for time integration. A linear reconstruction of the flux vector using the least-square method is applied to obtain the $2^{nd}$-order accuracy for the spatial derivatives. For a better resolution of the shock wave and slipline, the dynamic grid adaptation technique is adopted. The new concept of grid adaptation technique, which is much simpler than that of conventional techniques, is introduced for the current study. Three error indicators (divergence and curl of velocity, and gradient of density) are used for the grid adaptation procedure. Considering the quality of the solution and the numerical efficiency, the grid adaptation procedure was updated up to $2^{nd}$ level at every 20 time steps. For the convenience of comparison with other experimental and analytical results, the case of interaction between the straight incoming shock wave and a sharp wedge is simulated for various flow conditions. The numerical results show good agreement with other experimental and analytical results, in the shock wave reflecting structure, slipline, and the trajectory of the triple points. Some critical cases show disagreement with the analytical results, but these cases also have been proven to show hysteresis phenomena.

• Journal of Korean Port Research
• /
• v.12 no.2
• /
• pp.281-290
• /
• 1998

The propagation of water waves over irregular bottom bathymetry and around islands involves many process-shoaling, refraction, energy dissipation and diffraction. Numerical model in this study is developed with the mild slope equation to investigate wave transformation in water of varying depth and combined waves and a current. The method used is splitting method and minimax approximation. The numerical method used in this study is Crank-Nicolson scheme in the FDM. This model is applied to Vincent shoal and compared with laboratory experimental data. The results agreed well with laboratory data. Current effect is considered in this study. This model can be used for the estimation of rip current in the slowly varying topography.

• Current Optics and Photonics
• /
• v.8 no.1
• /
• pp.56-64
• /
• 2024

In response to the current situation where general methods cannot effectively compensate for the phase delay of ordinary optical mixers, a multi-layer spatial beam-splitting optical mixer is designed using total reflection triangular prisms and polarization beam splittings. The phase delay is generated by the wave plate, and the mixer can use the existing parallel plates in the structure to individually compensate for the phase of the four output beams. A mixer model is established based on the structure, and the influence of the position and orientation of the optical components on the phase delay is analyzed. The feasibility of the phase compensation method is simulated and analyzed. The results show that the mixer can effectively compensate for the four outputs of the optical mixer over a wide range. The mixer has a compact structure, good performance, and significant advantages in phase error control, production, and tuning, making it suitable for free-space coherent optical communication systems.

• Proceedings of the Optical Society of Korea Conference
• /
• 2000.02a
• /
• pp.124-125
• /
• 2000

We introduce a scheme with a maximized efficiency of detecting atoms passing through an optical standing-wave mode cavity. Consider a standing-wave optical cavity illuminated by a weak probe beam through one of its mirrors where the transmission through the other mirror is monitored by a photodetector. If an atom is put in the cavity, the atom-cavity coupling shifts the resonance frequency of the system via the so-called normal mode splitting, and thereby the transmission power will drop. In fact, this type of atom detection scheme has been used in recent single atom trap experiments In practice, however, the field in a standing-wave mode will have a geometrical structure having nodes and antinodes that when the atom traverses the cavity through one of the nodes, there will be no such effect of atom-field interaction. (omitted)

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.29-35
• /
• 1999

A two-dimensional Euler code based on flux vector splitting scheme has been developed to simulate the behavior of supersonic shock wave over the cylinder. AF+ADI scheme was used for time integration. The sliding multiblock technique was implemented to handle the relative motion of the moving cylinder and the stationary tunnel. The code is validated with a problem of subsonic flow around a Naca-0012 airfoil. The Computation results show complex phenomena of the propagation of shock waves and the reflection as expansion wave at tunnel exit.

• Journal of the Korean Magnetics Society
• /
• v.9 no.6
• /
• pp.271-277
• /
• 1999

Amorphous $Fe_{83}B_9Nb_7Cu_1$ alloy has been studied by M$\"{o}$ssbauer spectroscopy. Revised Vincze method was used and distributions of hyperfine field, isomer shift, and quadrupole line broadening of the sample at various temperatures have been evaluated and Curie temperature and $H_{hf}\;(0)$ were calculated to be 393 K and 231 kOe, respectively. Temperature variation of reduced average hyperfine field shows a flattered curvein comparison with the Brillouin curve for S=1. This behavior can be explained on the basis of Handrich molecular field model, in which the parameter Δ, which is a measure of fluctuation in exchange interactions, is assumed to have the temperature dependence ${Delta}=0.75-0.64{\tau}+0.47{\tau}^2$ where $\tau$ is $T/T_C$. At low temperature, the average hyperfine field can be fitted to $H_{hf}\;(T)=H_{hf}\;(0)\;[1-0.44\;(T/T_C)^{3/2}-0.28(T/T_C)^{5/2}-… ]$, which indicates the presence long wave length spin wave excitations. At temperature near TC, reduced average hyperfine field varies as $1.00\;[1-T/T_C]^{0.39}$. It is also found that half-width of the hyperfine field distribution was 102 kOe (3.29 mm/s) at 13 K and decreased monotonically as temperature increased. Above the Curie temperature, an average quadrupole splitting value of 0.43 mm/s was found. Average line broadening due to quadrupole splitting distribution was 0.31 mm/s at 13 K and decreases monotonically to 0.23 mm/s at 320 K, whereas that due to the isomer shift distribution is 0.1 mm/s at 13 K and 0.072 mm/s at 320 K, which is much smaller than that of both hyperfine field and quadrupole splitting. The temperature dependence of the isomer shift can be fitted within the harmonic approximation to a Deybe model with a Debye temperature ${Theta}_D=424{\pm}5K$.TEX>.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.4
• /
• pp.265-273
• /
• 2014

Wave-body interaction is simulated using a developed code based on the flux-difference splitting scheme for immiscible and incompressible fluids and the hybrid Cartesian/immersed boundary method. A free surface is captured as a moving contact discontinuity within a fluid domain and an approximated Riemann solver is used to estimate the inviscid flux across the discontinuity. Immersed boundary nodes are identified inside an instantaneous fluid domain near a moving body, then dependent variables are reconstructed at those immersed boundary nodes based on interpolation along local normal lines to the boundary. Free surface flows around an oscillating cylinder are simulated and the computed wave elevations are compared with other reported results. The generation of a solitary wave by a moving wave-maker is simulated and the time histories of wave elevations at two different points are compared with other results. The developed code is applied to simulate body motion of an elastically mounted circular cylinder as a solitary wave passes the body. The force acting on an elastically mounted cylinder is compared with the force acting on a fixed cylinder. Grid independency of the computed body motion is established based on a comparison of results using three different-size grids.