• 대한수학회지
• /
• 제50권1호
• /
• pp.81-93
• /
• 2013

We consider a distributed optimal flux control problem: finding the potential of which gradient approximates the target vector field under an elliptic constraint. Introducing the Lagrange multiplier and a change of variables the Euler-Lagrange equation turns into a coupled equation of an elliptic equation and a reaction diffusion equation. The change of variables reduces iteration steps dramatically when the Gauss-Seidel iteration is considered as a solution method. For the elliptic equation solver we consider the Cell Boundary Element (CBE) method, which is the finite element type flux preserving methods.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2001년도 춘계 학술대회논문집
• /
• pp.38-44
• /
• 2001

An efficient Gauss-Seidel time integration scheme is developed for solving the Euler and Navier-Stokes equations on unstructured meshes. Roe's FDS is used for the explicit residual computations and van Leer's FVS for evaluating implicit flux Jacobian. To reduce the memory requirement to a minimum level, off-diagonal flux Jacobian contributions are repeatedly calculated during the Gauss-Seidel sub-iteration process. Computational results based on the present scheme show that approximately $15\%$ of CPU time reduction is achieved while maintaining the memory requirement level to $50-60\%$ of the original Gauss-Seidel scheme.

• 대한인간공학회:학술대회논문집
• /
• 대한인간공학회 1997년도 추계학술대회논문집
• /
• pp.463-467
• /
• 1997

인간에게 운동감을 적절히 제시해주기 위하여는 Newton에 의한 운동의 세가지 법칙뿐만 아니라 EInstein의 상대성이론이 첨가되어야 한다. 즉, Newton운동의 제1법칙에 의하여 피실험자가 외력을 받지 않으면 등속운동 또는 정지상태를 계속 유지하게 되어 자신이 등속좌표계에 고정되어있기 때문에 시각적 인 정보가 없으면 어떠한 운동감도 못 느낀다. 이때 피실험자에게 정지해있는 기준좌표계에 대하여 등속 으로 움직이는 것을 인식시켜주기 위하여 피실험자에 대한 기준좌표계의 상대속도를 시각정보로 제공해 주어야 한다. 또한 Newton운동의 제2법칙에 의하여 똑같은 힘이 외력으로 작용하더라도 피실험자의 질량과 가속도는 서로 반비례하므로 화면이동속도변화를 피실험자의 질량에 반비례하도록 제시해 주어야 한다(김 정흠, 1982). 본 연구에서는 이러한 개념에 근거하여, 체중이 다른 여섯 피실험자들로 구성된 시스템에 대해서 각 피실험자에게 서로 다른 변위를 주고자할 때, 여섯가지 외력에 요구되는 작용시간을 Jacobi Iteration 방법과 Gauss-Seidel Iteration 방법으로 구하는 알고리즘을 제시하였다(D.V. Griffiths and I.M. Smith, 1991).

• Nuclear Engineering and Technology
• /
• 제54권2호
• /
• pp.532-545
• /
• 2022

In-depth convergence analyses for neutronics/thermal-hydraulics (T/H) coupled calculations are performed to investigate the performance of nonlinear methods based on the Fixed-Point Iteration (FPI). A simplified neutronics-T/H coupled system consisting of a single fuel pin is derived to provide a testbed. The xenon equilibrium model is considered to investigate its impact during the nonlinear iteration. A problem set is organized to have a thousand different fuel temperature coefficients (FTC) and moderator temperature coefficients (MTC). The problem set is solved by the Jacobi and Gauss-Seidel (G-S) type FPI. The relaxation scheme and the Anderson acceleration are applied to improve the convergence rate of FPI. The performances of solution schemes are evaluated by comparing the number of iterations and the error reduction behavior. From those numerical investigations, it is demonstrated that the number of FPIs is increased as the feedback is stronger regardless of its sign. In addition, the Jacobi type FPIs generally shows a slower convergence rate than the G-S type FPI. It also turns out that the xenon equilibrium model can cause numerical instability for certain conditions. Lastly, it is figured out that the Anderson acceleration can effectively improve the convergence behaviors of FPI, compared to the conventional relaxation scheme.

• 한국전자파학회논문지
• /
• 제25권5호
• /
• pp.601-606
• /
• 2014

IPO(Iterative Physical Optics) 방법은 대규모 물체의 산란파를 효과적으로 계산하는 고주파 근사 방법 중 하나인 PO(Physical Optics) 방법을 반복적으로 적용하는 계산방법이다. IPO 방법은 일차(first-order) PO 방법에서는 고려하지 못하는 다중 반사를 고려할 수 있어, 산란체 표면에 여기되는 전류의 정확도를 높일 수 있다. 그러므로 산란체의 RCS(Radar Cross Section)를 보다 정확하게 예측할 수 있다. 그러나 IPO 방법은 필요한 적분방정식을 정확하게 풀지 않아 수렴성에 문제가 생긴다. 그러므로 본 논문에서는 IPO 방법의 수렴성을 조절하기 위해, 행렬연산에 사용하는 Jacobi, Gauss-Seidel, SOR(Successive Over Relaxation) 그리고 Richardson 방법을 IPO 방법에 적용하였다. 그러므로 대규모 물체의 RCS 계산을 제안된 IPO 방법을 사용하여 효율적으로 계산할 수 있다. 또, 이들의 정확도를 시뮬레이션을 통해 검증하였다.

