• Ocean Systems Engineering
• /
• v.11 no.4
• /
• pp.313-330
• /
• 2021

In this study, the dynamic interaction between an Arctic Spar and drifting level ice is examined in time domain using the newly developed ice-hull-mooring coupled dynamics program. The in-house program, CHARM3D, which is the hull-riser-mooring coupled dynamic simulator is extended by coupling with the open-source discrete element method (DEM) simulator, LIGGGHTS. In the LIGGGHTS module, the parallel-bonding method is implemented to model the level ice using an assembly of multiple bonded spherical particles. As a case study, a spread-moored Artic Spar platform, whose hull surface near waterline is the inverted conical shape, is chosen. To determine the breaking-related DEM parameter (the critical bonding strength), the four-point numerical bending test is used. A series of numerical simulations is systematically performed under the various ice conditions including ice drift velocity, flexural strength, and thickness. Then, the effects of these parameters on the ice force, platform motions, and mooring tensions are discussed. The simulations reveal various features of dynamic interactions between the drifting ice and moored platform for various ice conditions including the novel synchronous resonance at low ice speed. The newly developed simulator is promising and can repeatedly be used for the future design and analysis including ice-floater-mooring coupled dynamics.

• Nuclear Engineering and Technology
• /
• v.56 no.3
• /
• pp.1052-1065
• /
• 2024

Advanced Pressurized Water Reactors (APWRs) and Boiling Water Reactors (BWRs) employ a suppression pool as a heat sink to prevent containment overpressure. Steam can be discharged into the pool through multi-hole spargers or blowdown pipes in both normal and accident conditions. Direct Contact Condensation (DCC) creates sources of momentum and heat. The competition between these two sources determines the development of thermal stratification or mixing of the pool. Thermal stratification is of safety concern as it reduces the cooling capability compared to a completely mixed pool condition. In this work we develop a scaling approach to prediction of the thermal stratification in a water pool induced by steam injection through spargers. Experimental data obtained from large-scale pool tests conducted in the PPOOLEX and PANDA facilities, as well as simulation results obtained using validated codes are used to develop the scaling. Two injection orientations, namely radial injection through multi-hole Sparger Head (SH) and vertical injection through Load Reduction Ring (LRR), are considered. We show that the erosion rate of the cold layer can be estimated using the Richardson number. In this work, scaling laws are proposed to estimate both the (i) transient erosion velocity and (ii) the stable position of the thermocline. These scaling laws are then implemented into a 1D model to simulate the thermal behavior of the pool during steam injection through the sparger.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.3
• /
• pp.418-424
• /
• 2015

The chemical response of energetic materials is analyzed in terms of 1) the thermal decomposition under the thermal stimulus and 2) the reactive flow upon the mechanical impact, both of which give rise to an exothermic thermal runaway or an explosion. The present study aims at building a set of chemical kinetics that can precisely model both thermal and impact initiation of a heavily aluminized cyclotrimethylene-trinitramine (RDX) which contains 35% of aluminum. For a thermal decomposition model, the differential scanning calorimetry (DSC) measurement is used together with the Friedman isoconversional method for defining the frequency factor and activation energy in the form of Arrhenius rate law that are extracted from the evolution of product mass fraction. As for modelling the impact response, a series of unconfined rate stick data are used to construct the size effect curve which represents the relationship between detonation velocity and inverse radius of the sample. For validation of the modeled results, a cook-off test and a pressure chamber test are used to compare the predicted chemical response of the aluminized RDX that is either thermally or mechanically loaded.

• Journal of the Korean Institute of Gas
• /
• v.20 no.1
• /
• pp.7-12
• /
• 2016

Natural gas liquefaction process is essential to transport natural gas for long distances. Lots of compressors in this process are needed and the energy for these compressors can be supplied by drivers. Total driver cost can be changed by selecting various drivers. This study focused on the minimization of the driver cost to provide the energy to the compressors. Moreover, scenarios, extracting velocity is changed during whole operating period, are set with considering gas well capacity. The mathematical model was established by considering trade off relationship between the capital cost and the operating cost of the turbines. The model also considers the life time of the driver equipments. As the result, the driver cost of the optimized case was reduced by 6.4% than the base case.

• Nuclear Engineering and Technology
• /
• v.22 no.1
• /
• pp.20-28
• /
• 1990

In this study the migration experiment using packed column with crushed tuff was conducted as a basic research to develop migration model of radionuclides through geologic media. The main emphasis was put on evaluating the validity of migration models. For this, two models were introduced: one is the model which is based on the assumption of instantaneous equilibrium reaction and the other the model based on kinetic process such as intraparticle diffusion. The coefficient of hydrodynamic dispersion in packed column was determined using iodine as nonsorbing tracer. The hydrodynamic dispersion coefficient, D$_{L}$ was shown to be 0.11$\times$10$^{-2}$ $\textrm{cm}^2$/min under the condition of the column porosity of 0.483 and the average water velocity of 0.915$\times$10$^{-2}$ cm/min. The distribution coefficient, Kd of Cs-137 on crushed tuff was 11.3 cc/g at the concentration of 2$\times$10$^{-6}$ M and the temperature of 2$0^{\circ}C$. The breakthrough curve of Cs-137 through packed column was shown to have an asymmetric curve in which long trailing tail appears at the end part of the curve. The results obtained from the comparison of introduced models with experimental data indicated that the mass transfer model with intraparticle diffusion as rate-controlling step simulated the behaviors of Cs-137 migration more adequately, when compared with the bulk reaction model in which the assumption of instantaneous equilibrium reaction was maded. Consequently, the intraparticle diffusion was found to be an important factor in the migration of Cs-137 through packed column.n.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.6
• /
• pp.1753-1764
• /
• 2014

