• Proceedings of the Korean Nuclear Society Conference
• /
• 1997.10a
• /
• pp.617-624
• /
• 1997

This Paper Presents an effective and simple procedure for the simulation of the motion of the solid-liquid interfacial boundary and the transient temperature field during phase change process. To accomplish this purpose, an iterative implicit solution algorithm has been developed by employing the dual reciprocity boundary element method. The dual reciprocity boundary element approach provided in this paper is much simpler than the usual boundary element method applying a reciprocity principle and an available technique for dealing with domain integral of boundary element formulation simultaneously. The effectiveness of the present analysis method have been illustrated through comparisons of the calculation results of an example with its semi-analytical or other numerical solutions where available.

• The Transactions of the Korean Institute of Electrical Engineers B
• /
• v.49 no.7
• /
• pp.457-461
• /
• 2000

The finite element method (FEM) is suitable for the analysis of a complicated region that includes nonlinear materials, whereas the boundary element method (BEM) is naturally effective for analyzing a very large region with linear characteristics. Therefore, considering the advantages in both methods, a novel algorithm for the alternate application of the FEM and BEM to magnetic field problems with the open boundary is presented. This approach avoids the disadvantages of the typical numerical methods with the open boundary problem such as a great number of unknown values for the FEM and non-symmetric matrix for the Hybrid FE-BE method. The solution of the overall problems is obtained by iterative calculations accompanied with the new acceleration method.

• Proceedings of the KIEE Conference
• /
• 1999.07a
• /
• pp.130-132
• /
• 1999

The finite element method (FEM) is suitable for the analysis of a complicated region that includes nonlinear materials, whereas the boundary element method (BEM) is naturally effective for analyzing a very large region with linear characteristics. Therefore, considering the advantages in both methods, a novel algorithm for the alternate application of the FEM and BEM to magnetic field problems with the open boundary is presented. This approach avoids the disadvantages of the typical numerical methods with the open boundary problem such as a great number of unknown values for the FEM and non-symmetric matrix for the Hybrid FE-BE method. The solution of the overall problems is obtained by iterative calculations accompanied with the new acceleration method.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2011.10a
• /
• pp.279-284
• /
• 2011

The calculation of Zwicker's loudness which is needed for multiple frequency response with a fine frequency resolution using the finite element (FE) procedure usually requires significant computation time since a numerical solution must be obtained for each considered frequency. Furthermore, if the analysis is the basis for an iterative optimization procedure this approach imposes high computational cost. In this work, we present an efficient approach for obtaining Zwicker's loudness via the Pad$\acute{e}$ approximants and applying in an acoustical topology optimization procedure. The paper is focused on an efficient and accurate calculation of Zwicker's loudness, design sensitivity analysis, and the acoustical topology optimization method by using Pad$\acute{e}$ approximants. The paper compares the efficient algorithm to results obtained by a standard FEM. Comparison are made both in terms of accuracy and in terms of CPU-times needed for the calculation.

• KSII Transactions on Internet and Information Systems (TIIS)
• /
• v.8 no.10
• /
• pp.3342-3360
• /
• 2014

In this paper, the capacity and BER performance of downlink distributed antenna systems (DAS) with transmit antenna selection and multiple receive antennas are investigated in MIMO composite channel, where path loss, Rayleigh fading and lognormal shadowing are all considered. Based on the performance analysis, using the probability density function (PDF) of the effective SNR and numerical integrations, tightly-approximate closed-form expressions of ergodic capacity and average BER of DAS are derived, respectively. These expressions have more accuracy than the existing expressions, and can match the simulation well. Besides, the outage capacity of DAS is also analyzed, and a tightly-approximate closed-form expression of outage capacity probability is derived. Moreover, a practical iterative algorithm based on Newton's method for finding the outage capacity is proposed. To avoid iterative calculation, another approximate closed-form outage capacity is also derived by utilizing the Gaussian distribution approximation. With these theoretical expressions, the downlink capacity and BER performance of DAS can be effectively evaluated. Simulation results show that the theoretical analysis is valid, and consistent with the corresponding simulation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.6
• /
• pp.753-765
• /
• 2003

