• 제목/요약/키워드: ODE solver

검색결과 8건 처리시간 0.019초

Conservative Upwind Correction Method for Scalar Linear Hyperbolic Equations

  • Kim, Sang Dong;Lee, Yong Hun;Shin, Byeong Chun
    • Kyungpook Mathematical Journal
    • /
    • 제61권2호
    • /
    • pp.309-322
    • /
    • 2021
  • A conservative scheme for solving scalar hyperbolic equations is presented using a quadrature rule and an ODE solver. This numerical scheme consists of an upwind part, plus a correction part which is derived by introducing a new variable for the given hyperbolic equation. Furthermore, the stability and accuracy of the derived algorithm is shown with numerous computations.

인터벌 연산 기반의 곡면간 교차선 계산에 관한 연구 (A Study on Intersection Computation in Interval Arithmetic)

  • 고광희
    • 한국CDE학회논문집
    • /
    • 제15권3호
    • /
    • pp.178-188
    • /
    • 2010
  • This paper addresses the problem of determining if two surfaces intersect tangentially or transversally in a mathematically consistent manner and approximating an intersection curve. When floating point arithmetic is used in the computation, due to the limited precision, it often happens that the decision for tangential and transversal intersection is not clear cut. To handle this problem, in this paper, interval arithmetic is proposed to use, which provides a mathematically consistent way for such decision. After the decision, the intersection is traced using the validated ODE solver, which runs in interval arithmetic. Then an iterative method is used for computing the accurate intersection point at a given arc-length of the intersection curve. The computed intersection points are then approximated by using a B-spline curve, which is provided as one instance of intersection curve for further geometric processing. Examples are provided to demonstrate the proposed method.

토폴로지와 수치적 정확도를 통합한 기하모델링에 관한 연구: 곡면간 교차선 (A Study on Unifying Topology and Numerical Accuracy in Geometric Modeling: Surface to Surface Intersections)

  • 고광희
    • 한국CDE학회논문집
    • /
    • 제12권5호
    • /
    • pp.344-353
    • /
    • 2007
  • In this paper, we address the problem of robust geometric modeling with emphasis on surface to surface intersections. We consider the topology and the numerical accuracy of an intersection curve to find the best approximation to the exact one. First, we perform the topological configuration of intersection curves, from which we determine the starting and ending points of each monotonic intersection curve segment along with its topological structure. Next, we trace each monotonic intersection curve segment using a validated ODE solver, which provides the error bounds containing the topological structure of the intersection curve and enclosing the exact root without a numerical instance. Then, we choose one approximation curve and adjust it within the bounds by minimizing an objective function measuring the errors from the exact one. Using this process, we can obtain an approximate intersection curve which considers the topology and the numerical accuracy for robust geometric modeling.

강성 구속기계계의 해법에 관한 연구 (A Study on Numerical Solution Methods in Stiff Constrained Mechanical Systems)

  • 한형석;박태원
    • 대한기계학회논문집
    • /
    • 제19권2호
    • /
    • pp.374-381
    • /
    • 1995
  • In this paper, integration methods for stiff constrained mechanical systems are studied to analyze dynamic response of the stiff mechanical system efficiently. The stiff, non-stiff systems are identified by using eigenvalues of the Jacobian of the systems. To integrate both stiff system and non-stiff system efficiently a new switching method between the non-stiff differential equation solver and the stiff differential equation solver is presented.

THE BUCHSTAB'S FUNCTION AND THE OPERATIONAL TAU METHOD

  • Aliabadi, M.Hosseini
    • Journal of applied mathematics & informatics
    • /
    • 제7권3호
    • /
    • pp.905-915
    • /
    • 2000
  • In this article we discuss some aspects of operational Tau Method on delay differential equations and then we apply this method on the differential delay equation defined by $\omega(u)\;=\frac{1}{u}\;for\;1\lequ\leq2$ and $(u\omega(u))'\;=\omega(u-1)\;foru\geq2$, which was introduced by Buchstab. As Khajah et al.[1] applied the Recursive Tau Method on this problem, they had to apply that Method under the Mathematica software to get reasonable accuracy. We present very good results obtained just by applying the Operational Tau Method using a Fortran code. The results show that we can obtain as much accuracy as is allowed by the Fortran compiler and the machine-limitations. The easy applications and reported results concerning the Operational Tau are again confirming the numerical capabilities of this Method to handle problems in different applications.

