• Fire Science and Engineering
• /
• v.10 no.3
• /
• pp.10-18
• /
• 1996

In foreign country such as USA and Japan, considerable research has been done regarding the spread of smoke in room of fire involvement by using computer. but, in our country it has not been, so, this paper presents a detailed qualitative description of phenomena which occure during typical fire scenarios through numberical analysis. The governing equations are solved by using FVM method with non-staggered grid. The SIMPLE method for pressure-velocity couple and power-law scheme for convention terms are used. It shows that a plume is formed, hot plume gases impinge on the ceiling and they spread across it. then, it eventually reaches the bounding walls of the enclosure. It takes 20s for smoke to fill the enclosure.

• Ocean Systems Engineering
• /
• v.1 no.1
• /
• pp.31-57
• /
• 2011

Impact pressure due to sloshing is of great concern for the ship owners, designers and builders of the LNG carriers regarding the safety of LNG containment system and hull structure. Sloshing of LNG in partially filled tank has been an active area of research with numerous experimental and numerical investigations over the past decade. In order to accurately predict the sloshing impact load, a new numerical method was developed for accurate resolution of violent sloshing flow inside a three-dimensional LNG tank including wave breaking, jet formation, gas entrapping and liquid-gas interaction. The sloshing flow inside a membrane-type LNG tank is simulated numerically using the Finite-Analytic Navier-Stokes (FANS) method. The governing equations for two-phase air and water flows are formulated in curvilinear coordinate system and discretized using the finite-analytic method on a non-staggered grid. Simulations were performed for LNG tank in transverse and longitudinal motions including horizontal, vertical, and rotational motions. The predicted impact pressures were compared with the corresponding experimental data. The validation results clearly illustrate the capability of the present two-phase FANS method for accurate prediction of impact pressure in sloshing LNG tank including violent free surface motion, three-dimensional instability and air trapping effects.

• Nuclear Engineering and Technology
• /
• v.41 no.10
• /
• pp.1347-1360
• /
• 2009

A multi-dimensional component for the thermal-hydraulic system analysis code, MARS, was developed for a more realistic three-dimensional analysis of nuclear systems. A three-dimensional and two-fluid model for a two-phase flow in Cartesian and cylindrical coordinates was employed. The governing equations and physical constitutive relationships were extended from those of a one-dimensional version. The numerical solution method adopted a semi-implicit and finite-difference method based on a staggered-grid mesh and a donor-cell scheme. The relevant length scale was very coarse compared to commercial computational fluid dynamics tools. Thus a simple Prandtl's mixing length turbulence model was applied to interpret the turbulent induced momentum and energy diffusivity. Non drag interfacial forces were not considered as in the general nuclear system codes. Several conceptual cases with analytic solutions were chosen and analyzed to assess the fundamental terms. RPI air-water and UPTF 7 tests were simulated and compared to the experimental data. The simulation results for the RPI air-water two-phase flow experiment showed good agreement with the measured void fraction. The simulation results for the UPTF downcomer test 7 were compared to the experiment data and the results from other multi-dimensional system codes for the ECC delivery flow.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.2
• /
• pp.188-195
• /
• 2010

Green water load is an important parameter to be considered in designing a modern ship or offshore structures like FPSO and FSRU. In this research, a numerical simulation method for green water phenomenon is introduced. The Navier-Stokes equations and the continuity equation are used as governing equations. The equations are calculated using Finite Difference Method(FDM) in rectangular staggered grid system. To increase the numerical accuracy near the body, the Cartesian cut cell method is employed. The nonlinear free-surface during green water incident is defined by Marker-density method. The green waters on a box in regular waves are simulated. The simulation results are compared with other experimental and computational results for verification. To check the applicability to moving ship, the green water of the ship which is towed by uniform force in regular wave, is simulated. The ship is set free to heave and to surge.

• Journal of Korea Water Resources Association
• /
• v.38 no.11
• /
• pp.973-984
• /
• 2005

