• Nuclear Engineering and Technology
• /
• v.47 no.5
• /
• pp.546-558
• /
• 2015

Chord length sampling method in Monte Carlo simulations is a method used to model spherical particles with random sampling technique in a stochastic media. It has received attention due to the high calculation efficiency as well as user convenience; however, a technical issue regarding boundary effect has been noted. In this study, after analyzing the distribution characteristics of spherical particles using an explicit method, an alternative chord length sampling method is proposed. In addition, for modeling in finite media, a correction method of the boundary effect is proposed. Using the proposed method, sample probability distributions and relative errors were estimated and compared with those calculated by the explicit method. The results show that the reconstruction ability and modeling accuracy of the particle probability distribution with the proposed method were considerably high. Also, from the local packing fraction results, the proposed method can successfully solve the boundary effect problem. It is expected that the proposed method can contribute to the increasing of the modeling accuracy in stochastic media.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.2
• /
• pp.93-114
• /
• 2019

In this article, we introduce a finite difference method for solving the Navier-Stokes equations in rectangular domains. The method is proved to be energy stable and shown to be second-order accurate in several benchmark problems. Due to the guaranteed stability and the second order accuracy, the method can be a reliable tool in real-time simulations and physics-based animations with very dynamic fluid motion. We first discuss a simple convection equation, on which many standard explicit methods fail to be energy stable. Our method is an implicit Runge-Kutta method that preserves the energy for inviscid fluid and does not increase the energy for viscous fluid. Integration-by-parts in space is essential to achieve the energy stability, and we could achieve the integration-by-parts in discrete level by using the Marker-And-Cell configuration and central finite differences. The method, which is implicit and second-order accurate, extends our previous method [1] that was explicit and first-order accurate. It satisfies the energy stability and assumes rectangular domains. We acknowledge that the assumption on domains is restrictive, but the method is one of the few methods that are fully stable and second-order accurate.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1993.10a
• /
• pp.248-255
• /
• 1993

Recently, due to the increasing attention to the structural safety under uncertain environments, many researches on the structural reliability analysis have been peformed. Some useful methods are available to evaluate performance reliability of structures with explicit limit states. However, for large structures, in which structural behaviors can be analyzed with finite element models and the limit states are only expressed implicitly, Monte-Carlo simulation method has been mainly used. However, Monte-Carlo simulation method spends too much computational time on repetitive structural analysis. Many alternative methods are suggested to reduce the computational work required in Monte-Carlo simulation. Response surface method is widely used to improve the efficiency of structural reliability analysis. Response surface method is based on the concept of approximating simple polynomial function of basic random variables for the limit state which is not easily expressed in explicit forms of design random variables. The response surface method has simple algorithm. However, the accuracy of results highly depends on how properly the stochastic characteristics of the original limit state has been represented by approximated function, In this study, an improved response surface method is proposed in which the sampling points for creating response surface are modified to represent the failure surface more adequately and the combined use of a linear response surface function and Rackwitz-Fiessler method has been employed. The method is found to be more effective and efficient than previous response surface methods. In addition more consistent convergence is achieved, Accuracy of the proposed method has been investigated through example.

• Nuclear Engineering and Technology
• /
• v.43 no.3
• /
• pp.271-278
• /
• 2011

When performing fire PSA in a nuclear power plant, an event mapping method, using an internal event PSA model, is widely used to reduce the resources used by fire PSA model development. Feasible initiating events and component failure events due to fire are identified to transform the fault tree (FT) for an internal event PSA into one for a fire PSA using the event mapping method. A surrogate event or damage term method is used to condition the FT of the internal PSA. The surrogate event or the damage term plays the role of flagging whether the system/component in a fire compartment is damaged or not, depending on the fire being initiated from a specified compartment. These methods usually require explicit states of all compartments to be modeled in a fire area. Fire event scenarios, when using explicit identification, such as surrogate or damage terms, have two problems: (1) there is no consideration of multiple fire propagation beyond a single propagation to an adjacent compartment, and (2) there is no consideration of simultaneous fire propagations in which an initiating fire event is propagated to multiple paths simultaneously. The present paper suggests a fire propagation equation to identify all possible fire event scenarios for an explicitly treated fire event scenario in the fire PSA. Also, a method for separating fire events was developed to make all fire events a set of mutually exclusive events, which can facilitate arithmetic summation in fire risk quantification. A simple example is given to confirm the applicability of the present method for a $2{\times}3$ rectangular fire area. Also, a feasible asymptotic approach is discussed to reduce the computational burden for fire risk quantification.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.12
• /
• pp.1153-1161
• /
• 2010

An efficient solution method for harmonic balance techniques with Fourier transform is presented for periodic unsteady flow problems. The present partially-implicit harmonic balance treats the flux terms implicitly and the harmonic source term is solved explicitly. The convergence of the partially Implicit method is much faster than the explicit Runge-Kutta harmonic balance method. The method does not need to compute the additional flux Jacobian matrix from the implicit harmonic source term. Compared with fully implicit harmonic balance method, this partial approach turns out to have good convergence property. Oscillating flows over NACA0012 airfoil are considered to verify the method and to compare with results of explicit R-K(Runge-Kutta) and dual time stepping methods.

