• Proceedings of the KIEE Conference
• /
• summer
• /
• pp.666-670
• /
• 2004

This paper presents a new approach for determining an accurate incentive levels of Direct Load Control (DLC) program using sequential Monte Carlo Simulation (MCS) techniques. The economic analysis of DLC resources needs to identify the hourly-by-hourly expected energy-not-served resulting from the random outage characteristics of generators as well as to reflect the availability and duration of DLC resources, which results the computational explosion. Therefore, the conventional methods are based on the scenario approaches to reduce the computation time as well as to avoid the complexity of economic studies. In this paper, we have developed a new technique based on the sequential MCS to evaluate the required expected load control amount in each hour and to decide the incentive level satisfying the economic constraints. And also the proposed approach has been considered multi-state as well as two-state of the generating units. In addition, we have applied the variance reduction technique to enhance the efficiency of the simulation. To show the efficiency and effectiveness of the suggested method the numerical studies have been performed for the modified IEEE reliability test system.

• Journal of computational fluids engineering
• /
• v.21 no.4
• /
• pp.19-27
• /
• 2016

The flow analyses through a porous hydraulic fractures is among the most important tasks in recently developed shale reservoirs but is rendered difficult by non-Darcy effects and geometric changes in the hydraulic fractures during production. In this study, several Computational Fluid Dynamics(CFD) models of hydraulic fractures, with a simple shape such as that of parallel plates, filled with proppants were built. Direct Numerical Simulation(DNS) analyses were then carried out to examine the flow loss characteristics of the fractures. The hydraulic diameters for the simulation models were calculated using the DNS results, and then they were compared with the results from Kozeny's definition of hydraulic diameter which is most widely used in the flow analysis field. Also, the characteristic parameters based on both hydraulic diameters were estimated for the investigation of the flow loss variation features. Consequently, it was checked in this study that the hydraulic diameter based on Kozeny's definition is not accordant to the results from the DNS analyses, and the case using the CFD results exhibits f Re robustness like general pipe flows, whereas the other case using Kozeny's definition doesn't. Ultimately, it is expected that discoveries reported in this study would help further porous flow analyses such as hydraulic fracture flows.

• Nuclear Engineering and Technology
• /
• v.48 no.4
• /
• pp.859-869
• /
• 2016

Component-scale modeling of boiling is predominantly based on the Eulerian-Eulerian two-fluid approach. Within this framework, wall boiling is accounted for via the Rensselaer Polytechnic Institute (RPI) model and, within this model, the bubble is characterized using three main parameters: departure diameter (D), nucleation site density (N), and departure frequency (f). Typically, the magnitudes of these three parameters are obtained from empirical correlations. However, in recent years, efforts have been directed toward mechanistic modeling of the boiling process. Of the three parameters mentioned above, the departure diameter (D) is least affected by the intrinsic uncertainties of the nucleate boiling process. This feature, along with its prominence within the RPI boiling model, has made it the primary candidate for mechanistic modeling ventures. Mechanistic modeling of D is mostly carried out through solving of force balance equations on the bubble. Forces incorporated in these equations are formulated as functions of the radius of the bubble and have been developed for, and applied to, low-pressure conditions only. Conversely, for high-pressure conditions, no mechanistic information is available regarding the growth rates of bubbles and the forces acting on them. In this study, we use direct numerical simulation coupled with an interface tracking method to simulate bubble growth under high (up to 45 bar) pressure, to obtain the kind of mechanistic information required for an RPI-type approach. In this study, we compare the resulting bubble growth rate curves with predictions made with existing experimental data.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.1
• /
• pp.37-50
• /
• 2016

A large eddy simulation (LES) of a turbulent channel flow is performed by using the third order low-storage Runge-Kutta method in time and second order finite difference formulation in space with staggered grid at a Reynolds number, $Re_{\tau}=590$ based on the channel half width, ${\delta}$ and wall shear velocity, $u_{\tau}$. To reduce the calculation cost of LES, algebraic wall model (AWM) is applied to approximate the near-wall region. The computation is performed in a domain of $2{\pi}{\delta}{\times}2{\delta}{\times}{\pi}{\delta}$ with $32{\times}20{\times}32$ grid points. Standard Smagorinsky model is used for subgrid-scale (SGS) modeling. Essential turbulence statistics of the flow field are computed and compared with Direct Numerical Simulation (DNS) data and LES data using no wall model. Agreements as well as discrepancies are discussed. The flow structures in the computed flow field have also been discussed and compared with LES data using no wall model.

• Journal of Korean Society on Water Environment
• /
• v.25 no.2
• /
• pp.181-192
• /
• 2009