A three-dimensional thermal hydrodynamic model is developed for carrying out unsteady simulation of the selective withdrawal of the stably stratified flow in a geometrically complex, natural reservoir The governing equations are discretized on a non-staggered grid using a second-order accurate, finite-volume scheme. The numerical model is validated by applying it to simulate three-dimensional, turbulent, stratified, shear-layer flow case. The numerical predictions appear to capture reasonably well the general shape of velocity and temperature profiles observed in the laboratory experiments, while significant overestimation of the magnitude of velocity profiles is observed in the application to the flow in a natural reservoir. The physics of selective withdrawal as emerge from the numerical simulations are also discussed.

• Ocean Systems Engineering
• /
• v.10 no.2
• /
• pp.201-226
• /
• 2020

A systematic numerical comparative study of the performance of semicircular and rectangular submerged breakwaters interacting with solitary waves is the basis of this paper. To accomplish this task, Nwogu's extended Boussinesq model equations are employed to simulate the interaction of the wave with breakwaters. The finite difference technique has been used to discretize the spatial terms while a fourth-order predictor-corrector method is employed for time discretization in our numerical model. The proposed computational scheme uses a staggered-grid system where the first-order spatial derivatives have been discretized with fourth-order accuracy. For validation purposes, five test cases are considered and numerical results have been successfully compared with the existing analytical and experimental results. The performances of the rectangular and semicircular breakwaters have been examined in terms of the wave reflection, transmission, and dissipation coefficients (RTD coefficients) denoted by K_R, K_T, K_D. The latter coefficient K_D emerges due to the non-energy conserving K_R and K_T. Our computational results and graphical illustrations show that the rectangular breakwater has higher reflection coefficients than semicircular breakwater for a fixed crest height, but as the wave height increases, the two reflection coefficients approach each other. un the other hand, the rectangular breakwater has larger dissipation coefficients compared to that of the semicircular breakwater and the difference between them increases as the height of the crest increases. However, the transmission coefficient for the semicircular breakwater is greater than that of the rectangular breakwater and the difference in their transmission coefficients increases with the crest height. Quantitatively, for rectangular breakwaters the reflection coefficients K_R are 5-15% higher while the diffusion coefficients K_D are 3-23% higher than that for the semicircular breakwaters, respectively. The transmission coefficients K_T for rectangular breakwater shows the better performance up to 2.47% than that for the semicircular breakwaters. Based on our computational results, one may conclude that the rectangular breakwater has a better overall performance than the semicircular breakwater. Although the model equations are non-dissipative, the non-energy conserving transmission and reflection coefficients due to wave-breakwater interactions lead to dissipation type contribution.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.5
• /
• pp.335-344
• /
• 2011

A curvilinear non-hydrostatic free surface model is developed to investigate nonlinear wave interactions in a circular channel. The proposed model solves the unsteady Navier-Stokes equations in a three-dimensional domain with a pressure correction method, which is one of fractional step methods. A hybrid staggered-grid layout in the vertical direction is implemented, which renders relatively simple resulting pressure equation as well as free surface closure. Numerical accuracy with respect to wave nonlinearity is tested against the fifth-order Stokes solution in a two-dimensional numerical wave tank. Numerical applications center on the evolution of nonlinear waves including diffraction and reflection affected by the curvature of side wall in a circular channel comparing with linear waves. Except for a highly nonlinear bichrmatic wave, the model's results are in good agreement with superimposed analytical solution that neglects nonlinear effects. Through the numerical simulation of the highly nonlinear bichramatic wave, the model shows its capability to investigate the evolution of nonlinear wave groups in a circular channel.

• 한국지구물리탐사학회:학술대회논문집
• /
• 2006.06a
• /
• pp.181-186
• /
• 2006

An imaging method of seismic sources using time-reversal wave propagation is presented. The method is based on the time-reversal invariance and the spatial reciprocity of the wave equation. Time-reversal wave propagation has been used to image anomalous features of a midium in medical imaging, non destructive testing and waveform tomography. Seismogram is the record whose energy is propagated from the seismic source. If time-reversed seismogram propagates back into the medium, seismic energy is concentrated at the origin time of the event and at the source location. In this work, a staggered-grid finite-difference method of the elastic wave equation is parallelized for 3-D wave propagation simulation. With numerical experiments, we show that the time-reversal imaging will enable us to explore the spatio-temporal history of complex earthquake.