• Title/Summary/Keyword: SOR Method

Search Result 58, Processing Time 0.029 seconds

Errors of the conformal mapping method in the analysis of electrodes (전극 해석에서 등각사상법의 오차 평가)

  • 한영진;김창민
    • Journal of the Korean Institute of Telematics and Electronics A
    • /
    • v.32A no.8
    • /
    • pp.71-77
    • /
    • 1995
  • Evaluated are errors of the conformal mapping method in the analysis of a coplanar waveguides and asmmetric coplanar strips. The errors are mainly due to the negligence of fringe field effect after electrode structures are transformed to parallel electrodes. For the purpose of comparison, results of the conformal mapping(CM) method are compared with those of the successive over relaxation(SOR) method.

  • PDF

AN ASSESSMENT OF PARALLEL PRECONDITIONERS FOR THE INTERIOR SPARSE GENERALIZED EIGENVALUE PROBLEMS BY CG-TYPE METHODS ON AN IBM REGATTA MACHINE

  • Ma, Sang-Back;Jang, Ho-Jong
    • Journal of applied mathematics & informatics
    • /
    • v.25 no.1_2
    • /
    • pp.435-443
    • /
    • 2007
  • Computing the interior spectrum of large sparse generalized eigenvalue problems $Ax\;=\;{\lambda}Bx$, where A and b are large sparse and SPD(Symmetric Positive Definite), is often required in areas such as structural mechanics and quantum chemistry, to name a few. Recently, CG-type methods have been found useful and hence, very amenable to parallel computation for very large problems. Also, as in the case of linear systems proper choice of preconditioning is known to accelerate the rate of convergence. After the smallest eigenpair is found we use the orthogonal deflation technique to find the next m-1 eigenvalues, which is also suitable for parallelization. This offers advantages over Jacobi-Davidson methods with partial shifts, which requires re-computation of preconditioner matrx with new shifts. We consider as preconditioners Incomplete LU(ILU)(0) in two variants, ever-relaxation(SOR), and Point-symmetric SOR(SSOR). We set m to be 5. We conducted our experiments on matrices from discretizations of partial differential equations by finite difference method. The generated matrices has dimensions up to 4 million and total number of processors are 32. MPI(Message Passing Interface) library was used for interprocessor communications. Our results show that in general the Multi-Color ILU(0) gives the best performance.

New algorithm for simulating heat transfer in a complex CPFS (Cable Penetration Fire Stop)

  • Yun, Jong-Pil;Kwon, Seong-Pil;Cho, Jae-Kyu;Yoon, En-Sup
    • 제어로봇시스템학회:학술대회논문집
    • /
    • 2003.10a
    • /
    • pp.1798-1803
    • /
    • 2003
  • In this work the dynamic heat transfer occurring in a cable penetration fire stop system built in the firewall of nuclear power plants is three-dimensionally investigated to develop a test-simulator that can be used to verify effectiveness of the sealants. The dynamic heat transfer can be described by a partial differential equation (PDE) and its initial and boundary conditions. For the shake of simplicity PDE is divided into two parts; one corresponding to the heat transfer in the axial direction and the other corresponding to the heat transfer on the vertical layers. Two numerical methods, SOR (Sequential Over-Relaxation) and FEM (Finite Element Method), are implemented to solve these equations respectively. The axial line is discretized, and SOR is applied. Similarly, all the layers are separated into finite elements, where the time and spatial functions are assumed to be of orthogonal collocation state at each element. The heat fluxes on the layers are calculated by FEM. It is shown that the penetration cable influences the temperature distribution of the fire stop system very significantly. The simulation results are shown in the three-dimensional graphics for the understanding of the transient temperature distribution in the fire stop system.

  • PDF

Effect of filters and reconstruction method on Cu-64 PET image

  • Lee, Seonhwa;Kim, Jung min;Kim, Jung Young;Kim, Jin Su
    • Journal of Radiopharmaceuticals and Molecular Probes
    • /
    • v.3 no.2
    • /
    • pp.65-71
    • /
    • 2017
  • To assess the effects of filter and reconstruction of Cu-64 PET data on Siemens scanner, the various reconstruction algorithm with various filters were assessed in terms of spatial resolution, non-uniformity (NU), recovery coefficient (RC), and spillover ratio (SOR). Image reconstruction was performed using filtered backprojection (FBP), 2D ordered subset expectation maximization (OSEM), 3D reprojection algorithm (3DRP), and maximum a posteriori algorithms (MAP). For the FBP reconstruction, ramp, butterworth, hamming, hanning, or parzen filters were used. Attenuation or scatter correction were performed to assess the effect of attenuation and scatter correction. Regarding spatial resolution, highest achievable volumetric resolution was $3.08mm^3$ at the center of FOV when MAP (${\beta}=0.1$) reconstruction method was used. SOR was below 4% for FBP when ramp, Hamming, Hanning, or Shepp-logan filter were used. The lowest NU (highest uniform) after attenuation & scatter correction was 5.39% when FBP (parzen filter) was used. Regarding RC, 0.9 < RC < 1.1 was obtained when OSEM (iteration: 10) was used when attenuation and scatter correction were applied. In this study, image quality of Cu-64 on Siemens Inveon PET was investigated. This data will helpful for the quantification of Cu-64 PET data.

Efficient Iterative Physical Optics(IPO) Algorithms for Calculation of RCS (RCS 계산을 위한 효율적인 IPO 계산 방법)

  • Lee, Hyunsoo;Jung, Ki-Hwan;Chae, Dae-Young;Koh, Il-Suek
    • 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.