• International Journal of CAD/CAM
• /
• v.6 no.1
• /
• pp.139-148
• /
• 2006

Catmull-Clark subdivision scheme provides a powerful method for building smooth and complex surfaces. But the number of faces in the uniformly refined meshes increases exponentially with respect to subdivision depth. Adaptive tessellation reduces the number of faces needed to yield a smooth approximation to the limit surface and, consequently, makes the rendering process more efficient. In this paper, we present a new adaptive tessellation method for general Catmull-Clark subdivision surfaces. Different from previous control mesh refinement based approaches, which generate approximate meshes that usually do not interpolate the limit surface, the new method is based on direct evaluation of the limit surface to generate an inscribed polyhedron of the limit surface. With explicit evaluation of general Catmull-Clark subdivision surfaces becoming available, the new adaptive tessellation method can precisely measure error for every point of the limit surface. Hence, it has complete control of the accuracy of the tessellation result. Cracks are avoided by using a recursive color marking process to ensure that adjacent patches or subpatches use the same limit surface points in the construction of the shared boundary. The new method performs limit surface evaluation only at points that are needed for the final rendering process. Therefore it is very fast and memory efficient. The new method is presented for the general Catmull-Clark subdivision scheme. But it can be used for any subdivision scheme that has an explicit evaluation method for its limit surface.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.2
• /
• pp.87-95
• /
• 2018

This paper introduces a new explicit spectral element method for the simulation of SH-waves in semi-infinite domains. To simulate the wave motion in unbounded domains, it is necessary to reduce the infinite extent to a finite computational domain of interest. To prevent the wave reflection from the trunctated boundaries, perfectly matched layer(PML) wave-absorbing boundary is introduced. The forward problem for simulating SH-waves in PML-truncated domains can be formulated as second-order PDEs. The second-order semi-discrete form of the governing PDEs is constructed by using a mixed spectral elements with Legendre-gauss-Lobatto quadrature method, which results in a diagonalized mass matrix. Then the second-order semi-discrete form is transformed to a first-order, whose solutions are calculated by the fourth-order Runge-Kutta method. Numerical examples showed that solutions of SH-wave in the two-dimensional analysis domain resulted in stable and accurate, and reflections from truncated boundaries could be reduced by using PML boundaries. Elastic wave propagation analysis using explicit time integration method may be apt for solving larger domain problems such as three-dimensional elastic wave problem more efficiently.

• Journal of the Korean Society for Library and Information Science
• /
• v.52 no.4
• /
• pp.91-110
• /
• 2018

We proposed a video summarization algorithm based on an ERP (Event Related Potentials)-based topic relevance model, a MMR (Maximal Marginal Relevance), and discriminant analysis to generate a semantically meaningful video skim. We then conducted implicit and explicit evaluations to evaluate our proposed ERP/MMR-based method. The results showed that in the implicit and explicit evaluations, the average scores of the ERP / MMR methods were statistically higher than the average score of the SBD (Shot Boundary Detection) method used as a competitive baseline, respectively. However, there was no statistically significant difference between the average score of ERP/MMR (${\lambda}=0.6$) method and that of ERP/MMR (${\lambda}=1.0$) method in both assessments.

• Structural Engineering and Mechanics
• /
• v.38 no.2
• /
• pp.211-229
• /
• 2011

In this paper, we present a new explicit procedure using periodic cubic B-spline interpolation polynomials to solve linear and nonlinear dynamic equation of motion governing single degree of freedom (SDOF) systems. In the proposed approach, a straightforward formulation was derived from the approximation of displacement with B-spline basis in a fluent manner. In this way, there is no need to use a special pre-starting procedure to commence solving the problem. Actually, this method lies in the case of conditionally stable methods. A simple step-by-step algorithm is implemented and presented to calculate dynamic response of SDOF systems. The validity and effectiveness of the proposed method is demonstrated with four examples. The results were compared with those from the numerical methods such as Duhamel integration, Linear Acceleration and also Exact method. The comparison shows that the proposed method is a fast and simple procedure with trivial computational effort and acceptable accuracy exactly like the Linear Acceleration method. But its power point is that its time consumption is notably less than the Linear Acceleration method especially in the nonlinear analysis.

• Atmosphere
• /
• v.23 no.1
• /
• pp.33-37
• /
• 2013

In meteorology, exploitation of Kalman filter as a data assimilation system is virtually impossible due to simultaneous requirements of adjoint model and large computer resource. The other substitute of utilizing ensemble Kalman filter is only affordable by compensating an enormous usage of computing resource. Furthermore, the latter employs ensemble integration sets for evolving the background error covariance matrix by compensating the dynamical feature of the temporal evolution of weather conditions. We propose a new implementation method that works without the adjoint model by utilizing the explicit expression of the background error covariance matrix in backward evolution. It will also break a barrier in the evolution of the covariance matrix. The method may be applied with a slight modification to the real time assimilation or the retrospective analysis.