계면마찰항을 고려한 이상유동에서 파동전파에 대한 수치적 연구 (TWO-PHASE WAVE PROPAGATIONS PREDICTED BY HLL SCHEME WITH INTERFACIAL FRICTION TERMS)

  • 염금수;장근식;정문선
    • 한국전산유체공학회:학술대회논문집
    • /
    • 한국전산유체공학회 2009년 추계학술대회논문집
    • /
    • pp.115-119
    • /
    • 2009
  • We numerically investigated propagation of various waves in the two-phase flows such as sound wave, shock wave, rarefaction wave, and contact discontinuity in terms of pressure, void fraction, velocity and density of the two phases. The waves have been generated by a hydrodynamic shock tube, a pair of symmetric impulsive expansion, impulsive pressure and impulsive void waves. The six compressible two-fluid two-phase conservation laws with interfacial friction terms have been solved in two fractional steps. The first PDE Operator is solved by the HLL scheme and the second Source Operator by the semi-implicit stiff ODE solver. In the HLL scheme, the fastest wave speeds were estimated by the analytic eigenvalues of an approximate Jacobian matrix. We have discussed how the interfacial friction terms affect the wave structures in the numerical solution.

  • PDF

CPFS 내에서 일어나는 동적 열전도 현상을 해석하기 위한 수식 및 혼합알고리즘 개발 (Development of the model and the hybrid algorithm toy analyzing the dynamic heat conduction in the CPES system)

  • 윤종필;권성필;윤인섭
    • 한국가스학회:학술대회논문집
    • /
    • 한국가스학회 2003년도 추계학술발표회 논문집
    • /
    • pp.120-125
    • /
    • 2003
  • 본 연구는 원자력 발전소에 있는 방화벽의 케이블 관통부위에 설치된 CPFS(Cable Penetration Fire Stop)시스템 안에서 일어나는 동적열전달 현상을 3 차원으로 나타낼 수 있는 시험시뮬레이터에 사용될 수학적 모델과 수치계산 알고리즘의 개발에 관한 것이다. CPFS 내에서 일어나는 열전도 현상을 나타내는 지배방정식은 주어진 조건들 하에서 포물선형 편미분방정식(Parabolic PDE)으로 나타난다. 문제를 단순화하기 위해 열의 흐름을 두 성분으로 나누었다 즉, 케이블과 평행한 선을 따라서 일어나는 열전도와 벽면과 평행한 평면 위에서 일어나는 열전도로 나누었다. 먼저 선을 따라 일어나는 동적 열전도 현상을 나타내는 PDE를 연속과완화(SOR: Successive Over-Relaxation)를 적용하여 유한한 불연속점들에 대한 연립 상미분방정식(ODE)으로 전환했고, 그 연립방정식은 ODE Solver 를 이용하여 풀 수 있었다. 둘째로, 각 불연속 점에 위치한 평면 위에서 일어나는 열전도를 계산하기 위해서, 유한요소의 합을 근사식으로 이용하여 PDE를 ODE로 전환해서 계산하는 유한요소법(Finite Element Method)이 이용된다. 여기서 시간과 공간의 함수 T(x, y, z, t)인 온도는 각 선의 점들과 각 평면의 요소들에 대해서 일정한 시간간격으로 초기온도와 경계온도를 업데이트하여 계산을 반복한다. 이러한 일련의 계산결과를 바탕으로 CPFS 시스템 내에서의 온도분포의 동적인 변화를 해석한다. 결론적으로 관통하는 케이블이 CPFS 시스템의 온도분포에 매우 중요한 역할을 한다는 것을 알 수 있다. 시뮬레이션 결과는 CPFS 내의 온도분포를 쉽게 이해할 수 있도록 3 차원 그래픽으로 나타냈으며, 상용소프트웨어 FEMLAB 으로 계산한 결과와 비교해서 개발된 모델과 계산 알고리즘의 정당성을 보였다. 맞이하고 있음을 볼 수 있다. 국내광업이 21C 급변하는 산업환경에 적응하여 생존하기 위해서는 각종 첨단산업에서 요구하는 소량 다품종의 원료광물을 적기에 공급 할 수 있는 전문화된 기술력을 하루속히 확보해야 하며, 이를 위해 고품위의 원료광물 확보를 위한 탐사 및 개발을 적극 추진하고 가공기술의 선진화를 위해 선진국과의 기술제휴 등 자원산업 글로벌화 정책이 절실히 요구되고 있음을 알 수 있다. 또한 삶의 질을 향상시키려는 현대인의 가치관에 부합하기 위해서는 각종 소비제품의 원료를 제공하는 광업의 본래 목적 이외에도 자연환경 훼손을 최소화하며 개발 할 수밖에 없는 구조적인 어려움에 직면할 수밖에 없다. 이처럼 국내광업이 안고 있는 여러 가지 난제들을 극복하기 위해서는 업계와 정부가 합심하여 국내광업 육성의 중요성을 재인식하고 새로운 마음가짐으로 관련 정책을 수립 일관성 있게 추진해 나가야 할 것으로 보인다.의 연구 결과를 요약하면 다음과 같다. 첫째, 브랜드 이미지와 서비스 품질과의 관계에서 브랜드이미지는 서비스 품질의 선행변수가 될 수 있음을 증명하였으며 4개 요인의 이미지 중 사풍이미지를 제외한 영업 이미지, 제품 이미지, 마케팅 이미지가 서비스 품질에 영향을 미치고 있음을 알 수 있다. 둘째, 지각된 서비스 품질과 가격 수용성과의 관계에서, 서비스 품질은 최소 가격에 신뢰서비스 요인에서 정의 영향을 미치고 있으나 부가서비스, 환경서비스에서는 역의 영향을 미침을 알수 있고, 최대 가격에 있어서는 욕구서비스 요인은 정의 영향을 미치지만 부가서비스의 경우에는 역의 영향을 미치고 있음을 알 수 있다. 셋째, 서비스품질과 재 방문 의도와의 관계에 있어서 서비스품질은 재 방문 의도에 영향을 미침을 알 수 있다. 따라서 브랜드 이미지는 서비스품질의 선행변수가 될 수 있으며, 서비스품질은 가격 수용성과 재방문 의도에 영향을 미치고 있음을 알 수

  • PDF