A finite element code for the numerical solution of the Navier-Stokes equation is parallelized by vertex-oriented domain decomposition. To accelerate the convergence of iterative solvers like conjugate gradient method, parallel block ILU, iterative block ILU, and distributed ILU methods are tested as parallel preconditioners. The effectiveness of the algorithms has been investigated when P1P1 finite element discretization is used for the parallel solution of the Navier-Stokes equation. Two-dimensional and three-dimensional Laplace equations are calculated to estimate the speedup of the preconditioners. Calculation domain is partitioned by one- and multi-dimensional partitioning methods in structured grid and by METIS library in unstructured grid. For the domain-decomposed parallel computation of the Navier-Stokes equation, we have solved three-dimensional lid-driven cavity and natural convection problems in a cube as benchmark problems using a parallelized fractional 4-step finite element method. The speedup for each parallel preconditioning method is to be compared using upto 64 processors.

• Journal of the Society of Naval Architects of Korea
• /
• v.53 no.6
• /
• pp.447-455
• /
• 2016

Precise propulsion shafting alignment of ships is very important to prevent damage of its support bearings due to excessive reaction forces caused by hull deflection, forces acted on propeller and crankshaft, and so forth. In this paper, a new iterative shafting alignment calculation procedure considering the interaction between shaft deflection and oil film pressure of Sterntube Journal Bearing (SJB) bush with single or multiple slopes is proposed. The procedure is based on a pressure analysis to evaluate distributed equivalent support stiffness of SJB by solving Reynolds equation and a deflection analysis of shafting system by a finite element method based on Timoshenko beam theory. SJB is approximated with multi-point biaxial elastic supports equally distributed to its length. Their initial stiffness values are estimated from dynamic reaction force calculated by assuming SJB as single rigid support. Then, the shaft deflection and the support stiffness of SJB are sequentially and iteratively calculated by applying a criteria on deflection variation between sequential calculation results. To demonstrate validity and applicability of the proposed procedure for optimal slope design of SJB, numerical analysis results for a shafting system are described.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.4
• /
• pp.53-60
• /
• 1997

Nonlinear free surface flow phenomena have teen studied by several kinds of numerical methods, of which boundary element method has been known as most promising one. There, however, remain many difficulties to be solved in the numerical procedures by boundary element analysis. In this paper, an efficient calculation of elemental integrals and iterative solution algorithm for the resulting system of equations were thoroughly investigated in order to enhance the procedure of the boundary element analysis. Advantages of the herein developed boundary element analysis code are demonstrated in terms of accuracy and convergence for typical boundary-value problems with free surface.

• KSII Transactions on Internet and Information Systems (TIIS)
• /
• v.11 no.1
• /
• pp.253-271
• /
• 2017

In this paper, considering that perfect channel state information (CSI) is hard to achieve in practice, the downlink capacity of distributed antenna systems (DAS) with imperfect CSI and multiple receive antennas is investigated over composite Rayleigh fading channel. According to the performance analysis, using the numerical calculation, the probability density function (PDF) of the effective output SNR is derived. With this PDF, accurate closed-form expressions of ergodic capacity and outage probability of DAS with imperfect CSI are, respectively, obtained, and they include the ones under perfect CSI as special cases. Besides, the outage capacity of DAS in the presence of imperfect CSI is also derived, and a Newton's method based practical iterative algorithm is proposed to find the accurate outage capacity. By utilizing the Gaussian distribution approximation, another approximate closed-form expression of outage capacity is also derived, and it may simplify the calculation of accurate outage capacity. These theoretical expressions can provide good performance evaluation for downlink DAS for both perfect and imperfect CSI. Simulation results verify the effectiveness of the theoretical analysis, and the system capacity can be improved by increasing the receive antennas, and decreasing the estimation error or path loss. Moreover, the system can tolerate the estimation error variance up to about 0.01 with a slight degradation in the capacity.

• Computers and Concrete
• /
• v.11 no.5
• /
• pp.441-462
• /
• 2013

When using the conventional finite element method, a great number of grid nodes are necessary to describe the large and uneven temperature gradients in the concrete around cooling pipes when calculating the temperature field of mass concrete with cooling pipes. In this paper, the temperature gradient properties of the concrete around a pipe were studied. A new calculation method was developed based on these properties and an explicit iterative algorithm. With a small number of grid nodes, both the temperature distribution along the cooling pipe and the temperature field of the concrete around the water pipe can be correctly calculated with this new method. In conventional computing models, the cooling pipes are regarded as the third boundary condition when solving a model of concrete with plastic pipes, which is an approximate way. At the same time, the corresponding parameters have to be got by expensive experiments and inversion. But in the proposed method, the boundary condition is described strictly, and thus is more reliable and economical. And numerical examples were used to illustrate that this method is accurate, efficient and applicable to the actual engineering.