• Journal of the Korea Computer Graphics Society
• /
• v.16 no.1
• /
• pp.9-20
• /
• 2010

In order to visualize participating media in 3D space, they usually calculate the incoming radiance by subdividing the ray path into small subintervals, and accumulating their respective light energy due to direct illumination, scattering, absorption, and emission. Among these light phenomena, scattering behaves in very complicated manner in 3D space, often requiring a great deal of simulation efforts. To effectively simulate the light scattering effect, several approximation techniques have been proposed. Volume photon mapping takes a simple approach where the light scattering phenomenon is represented in volume photon map through a stochastic simulation, and the stored information is explored in the rendering stage. While effective, this method has a problem that the number of necessary photons increases very fast when a higher variance reduction is needed. In an attempt to resolve such problem, we propose a different approach for rendering particle-based volume data where kernel smoothing, one of several density estimation methods, is explored to represent and reconstruct the light in-scattering effect. The effectiveness of the presented technique is demonstrated with several examples of volume data.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.32 no.4B
• /
• pp.243-251
• /
• 2012

To investigate the characteristics of rip current occurring at Haeundae beach, observations obtained from a buoy and a CCTV were analyzed and numerical experiments were conducted. During observed rip-current events, the CCTV images showed that a couple of wave-trains, which are close to regular waves with slightly different directions, propagated to the beach, and wavelet analyses of data from the buoy showed very narrow-banded spectra with a peak frequency. From the evidences, it was inferred that a known mechanism of generating rip current due to the nodal line area of honeycomb-patterned wave crest was one of the significant factors of rip current occurrences of Haeundae beach. The mechanism has been explained by the following: When two wave-trains with slightly different directions propagate to a beach, wave crests of the incident wave-trains form honeycomb pattern due to nonlinear interaction. The nodal lines of honeycomb pattern are developed in the cross-shore direction. And longshore currents flow toward the nodal line area which has very low wave energy. Consequently their mass flux is expelled through the area toward the sea direction. To confirm the generation, numerical experiments were performed using a nonlinear Boussinesq equation model. In the cases with two incident wave-trains with slightly different directions and with a monochromatic wave propagating over a submerged shoal, it was seen that the honeycomb pattern of wave crests was well developed, and thus rip currents were evolved along the nodal lines.

• The Journal of the Acoustical Society of Korea
• /
• v.36 no.5
• /
• pp.314-320
• /
• 2017

A pressure relief valve is generally used to prevent piping systems from being broken due to high pressure gas flows. However, the sudden pressure drop caused by the pressure relief valve produces high acoustic energy which propagates in the form of compressible acoustic waves in the pipe and sometimes causes severe vibration of the pipe structure, thereby resulting in its failure. In this study, internal aerodynamic noise due to valve flow is estimated for a simple contraction-expansion pipe by combining the LES (Large-Eddy Simulation) technique with the wavenumber-frequency analysis, which allows the decomposition of fluctuating pressure into incompressible hydrodynamic pressure and compressible acoustic pressure. In order to increase the convergence, the steady Reynolds-Averaged Navier-Stokes equations are numerically solved. And then, for the unsteady flow analysis with high accuracy, the unsteady LES is performed with the steady result as the initial value. The wavenumber-frequency analysis is finally performed using the unsteady flow simulation results. The wavenumber-frequency analysis is shown to separate the compressible pressure fluctuation in the flow field from the incompressible one. This result can provide the accurate information for the source causing so-called acoustic-induced-vibration of a piping system.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.4
• /
• pp.401-412
• /
• 2008

Simulations of three-dimensional numerical wave tank are performed to investigate wave force acting on a large cylindrical structure and consequent wave deformation, which are induced by bore after breaking waves. The numerical model is based on the three-dimensional Navier-Stokes equations with a finite-difference method combined with a volume of fluid(VOF) method, which is capable of tracking the complex free surface, including wave breaking. In order to promote wave breaking of the incident wave, the approach slope was built seaward of the structure with a constant slope and a large cylindrical structure was installed on a flat bed. The incident waves were broken on the approach slope or flat bed by its wave height. In the present study, all waves acting on the large cylindrical structure were limited to breaking bore after wave breaking. The effects of the position of the structure and the incident wave height on the wave force and wave transformations were mainly investigated with the concern of wave breaking. Further, the relations between the variation of wave energy by wave propagation after wave breaking and wave force acting on the structure were discussed to give the understanding of the full-linear wave-structure interactions in three-dimensional wave fields.