새로운 혼합알고리즘을 이용한 CPFS 내에서의 일어나는 동적 열전달의 수식화 및 해석 (Modeling and analysis of dynamic heat transfer in the cable penetration fire stop system by using a new hybrid algorithm)

  • 윤인섭;윤종필;권성필
    • 한국가스학회지
    • /
    • 제7권4호
    • /
    • pp.44-52
    • /
    • 2003
  • 본 논문은 원자력발전소 방화벽에 설치된 케이블관통부 충전시스템(CPFS: Cable penetration fire stop) 안에서 일어나는 동적 열전달 현상을 수식화하고, 새로운 혼합알고리즘을 이용해서 수치적으로 계산하여, 3차원 그래픽으로 나타내는 작업에 관한 연구이다. CPFS 내에서의 열전도 현상을 주어진 초기조건과 경계조건하에서 포물선 편미분방정식(Parabolic PDE)으로 수식화하였다. 계산을 단순화하기 위하여 전체 열 흐름을 z-축직선상에서의 일어나는 열전도 성분과 x-y-좌표 평면상에서 일어나는 열전도 성분으로 나누었다. z-축과 평행한 직선상에서 일어나는 열전도를 나타내는 PDE는 연속과완화법(SOR: Sequential over-relaxation)을 이용하여 유한불연속 점들에 대한 연립상미분방정식(ODE)으로 만들어서 풀었고, x-y-좌표 평면상에서 일어나는 열전도에 관한 PDE는 Galerkin 유한요소법(FEM: Finite element method)을 적용하여 ODE로 전환해서 풀었다. 여기서 시간과 공간의 함수인 온도는 각 직선상의 점들과 각 평면상의 요소절점들에 대해서 일정한 시간간격으로 초기온도와 경계온도를 업데이트하여 번갈아 가며 계산한다. 이러한 일련의 계산결과를 바탕으로 CPFS시스템 내에서의 온도분포의 동적인 변화를 계산해 낼 수 있었다. 결론적으로 관통하는 케이블이 CPFS시스템의 온도분포에 매우 중요한 역할을 한다는 것을 알 수 있었다. 시뮬레이션 결과는 CPFS내의 온도분포를 쉽게 이해할 수 있도록 3차원 그래픽으로 나타냈으며, 관통하는 케이블이 방화시스템의 온도분포에 매우 중요한 영향을 끼친다는 것을 알 수 있었다. 마지막으로 계산결과를 실험결과와 직접 비교함으로써, 개발된 모델과 계산 알고리즘의 정당성을 보였다.

  • PDF