• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.18 no.9
• /
• pp.1023-1029
• /
• 2007

In this paper, we propose the method to solve the BLT equations using Laguerre polynomials in time domain. The solution of BLT equations is obtained by recursive, differential and integral properties of Laguerre polynomials. The verification of the proposed method is tested by applying it to the two-wired transmission line with resistors and capacitors, which is illuminated by the electromagnetic plane wave pulse. And the result is compared with the corresponding transient responses obtained from inverse fast Fourier transform(IFFT) of the frequency domain solutions of BLT equations.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.23 no.2
• /
• pp.169-176
• /
• 2012

This paper uses loop-star basis functions to overcome the low frequency breakdown problem in method of moments (MoM) based on electric field integral equation(EFIE). In addition, p-Type Multiplicative Schwarz preconditioner (p-MUS) technique is employed to reduce the number of iterations required for the conjugate gradient method(CGM). Low frequency instability with Rao Wilton Glisson(RWG) basis functions in EFIE can be resolved using loop-start basis functions and frequency normalized techniques. However, loop-star basis functions, consisting of irrotational and solenoidal components of RWG basis functions, require a large number of iterations to calculate a solution through iterative methods, such as conjugate gradient method(CGM), due to high condition number. To circumvent this problem, in this paper, the pMUS preconditioner technique is proposed to reduce the number of iterations in CGM. Simulation results show that pMUS preconditioner is much faster than block diagonal preconditioner(BDP) when the sparsity of pMUS is the same as that of BDP.

• Nuclear Engineering and Technology
• /
• v.54 no.6
• /
• pp.2023-2030
• /
• 2022

A novel pulse fitting analysis (PFA) method is presented for the acquisition of nuclear spectra. The charging process of the feedback capacitor in the resistive feedback charge-sensitive preamplifier is equivalent to the impulsive pulse, and its impulse response function (IRF) can be obtained by non-linear fitting of the falling edge of the nuclear pulse. The integral of the IRF excluding the baseline represents the energy deposition of the particles in the detector. In addition, since the non-linear fitting process in PFA method is difficult to achieve in the conventional architecture of spectroscopy system, a new multi-channel analyzer (MCA) based on Zynq SoC is proposed, which transmits all the data of nuclear pulses from the programmable logic (PL) to the processing system (PS) by high-speed AXI-Stream in order to implement PFA method with precision. The linearity of new MCA has been tested. The spectrum of ¹³⁷Cs was obtained using LaBr₃(Ce) scintillator detector, and was compared with commercial MCA by ORTEC. The results of tests indicate that the MCA based on PFA method has the same performance as the commercial MCA based on pulse height analysis (PHA) method and excellent linearity for γ-rays with different energies, which infers that PFA method is an effective and promising method for the acquisition of spectra. Furthermore, it provides a new solution for nuclear pulse processing algorithms involving regression and iterative processes.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.2
• /
• pp.207-216
• /
• 2003

In this paper, the finite element method for the elastodynamic axisymmetric fracture analysis is presented in matrix form through the application of the Galerkin method to the time integral equations of motion with no inertia forces. Isoparametric quadratic quadrilateral element and triangular crack tip singular elements with one-quarter node are used in the mesh division of the finite element model. To show the validity and accuracy of the proposed method, the infinite elastic medium with the penny shaped crack is solved first and compared with the analytical solution and the numerical results by the finite difference method and the boundary element method existing in the published literatures, and then the dynamic stress intensity factors of solid and hollow cylinders of finite dimensions haying penny-shaped cracks and internal and external circumferential tracks are computed in detail.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.7 no.2
• /
• pp.156-162
• /
• 1995

A numerical method for a 2-D oil boom model considering the flexibility of skirt has been developed The neater is assumed rigid and the skirt is tensioned membrane having a point mass at its end The fluid motion is potential. The kinematic condition which demands the continuity of the displacement is imposed at the joint between the floater and the skirt. The dynamic condition for the point mass is imposed at the bottom end of the skirt. The numerical method is based on the Green's function method in the frame of linear potential theory. It finds it's solution simultaneously from the total system of three equations, integral equation, the equation of motion of the floater and the equilibrium equation of the deformation of the skirt. Integral equation is derived by applying the Green's theorem to radiation potential and Green's function. Proper descretization of those three equations leads to the system of a linear algebraic equation. Due to the flexibility of skirt the motion of floater can be diminished in some range of wave frequency and furthermore the mechanism of resonance of the oil boom can be changed. The motion responses of various oil booms have been compared varying the length of the skirt and the point mass. The numerical method has been validated indirectly from the good correspondence between the motion responses of the flexible skirt model and the rigid skirt model in low frequency limit.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.7 no.1
• /
• pp.11-22
• /
• 1987