Numerical flow simulation models in the riverine environments have been widely utilized for analyzing flow dynamics in various degrees in researches and practical applications. However, most of the simulated results have been validated based on the data from indoor experimental models or very limited in-situ measurements. Therefore, it has been required to more accurately validate the performance of the numerical models in terms of the detailed field observations. In particular, it was also hard to validate the performances of the existing numerical models in the real meandered river channels that encompass more sophisticated flow and geometric structures. Recently, advancements of the modern flow measuring instrumentations such as acoustic Doppler current profilers (ADCPs) enabled us to efficiently acquire the detailed flow field in the broad range of river channels, thus that it became to be possible to accurately validate any numerical models with the field observations. In this study, based on the detailed flow measurements in a actual meandered river channel using ADCP, we validated FaSTMECH model in iRIC in terms of water surface elevation, which is relatively new but began to get highlighted in the research areas. As the validation site, a meandering channel in River Experiment Center of KICT was chosen, which has 6.5 m of width, 0.38m of flow depth, 1.54 m3/s of flow discharge, 0.61 m/s of mean flow velocity, and 1.2 of sinuosity. As results, whereas the FaSTMECH precisely simulated water surface elevation, simulated velocity field in the bend did not match well with ADCP dataset.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.28 no.3
• /
• pp.146-159
• /
• 2016

This study used Navier-Stokes Solver(LES-WASS-2D) for analyzing hydrodynamic characteristics with high order in order to analyze self-burial mechanism of pipeline with spoiler under steady flow. For the validity and effectiveness of numerical model used, it was compared and analyzed with the experiment to show flow characteristics around the pipeline with and without the spoiler. And the hydraulic(flow, vortex, and pressure) and force characteristics were numerically analyzed around the pipeline according to the incident velocity, and shape and arrangement of spoiler. Primarily, if the spoiler is attached to the pipeline, the projected area is increased resulting in higher flow velocity toward the back and strong vortex caused by wake stream in the back. Secondly, the spoiler causes vertically asymmetric flow and vorticity fields and thus asymmetric pressure field. It increases the asymmetry of force on the pipe and thus develops large downward fluid force. Both of them are the causes of selfburying of the pipeline with spoiler.

• Journal of Energy Engineering
• /
• v.8 no.2
• /
• pp.224-232
• /
• 1999

Latent heat storage system using micro-encapsuled phase change material is effective method for floor heating of house and building. The temperature profile in capsule block and flow rate of hot water are important parameters for the development of heat storage system. In the present study, a mathematical model based on 3-D, non-steady state, Navier-Stokes equations, scalar conservation equations and turbulence model ($\kappa$-$\varepsilon$), is used to predict the temperature profiles in capsule and the velocity vectors in hot water pipe. The multi-block grids and fine grids embedding are used to join the circle in hot water pipe and square in capsule block. The phase change process of the capsule is quite complex not only because the size of phase change material is very small, but also because phase change material is mixed with the cement to form thermal storage block. In calculation, it's assumed that the phenomena of phase change is limited only the thermal properties of phase change material and the change of boundary is not happened in capsule. The purpose of this study is to calculate the temperature profiles in capsule block and velocity vectors in hot water pipe using the numerical calculation. Two kinds of thermal boundary condition were considered, the first (case 1) is the adiabatic condition for the both outside surfaces of the wall, the second (case 2) is the case in which one surface is natural convection with atmosphere and another surface is adaibatic. Calculation results are shown that the temperature profile in capsule block for case 1 is higher than that for case 2 due to less heat loss in adaibatic surface. Specially, in the domain of near Y=0, the difference of temperature is greater in case 1 than in case 2. The detailed experimental data of capsule block on the temperature profile and the thermal properties such as specific heat and coefficient of heat transfer with the various temperature are required to predict more exact phenomena of heat transfer.

• Journal of the Korean earth science society
• /
• v.32 no.6
• /
• pp.640-653
• /
• 2011

We simulate the propagation of earthquake waves in the continental margin of Antarctica using the elastic wave modeling algorithm, which is modified to be suitable for acoustic-elastic coupled media and earthquake source. To simulate the various types of earthquake source, the staggered-grid finite-difference method, which is composed of velocity-stress formulae, can be more appropriate to use than the conventional, displacement-based, finite-difference method. We simulate the elastic wave propagation generated by earthquakes combining 3D staggered-grid finite-difference algorithm composed of displacement-velocity-stress formulae with double couple mechanisms for earthquake source. Through numerical tests for left-lateral strike-slip fault, normal fault and reverse fault, we could confirm that the first arrival of P waves at the surface is in a good agreement with the theoretically-predicted results based on the focal mechanism of an earthquake. Numerical results for a model made after the subduction zone in the continental margin of Antarctica showed that earthquake waves, generated by the reverse fault and propagating through the continental crust, the oceanic crust and the ocean, are accurately described.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.50 no.8
• /
• pp.531-540
• /
• 2022

In this paper, we present an algorithm that identifies artillery type and rapidly predicts the impact point based on the IMM filter. The ballistic trajectory equation is used as a system model, and three models with different ballistic coefficient values are used. Acceleration was divided into three components of gravity, air resistance, and lift. And lift acceleration was added as a new state variable. The kinematic condition that the velocity vector and lift acceleration are perpendicular was used as a pseudo-measurement value. The impact point was predicted based on the state variable estimated through the IMM filter and the ballistic coefficient of the model with the highest mode probability. Instead of the commonly used Runge-Kutta numerical integration for impact point prediction, a semi-analytic method was used to predict impact point with a small amount of calculation. Finally, a state variable initialization method using the least-square method was proposed. An integrated algorithm including artillery type identification, impact point prediction and initialization was presented, and the validity of the proposed method was verified through simulation.