• 멤브레인
• /
• 제22권6호
• /
• pp.470-480
• /
• 2012

본 연구에서는 지름이 일정한 하나의 직선 관형 막을 가정하여 그 단위구간의 투과속도에 대한 물질수지를 세웠으며 이것을 기초로 하여 파이프의 정상상태 식과 함께 비선형 연립 2차 미분방정식을 이루었다. 이러한 관형 막의 압력손실을 표현한 연립방정식은 비선형이므로 Gauss-Seidel method와 같은 반복법에 의해서 해결될 수 있다. 이러한 수치해를 나타내기 위해 모사 알고리즘을 제시하였다. 또한 투과수는 운전 조건의 변화에 따라 변하므로 각 조건에서의 연립 방정식의 해를 수치적으로 적분하여 해결하였다. 모사의 결과를 검증, 해석하기 위해서 실제 중공사막과 유사한 관형 막을 사용하여 실험을 진행하였다. 본 연구에서는 중공사막을 유체가 흐를 때 발생하는 압력손실, 유량, 그리고 투과수의 관계를 분석적 방법을 이용해 제시하였고 이러한 이론적 기초를 바탕으로 실제 기공을 가진 막에 적용하여 그 정확성을 실험을 통하여 비교하였다.

• Journal of Biosystems Engineering
• /
• 제23권1호
• /
• pp.21-30
• /
• 1998

A theoretical model was developed to evaluate the effects of soil and airflow characteristics on the soil-air heat exchanger performances. The model, which includes three-dimensional transient energy and mass equilibrium-equation, was solved by using a computer program that uses Finite Difference Methods and Gauss-Seidel iteration computation. Energy gains, heat exchange efficiencies, and outlet air temperature are presented including the effects of soil moisture content, soil conductivity, soil thermal diffusivity, and soil initial temperature. Also, data related to the effects of airflow rate and inlet air temperature on the thermal performance of the system are presented. The results indicated that energy gains depend on soil conductivity, soil thermal diffusivity, and soil initial temperature. Heat exchange efficiencies relied on air mass flow rate and soil moisture content.

• KIEAE Journal
• /
• 제5권1호
• /
• pp.27-33
• /
• 2005

Global efforts have made to reduce energy consumption and $CO_2$ gas emission. One of the weakest parts for energy loss through the whole building components is building envelopes. Lots of technologies to increase the thermal performance of building envelopes have been introduced in recent year. Transparent Insulation Wall(TIW) is a new technology for building insulation and has been function both solar transmittance and thermal insulation. A mathematical model of a Transparent Insulation Wall equipped with south wall was proposed in order to predict thermal performance under varying climates(summer and winter). Unsteady state heat transfer equations were set up using an energy balance equation and solved using Gauss-Seidel iteration solution procedure. The thermal performance of the TIW determined from a wall surface and air layer temperature, non-airconditioned room temperature and air conditioning load. As a result, this numerical study shows that the TIW is effective in an air conditioning load reduction. Further experimental study is required to establish complete TIW system.

• Journal of Biosystems Engineering
• /
• 제22권4호
• /
• pp.459-468
• /
• 1997

To optimize the design and operation of a soil- air heat exchanger system, the effects of variables characterizing system design and operation on the performance of the system were analyzed by a theoretical model which included the three-dimensional transient heat conduction equation. The solution of the theoretical model was acquired by a computer program that uses Finite Difference Methods and Gauss-Seidel iteration computation, in which the time discretization scheme was an implicit difference appoximation. The computer program was validated first by comparison of the results for different grid sizes. Air outlet temperature, energy gain, and heat exchange efficiency of the system were analyzed based upon the tube diameter, tube length, tube thickness, and tube thermal diffusivity.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1995년도 가을 학술발표회 논문집
• /
• pp.266-273
• /
• 1995

An efficient dynamic analysis method is developed f3r building structures subjected to ground home loadings. The soil medium is modeled using the finite elements and infinite elements. Then, the dynamic stiffness of the soil medium is calculated at the interfacial nodes between the soil and the building foundation. The equivalent subway loading at the interfacial nodes are obtained from the wave propagation analysis of the subway loading through the soil medium. The dynamic response of the building Is computed using the mode superposition method equipped with gauss-seidel iteration technique. The analysis is carried out by the frequency domain and the time domain methods.