• Korean Journal of Food Science and Technology
• /
• v.24 no.3
• /
• pp.256-260
• /
• 1992

This study was conducted to calculate and predict the shelf-life of meat sausage through physicochemical and biological analysis. Judging from physico-chemical components, it can be found out the most effective indicator in meat sausage. And also, the materials used for packaging is cellulose #23 in inner-packaging and CN/HDPE laminating film in outer-packaging. The changes of the most effective indicator were discussed through the method of kinetic analysis. Judging from physico-chemical components, VBN was the most available component in quality judgement of meat sausage and their upper limiting contents were 20 mg%. It is possible to calculate and predict the shelf-life of meat sausage through the regression equation and $Q_{10}$ value. As a result, the shelf-life prediction was $58{\sim}63$ days at $10^{\circ}C$, $47{\sim}51$ days at $20^{\circ}C$ and 26 days at $40^{\circ}C$, respectively, but the difference between two methods showed about $4{\sim}5$ days. $Q_{10}$ value on the changes contents was 1.35 at acceralated temperature $40^{\circ}C$. The reaction rate of VBN contents could be interpreted as a first order reaction that divided into 2 periods with different reaction rate constants. The corresponding Arrhenius activation energies were 2.959 Kcal/mole and 3.632 Kcal/mole, respectively.

• International Journal of Highway Engineering
• /
• v.7 no.1 s.23
• /
• pp.49-61
• /
• 2005

The dynamic modulus of hot mix asphalt can be determined according to the different combinations of testing temperature and loading frequency. The superposition rule is adapted to get the master curve of dynamic modulus for each hot mix asphalt. There are couple of different methods to get the shift factor which is a key for making the master curve. In this paper, Arrehnius, 2002 AASHTO, and experimental method was employed to get the master curve. Evaluation of dynamic modulus for 25mm base course of hot mix asphalt with granite aggregate and two asphalt binders(AP-3 and AP-5) was carried out. Superpave Level 1 Mix Design with gyratory compactor was adopted to determine the optimum asphalt binder content(OAC) and the measured ranges of OAC were between 4.1% and 4.4%. UTM was used for laboratory test. The dynamic modulus and phase angle were determined by testing on UTM, with 5 different testing temperature(-10, 5, 20, 40, & $55^{\circ}C$) and 5 different loading frequencies(0.05, 0.1, 1, 10, 25 Hz). Using the measured dynamic modulus and phase angle, the input parameters of Sigmoidal function equation to represent the master curve were determined and these will be adopted in FEM analysis for asphalt pavements. The shift factor and activation energy for determination of master curve were calculated.

• Membrane Journal
• /
• v.20 no.1
• /
• pp.29-39
• /
• 2010

This study focuses on the modeling of DMFC to predict the characteristics and to improve its performance. This modeling requires deep understanding of the design and operating parameters that influence on the cell potential. Furthermore, the knowledge with reference to electrochemistry, transport phenomena and fluid dynamics should be employed for the duration of mathematical description of the given process. Considering the fact that MEA is the nucleus of DMFC, special attention was made to the development of mathematical model of MEA. Multiscale modeling is comprised of process modeling as well as a computational fluid dynamics (CFD) modeling. The CFD packages and process simulation tools are used in simulating the steady-state process. The process simulation tool calculates theelectrochemical kinetics as well as the change of fractions, and at the same time, CFD calculates various balance equations. The integrated simulation with multiscal modeling explains experimental observations of transparent DMFC.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.34 no.10
• /
• pp.1427-1435
• /
• 2010

A mixed finite element model was developed using the classical plate theory to analyze the nonlinear bending of a plate. The appropriate weight functions for the constraints integrated over the domain were determined by the Lagrange multiplier method by using the principle of minimum virtual energy; which provides the constitutive relations between force-like variables and strains. All of detail terms of element wise coefficient matrices and associate tangent matrices to be used in the Newton iterative method are presented. Then, the linear solutions of the current model and those of the traditional displacement model under the SS (simple support) boundary conditions were compared with the existing analytical solution. The post-processed images of the nonlinear results of the force-like variables are presented to show the continuity of the solutions at the joint of the element boundaries. Finally, the converged nonlinear finite element solutions of the current model are compared with those of existing traditional displacement model.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.55-55
• /
• 2011

