• Nuclear Engineering and Technology
• /
• v.49 no.5
• /
• pp.1019-1030
• /
• 2017

Diffusion is a crucial mechanism that regulates the migration of radioactive nuclides. In this study, an innovative numerical method was developed to simultaneously calculate the diffusion coefficient of both parent and, afterward, series daughter nuclides in a sequentially reactive through-diffusion model. Two constructed scenarios, a serial reaction (RN_1 ${\rightarrow}$ RN_2 ${\rightarrow}$ RN_3) and a parallel reaction (RN_1 ${\rightarrow}$ RN_2A + RN_2B), were proposed and calculated for verification. First, the accuracy of the proposed three-member reaction equations was validated using several default numerical experiments. Second, by applying the validated numerical experimental concentration variation data, the as-determined diffusion coefficient of the product nuclide was observed to be identical to the default data. The results demonstrate the validity of the proposed method. The significance of the proposed numerical method will be particularly powerful in determining the diffusion coefficients of systems with extremely thin specimens, long periods of diffusion time, and parent nuclides with fast decay constants.

• Journal of Ocean Engineering and Technology
• /
• v.34 no.6
• /
• pp.387-393
• /
• 2020

Frictional resistance comprises more than 60% of the total resistance for most merchant ships. Active and passive devices have been used to reduce frictional resistance, but the most effective and practical device is an air lubrication system. Such systems have been applied in several ships, and their effects have been verified in sea trials. On the other hand, there are some differences between the results predicted in model tests and those measured in sea trials. In this study, numerical analyses were carried out for a model and a full-scale ship. A new extrapolation method was proposed to improve the estimation of the full-scale resistance of a ship with an air lubrication system. The volume of fluid (VOF) method was considered for the numerical models of the air layer. The numerical method was validated by comparing the experimental data on the air layer pattern and the total resistance.

• Journal of the Korean Solar Energy Society
• /
• v.37 no.5
• /
• pp.77-84
• /
• 2017

The purpose of this paper is to propose the way of computing conduction time series factors (CTSF) using numerical method. After the accuracy of the numerical solution procedure being verified, the method is applied to the wall type 24 and roof type 14 of ASHARE to find the conduction time series coefficients, so called conduction time series factors. The results agree well with the values presented in the ASHRAE handbook. The method proposed can be easily applied to find unknown CTSF for more complex structures. It provides information about the temperature changes at a given location and time, thus validity of generated CTSF can be checked easily.

• Structural Engineering and Mechanics
• /
• v.68 no.6
• /
• pp.713-720
• /
• 2018

In this paper a new numerical method to determine the elastic curve of the simply supported beams of variable cross-section is demonstrated. In general case it needs to solve linear or small nonlinear second order differential equations with prescribed boundary conditions. For numerical solution the initial values of the slope and the deflection of the end cross-section of the beam is necessary. For obtaining the initial values a lively procedure is developed: it is a special application of the shooting method because boundary value problems can be transformed into initial value problems. As a result of these transformations the initial values of the differential equations are obtained with high accuracy. Procedure is applied for calculating of elastic curve of a simply supported beam of variable cross-section. Results of these numerical procedures, analytical solution of the linearized version and finite element method are compared. It is proved that the suggested procedure yields technically accurate results.

• Geomechanics and Engineering
• /
• v.25 no.4
• /
• pp.267-273
• /
• 2021

In this work, we compare the analytical solutions with the numerical solutions for photothermal interactions in semiconductor medium containing cylindrical cavity. This paper is devoted to a study of the photothermal interactions in semiconductor medium in the context of the coupled photo-thermal model. The basic equations are formulated in the domain of Laplace transform and the eigenvalue scheme are used to get the analytical solutions. The numerical solution is obtained by using the implicit finite difference method (IFDM). A comparison between the analytical solution and the numerical solutions are obtained. It is found that the implicit finite difference method (IFDM) is applicable, simple and efficient for such problems.

• Journal of Ocean Engineering and Technology
• /
• v.32 no.6
• /
• pp.447-457
• /
• 2018