A plane buoyant jet discharged vertically upward into a crossflow is analyzed by numerical solution of the governing equations of continuity, momentum and constituent transport. The turbulent transport is modelled by the Prandtl's mixing length theory. In the numerical solution procedure, the governing equations are transformed by stream function and vorticity transport, non-dimensionalyzed by discharge velocity, slot width, and parameters representing flow characteristics, and solved by Gauss-Seidel iteration method with successive underrelaxation. The numerical experiments were performed for the region of established flow of buoyant jet in the range of discharge densimetric Froude number of 4 to 32 and in the range of velocity ratio of 8 to 15, which is the ratio of discharge velocity to crossflow velocity. Variations of velocities and temperatures, flow patterns and vorticity patterns of receiving water due to buoyant jet were investigated. Also investigated are the effects of velocity ratio and discharge densimetric Froude number on the trajectories of buoyant jet. Computed are velocities, temperatures and local densimetric Froude numbers along the trajectory of the buoyant jet. Spreading rate and dispersion ratio were analyzed in terms of discharge densimetric Froude number, local densimetric Froude number and distance from the source along the jet trajectory. It was noted that the similarity law holds in both the profiles of velocity and temperatures across the jet trajectory and the integral type analysis of Gaussian distribution is applicable.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.3
• /
• pp.25-34
• /
• 1993

The long-term behavior, such as in excavation problems of weak medium, can be dealt with by the elasto-viscoplasticity models. In this paper, a combined formulation of elasto-viscoplasticity using boundary elements and finite elements without using internal cells is presented. The domain integral introduced due to the viscoplastic stresses is transformed into a boundary integral applying direct integration in cylindrical coordinates. The results of the developed boundary element analysis are compared with those from the explicit solution and from the finite element analysis. It is observed that the boundary element analysis without internal cells results in some error because of its deficiency in handling the nonlinearity in local stress concentration. Therefore, a coupled analysis of boundary elements and finite elements, in which finite elements are used in the area of stress concentration, is developed. The coupled method is applied to a time dependent inelastic problem with semi-infinite boundaries. It results in reasonable solution compared with other methods where relatively higher degree of freedoms are employed. Thus, it is concluded that the combined analysis may be used for such problems in the effective manner.

• Geophysics and Geophysical Exploration
• /
• v.4 no.2
• /
• pp.25-33
• /
• 2001

Loop-loop electromagnetic (EM) survey in frequency domain has been carried out in order to provide basic solution to geotechnical applications. Source and receiver configuration may be horizontal co-planar (HCP) and/or vertical co-planar (VCP). Three quadrature components of mutual impedance ratio for each configuration are used to construct the subsurface image. For the purpose of obtaining the model response and validating the reasonable performance of the inversion, we obtained each responses of two-layered and three-layered earth models and two-dimensional (2-D) isolated anomalous body. The response of 2-D isolated anomalous body has been calculated using extended Born approximation for the solution of 2.5-D integral equation describing EM scattering problem. As a result of the least-squares inversion with variable Lagrangian multiplier, we could construct more resolvable image from HCP data than VCP data. Furthermore, joint inversion of HCP and VCP data made better stability and resolution of the inversion. Resistivity values, however, did not exactly match the true ones. Loop-loop EM field data was obtained with EM34-3XL system manufactured by Geonics Ltd. (Canada). Electrical resistivity survey was conducted on the same line for the comparison in advance. Since the constructed image from loop-loop EM data by 2-D inversion algorithm showed almost similar resistivity distribution to that from electrical resistivity one, we expect the developed 2.5-D loop-loop EM inversion program can be applied for the reconnaissance site survey.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.25 no.5
• /
• pp.601-606
• /
• 2014

The IPO(Iterative Physical Optics) method repeatedly applies the well-known PO(Physical Optics) approximation to calculate the scattered field by a large object. Thus, the IPO method can consider the multiple scattering in the object, which is ignored for the PO approximation. This kind of iteration can improve the final accuracy of the induced current on the scatterer, which can result in the enhancement of the accuracy of the RCS(Radar Cross Section) of the scatterer. Since the IPO method can not exactly but approximately solve the required integral equation, however, the convergence of the IPO solution can not be guaranteed. Hence, we apply the famous techniques used in the inversion of a matrix to the IPO method, which include Jacobi, Gauss-Seidel, SOR(Successive Over Relaxation) and Richardson methods. The proposed IPO methods can efficiently calculate the RCS of a large scatterer, and are numerically verified.

• Journal of Power Electronics
• /
• v.8 no.1
• /
• pp.10-22
• /
• 2008

The wind-driven doubly fed induction generator (DFIG) is currently under pressure to be more grid-compatible. The main concern is the fault ride-through (FRT) requirement to keep the generator connected to the grid during faults. In response to this, the paper introduces a novel model and new control scheme for the DFIG. The model provides a means of direct stator power control and considers the stator transients. On the basis of the derived model, a robust linear quadratic (LQ) controller is synthesized. The control law has proportional and integral actions and takes account of one sample delay in the input owing to the microprocessor's execution time. Further, the influence of the grid voltage imperfection is mitigated using frequency shaped cost functional method. Compensation of the rotor current pulsations is proposed to improve the FRT capability as well as the generator performance under grid voltage unbalance. As a consequence, the control system can achieve i) fast direct power control without instability risk, ii) alleviation of the problems associated with the DFIG operation under unbalanced grid voltage, and iii) high probability of successful grid FRT. The effectiveness of the proposed solution is confirmed through simulation studies on 2MW DFIG.