연안 및 해안공학의 발달과 더불어 부유식방파제의 기능적 효율성이 중요시 되고 있다. 흔히 사용되어오던 착저식방파제는 설치에 많은 시간과 경비가 소요되고 환경 및 생태계에 많은 변화를 줄 수 있으며, 설치 예정지의 수리학적 특성 등의 여건에 많은 제약을 받는 단점이 있다. 부유식방파제는 일본 등의 선진국을 중심으로 활용이 잦아지고 있는 방파제로서 수면 위에 설치되기 때문에 수중 생태계에 미치는 영향이 적은 친환경방파제이다. 또한 기존에 시공된 중력식방파제와는 달리 수심에 제한을 덜 받고, 공사기간이 짧기 때문에 경제적이다. 실제 시공사례로는 2007년 마산 원전항에 완공된 부유식방파제가 대표적이며, 지금까지도 부유식방파제에 대한 여러 연구자들의 관심이 증가하고 있는 추세이다. 방파제뿐만 아니라 우리나라처럼 국토의 면적이 작은 지역에서 증가하는 해상물동량을 소화하기 위해서 부유식방파제 등을 이용한 항만의 시공이 필요한 실정이다. 이러한 부유식방파제의 분석적인 측면에 있어서 수치해석은 파랑과 구조물의 상호작용을 해석하는 데 한계가 있으며, 부유식방파제 단면형상을 정확하게 재현할 수 없으므로, 수리모형실험을 통한 부유식방파제의 연구가 필요할 것으로 판단된다. 최근 기술의 발달로 인한 유동장 해명이 가능해 졌으며, PIV(Particle image velocimetry) 및 LDV시스템은 다양한 분야에서 응용되고 있다. 특히, LDV시스템은 측정하려는 한 지점에 대하여 레이저 빔을 단면(Cross-section)으로 만들고 입자의 산란광을 후방산란(Back scatter)으로 받아서 도플러 효과를 이용, 속도에 대한 주파수를 획득하며, 유속을 측정하는 장비로 매우 높은 정확도와 비접촉식 이라는 장점을 가지고 있다. 또한, PIV 시스템에 비하여 측정시간이 오래 걸리는 반면 데이터를 가공하지 않고 활용할 만큼 높은 정확성을 가지고 있다. 본 연구에서는 수리모형실험을 통하여 단독형, 2열형 및 3열형 부유식방파제의 형상, 흘수 및 거리를 변화시키며 유동장을 수집하였으며, 방파성능에 따른 와의 생성 및 소멸시점에서의 파랑변형과의 관계를 분석하였다. 방파제의 형상과 흘수를 달리하여 수리모형실험을 수행하였으며, 와류의 상관관계를 분석하였다. 또한, 연직 2차원 Navier-Stokes 방정식 모형을 이용하여 수치모형실험을 수행하였으며, 수치모형실험 결과와 수리모형실험 결과를 비교 분석하였다. 후방방파제에서 발생되는 파랑은 입사파의 주기가 길어질수록 상대적으로 커지는 현상을 보였으며, 흘수심이 깊어질수록 전방방파제 입사 면에서 자유 수면이 높게 관측되는 결과를 보였다. 또한, 비교적 장주기파랑에 해당하는 입사파랑의 경우 전달파고비 산정에 있어서 설계기준인 0.5를 대다수 초과하는 반면, 3열형 구조에서는 대부분이 0.5이하로 상당히 높은 방파성능 결과를 나타내었다.

• Journal of the Korean Institute of Gas
• /
• v.16 no.4
• /
• pp.38-43
• /
• 2012

The purpose of this study is to carry out the consequence analysis of an accident related to the release of Hydrogen chloride occurred in a petrochemical company in Korea and suggest the measures to prevent similar accidents from happening again. The total amount released through the safety valve of HCl Column was calculated based on the rated capacity of the safety valve, the ideal gas equation and mechanical energy balance, respectively. As a result of the calculation, we found that the amount of HCl released through the safety valve was at least 76.8 kg. Also, we predicted the dispersion concentration at the position of the injured workers(more than 350 m away from the accident location) using simulation programs such as PHAST. The results of ALOHA and K-CARM are 304 ppm and 1,700 ppm respectively. However, PHAST calculation indicated that the concentration is less than l ppm. From these results, we can understand that workers were injured by HCl gas released from the safe valve and the concentration of gas might be less than 1 ppm. Also, it is important for toxic gases such as HCl to be vented safely to the atmosphere after scrubbing.