In this study, a two-dimensional fully nonlinear transient wave numerical tank was developed using a desingularized indirect boundary integral equation method. The desingularized indirect boundary integral equation method is simpler and faster than the conventional boundary element method because special treatment is not required to compute the boundary integral. Numerical simulations were carried out in the time domain using the fourth order Runge-Kutta method. A mixed Eulerian-Lagrangian approach was adapted to reconstruct the free surface at each time step. A numerical damping zone was used to minimize the reflective wave in the downstream region. The interpolating method of a Gaussian radial basis function-type artificial neural network was used to calculate the gradient of the free surface elevation without element connectivity. The desingularized indirect boundary integral equation using an isolated point source and radial basis function has no need for information about the element connectivity and is a meshless method that is numerically more flexible. In order to validate the accuracy of the numerical wave tank based on the desingularized indirect boundary integral equation method and meshless technique, several numerical simulations were carried out. First, a comparison with numerical results according to the type of desingularized source was carried out and confirmed that continuous line sources can be replaced by simply isolated sources. In addition, a propagation simulation of a $2^{nd}$-order Stokes wave was carried out and compared with an analytical solution. Finally, simulations of propagating waves in shallow water and propagating waves over a submerged bar were also carried and compared with published data.

• International Journal of Automotive Technology
• /
• v.8 no.1
• /
• pp.33-38
• /
• 2007

Since it is difficult to analytically express the Melnikov function when a dynamic system possesses multiple saddle fixed points with homoclinic and/or heteroclinic orbits, this paper investigates a vehicle model with nonlinear suspension spring and hysteretic damping element, which exhibits multiple heteroclinic orbits in the unperturbed system. First, an algorithm for Melnikov integrals is developed based on the Melnikov method. And then the amplitude threshold of road excitation at the onset of chaos is determined. By numerical simulation, the existence of chaos in the present system is verified via time history curves, phase portrait plots and $Poincar{\acute{e}}$ maps. Finally, in order to further identify the chaotic motion of the nonlinear system, the maximal Lyapunov exponent is also adopted. The results indicate that the numerical method of estimating chaotic threshold is an effective one to complicated vehicle systems.

• Proceedings of the Korean Institute of Building Construction Conference
• /
• 2017.11a
• /
• pp.85-86
• /
• 2017

Fiber reinforced concrete as a construction material has been widely used. Fibers, as the reinforced component, the physical properties and the distribution influence the engineering properties of the composite. To illustrate the engineering properties, fiber distribution and orientation are necessary. Steel fibers can be easily captured by X-ray, but it is difficult them to express being numerical because they don't show as perfect circular shape on the grinding face. To get the more exact information for this, the numerical method for the orientation and distribution of fibers have to be more elaborately. This paper presents a possible method which makes the calculate for orientation possible.

• Bulletin of the Korean Mathematical Society
• /
• v.55 no.3
• /
• pp.881-898
• /
• 2018

This article deals with implementation of a high-order finite difference scheme for numerical solution of the incompressible Navier-Stokes equations on curvilinear grids. The numerical scheme is based on pseudo-compressibility approach. A fifth-order upwind compact scheme is used to approximate the inviscid fluxes while the discretization of metric and viscous terms is accomplished using sixth-order central compact scheme. An implicit Euler method is used for discretization of the pseudo-time derivative to obtain the steady-state solution. The resulting block tridiagonal matrix system is solved by approximate factorization based alternating direction implicit scheme (AF-ADI) which consists of an alternate sweep in each direction for every pseudo-time step. The convergence and efficiency of the method are evaluated by solving some 2D benchmark problems. Finally, computed results are compared with numerical results in the literature and a good agreement is observed.

• Journal of Power System Engineering
• /
• v.8 no.4
• /
• pp.49-56
• /
• 2004

To design very large ships, such as very large drillships, we have to estimate the hydroelastic responses of the very large ships in waves. A numerical procedure is described for estimating the hydroelastic responses of very large ships advancing with slow speed in waves. The developed numerical approach is based on a combination of the three-dimensional source distribution method and the finite element method, including fluid-structure interaction by regarding a very large ship as many hull elements connected with elastic beam elements. Numerical results are compared with experimental and numerical ones obtained in the literature. The results of comparison confirmed the validity of the proposed approach.