In this study, the hydrological properties of permeable pavement were analyzed by the experiment and the numerical simulation. The numerical model used was a modified SWMM especially for considering the hydrological response of permeable pavement. The parameters of modified SWMM were revised by the experimental results, and then the practicability was evaluated through the comparison of the experimental and numerical simulation results. In the experiments, three different rainfall intensities such as 65 mm/hr, 90 mm/hr, 95 mm/hr were supplied for 4 hrs, and the hydraulic properties including surface outflow, subsurface outflow, ground water level, soil water contents were measured for 10 hrs. The results showed rainfall intensity effected directly on surface outflow volume and subsurface outflow volume was more effected by ground water level than rainfall intensity. The ground water level and the soil water contents were under estimated as compared with the experimental data except the portion of occurring direct runoff. The surface and subsurface outflow discharge were simulated very well in comparison with the experimental data. Consequently, the modified SWMM could be used very effectively to evaluate the hydrological property of permeable pavement.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.33 no.6
• /
• pp.505-510
• /
• 2013

An encircling send-receive type pulsed eddy current (PEC) probe is designed for use in aluminum tube inspection. When bare receive coils located away from the exciter were used, the peak time of the signal did not change although the distance from the exciter increased. This is because the magnetic flux from the exciter coil directly affects the receive coil signal. Therefore, in this work, both the exciter and the sensor coils were shielded in order to reduce the influence of direct flux from the exciter coil. Numerical simulation with the designed shielded encircling PEC probe showed the corresponding increase of the peak time as the sensor distance increased. Ferrite and carbon steel shields were compared and results of the ferrite shielding showed a slightly stronger peak value and a quicker peak time than those of the carbon steel shielding. Simulation results showed that the peak value increased as the defect size (such as depth and length) increased regardless of the sensor location. To decide a proper sensor location, the sensitivity of the peak value to defect size variation was investigated and found that the normalized peak value was more sensitive to defect size variation when the sensor was located closer to the exciter.

• Structural Engineering and Mechanics
• /
• v.30 no.2
• /
• pp.177-190
• /
• 2008

Material failure behavior is generally dependent on loading rate. Especially in brittle and quasi-brittle materials, rate dependent material behavior can be significant. Empirical formulations are often used to predict the rate dependency, but such methods depend on extensive experimental works and are limited by practical constraints of physical testing. Numerical simulation can be an effective means for extracting knowledge about rate dependent behavior and for complementing the results obtained by testing. In this paper, the failure behavior of a brittle material under different loading rates is simulated by molecular dynamics analysis. A notched specimen is modeled by sub-million particles with a normalization scheme. Lennard-Jones potential is used to describe the interparticle force. Numerical simulations are performed with six different loading rates in a direct tensile test, where the loading velocity is normalized to the ratio of the pseudo-sonic speed. As a consequence, dynamic features are achieved from the numerical experiments. Remarkable failure characteristics, such as crack surface interaction/crack arrest, branching, and void nucleation, vary in case of the six loading cases. These characteristics are interpreted by the energy concept approach. This study provides insight into the change in dynamic failure mechanism under different loading rates.

• Ocean Systems Engineering
• /
• v.5 no.2
• /
• pp.109-123
• /
• 2015

In the present study, the coupled dynamic response of a Submerged Floating Tunnel (SFT) and mooring lines under regular waves is solved by using two independent numerical simulation methods, OrcaFlex and CHARM3D, in time domain. Variations of Buoyancy to Weight Ratio (BWR), wave steepness/period, and water/submergence depth are considered as design and environmental parameters in the study. Two different mooring-line configurations, vertical and inclined, are studied to find an optimum design in terms of limiting tunnel motions and minimizing mooring-line tension. The numerical results are successfully validated by direct comparison against published experimental data. The results show that tunnel motions and tether tensions grow with wave height and period and decrease with submergence depth. The inclined mooring system is more effective in restricting tunnel motions compared to the vertical mooring system. Overall, the present study demonstrates the feasibility of this type of structure as an alternative to traditional bridges or under-seabed tunnels.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.10
• /
• pp.1317-1325
• /
• 2000

A new efficient numerical method for computing three-dimensional, unsteady, incompressible flows is presented. To eliminate the restriction of CFL condition, a fully-implicit time advancement in which the Crank-Nicolson method is used for both the diffusion and convection terms, is adopted. Based on an approximate block LU decomposition method, the velocity -pressure decoupling is achieved. The additional decoupling of the intermediate velocity components in the convection term is made for the fully -implicit time advancement scheme. Since the iterative procedures for the momentum equations are not required, the velocity components decouplings bring forth the reduction of computational cost. The second-order accuracy in time of the present numerical algorithm is ascertained by computing decaying vortices. The present decoupling method is applied to minimal channel flow unit with DNS (Direct Numerical Simulation).

• Bulletin of the Korean Mathematical Society
• /
• v.50 no.6
• /
• pp.2035-2051
• /
• 2013

In this paper, we study numerical schemes for solving multi-dimensional option pricing problem. We compare the direct solving method and the Operator Splitting Method(OSM) by using finite difference approximations. By varying parameters of the Black-Scholes equations for the maximum on the call option problem, we observed that there is no significant difference between the two methods on the convergence criterion except a huge difference in computation cost. Therefore, the two methods are compatible in practice and one can improve the time efficiency by combining the OSM with parallel computation technique. We show numerical examples including the Equity-Linked Security(ELS) pricing based on either two assets or three assets by using the OSM with the Monte-Carlo Simulation as the benchmark.