A Development of 3-D Numeric Model for the Confined Flow and Discharge under Sea Ground (해저 지중 피압유체의 흐름과 양수량 산출에 관한 3차원 수치모형의 개발)

  • Kim, Sang-Jun
    • Journal of Korean Society of Coastal and Ocean Engineers
    • /
    • v.22 no.6
    • /
    • pp.387-396
    • /
    • 2010
  • A 3-D numeric model for the confined transient flow under sea ground have been developed. This is FDM model using Gauss-Seidel SOR (successive over-relaxation). This model shows the similar head distribution pattern to Theis analytic solution and MODFLOW simulation. The input flowrate to the aquifer and discharge of well have been compared. And it have been found that mass balance is influenced by the weight factor ${\alpha}$, i.e. fullyimplicit method (${\alpha}$=1) shows 5% error, but when ${\alpha}$ becomes to 0.5(Crank and Nicolson method) the mass balance becomes worse and the model result diverges. And the convergency of the model is not much different when $\lambda$ (over-relaxation factor)=0.8~1.5, but when $\lambda$>1.5, the model result diverges. The test-run shows that the well discharge becomes smaller when another well is near. This model can cover the isotropy$(Kx{\neq}Ky{\neq}Kz)$ and inhomogeneity, and can be used for the selection of well site, discharge calculation, and head prediction in case of the artificial recharge etc.

A Numerical Study on a Circulation Control Foil using Coanda Effect (코앤다 효과를 이용한 순환 제어 날개의 수치적 연구)

  • J.J. Park;S.H. Lee
    • Journal of the Society of Naval Architects of Korea
    • /
    • v.37 no.2
    • /
    • pp.70-76
    • /
    • 2000
  • A numerical study on the viscous flow around a 2-dimensional circulation control foil is carried out for application on the field of naval architecture and ocean engineering. The governing equations are the RANS and the continuity equations. The equations are discretized by finite difference method and MAC method and the pressure poisson equation is calculate by a SOR method and an O-type non-staggered boundary fitted coordinate system which is overlapped near the slot is used to improve the numerical accuracy. Turbulence is approximated by a modified Baldwin-Lomax turbulence model. In the present paper, the Coanda effect on a 2-dimensional foil of a 20% thickness ellipse with modified rounded trailing edge has been numerically studied. The change in drag and lift of the foil with various jet momentums are calculated and compared to the experimental results to show good agreements.

  • PDF

Mach-Zehnder Type $Ti:LiNbO_3$ Traveling-Wave Optical Modulator with m-Section Phase Reversal (m-section의 위상반전이 있는 Mach-Zehnder형 진행파 $Ti:LiNbO_3$ 광변조기)

  • Lee, Woo-Jin;Kim, Kyoung-Am;Kim, Woo-Kyung;Kim, Chang-Min
    • Journal of the Institute of Electronics Engineers of Korea SD
    • /
    • v.39 no.7
    • /
    • pp.26-36
    • /
    • 2002
  • Mach-Zehnder type traveling-wave optical modulators with 3-section and 5-section phase reversal were designed and fabricated on z-cut $LiNbO_3$ substrates. Optical waveguides were designed by means of the FDM(Finite Difference Method). Design of CPW traveling-wave electrodes were performed by the SOR(Successive Over Relaxation) in the active region and by the CMM(Conformal Mapping) in the input/output section. The optical response R(${\omega}$ calculated based on the measured S-parameters showed the bandwidth of 15GHz centered at 25GHz for the 3-section, and the bandwidth of 22GHz at 45GHz for the 5-section.

Numerical Solution of Nonlinear Diffusion in One Dimensional Porous Medium Using Hybrid SOR Method

  • Jackel Vui Lung, Chew;Elayaraja, Aruchunan;Andang, Sunarto;Jumat, Sulaiman
    • Kyungpook Mathematical Journal
    • /
    • v.62 no.4
    • /
    • pp.699-713
    • /
    • 2022
  • This paper proposes a hybrid successive over-relaxation iterative method for the numerical solution of a nonlinear diffusion in a one-dimensional porous medium. The considered mathematical model is discretized using a computational complexity reduction scheme called half-sweep finite differences. The local truncation error and the analysis of the stability of the scheme are discussed. The proposed iterative method, which uses explicit group technique and modified successive over-relaxation, is formulated systematically. This method improves the efficiency of obtaining the solution in terms of total iterations and program elapsed time. The accuracy of the proposed method, which is measured using the magnitude of absolute errors, is promising. Numerical convergence tests of the proposed method are also provided. Some numerical experiments are delivered using initial-boundary value problems to show the superiority of the proposed method against some existing numerical methods.

comparison of Numercal Methods for Obtaining 2-D Impurity Profile in Semiconductor (반도체 내에서의 2차원 불순물 분포를 얻기 위한 수치해법의 비교)

  • Yang, Yeong-Il;Gyeong, Jong-Min;O, Hyeong-Cheol
    • Journal of the Korean Institute of Telematics and Electronics
    • /
    • v.22 no.3
    • /
    • pp.95-102
    • /
    • 1985
  • An efficient numerical scheme for assessing the two-dimensional diffusion problem for modelling impurity profile in semiconductor is described. 4 unique combination of ADI (Al-ternating Direction Bmplicit) method and Gauss Elimination has resulted in a reduction of CPU time for most diffusion processes by a factor of 3, compared to other iteration schemes such as SOR (Successive Over-Relaxation) or Stone's iterative method without additional storage re-quirement. Various numerical schemes were compared for 2-D as well as 1-0 diffusion profile in terms of their CPU time while retaining the magnitude of relative error within 0.001%. good agree-ment between 1-D and 2-D simulation profile as well as between 1-D simulation profile and experiment has been obtained.

  • PDF