• 제목/요약/키워드: partial differential equations

검색결과 518건 처리시간 0.029초

이동최소제곱 유한차분법을 이용한 응력집중문제 해석(I) : 고체문제의 정식화 (Analysis of Stress Concentration Problems Using Moving Least Squares Finite Difference Method(I) : Formulation for Solid Mechanics Problem)

  • 윤영철;김효진;김동조;윙 캠 리우;테드 벨리치코;이상호
    • 한국전산구조공학회논문집
    • /
    • 제20권4호
    • /
    • pp.493-499
    • /
    • 2007
  • 본 연구에서는 미분 가능한 함수가 Taylor 전개로 표현되고 그 계수들은 주어진 함수와 미분에 대한 근사값을 제공할 수 있다는 점에 착안하여 m차 Taylor 다항식을 구성하고 이동최소제곱법을 이용하여 그 계수들을 구했다. 계산된 근사함수와 미분을 콜로케이션 개념을 바탕으로 균열 문제를 포함하는 고체문제에 대한 지배 미분방정식에 적용하여 차분식 형태의 이산화된 계방정식을 구성하였다. 본 연구의 해석기법은 격자망(grid)에 의존적이고 근사함수가 없는 유한차분법과 형상함수의 미분과 약형식의 적분산정, 필수경계조건 처리가 어려운 Galerkin법 기반의 무요소법의 단점을 효과적으로 극복한 새로운 수치기법이다.

전산유체역학 응용에서의 효율적인 최적 2차 변수 계산 경로 추정 기법 (An Efficient Method for Estimating Optimal Path of Secondary Variable Calculation on CFD Applications)

  • 이중연;김민아;허영주
    • 한국콘텐츠학회논문지
    • /
    • 제16권12호
    • /
    • pp.1-9
    • /
    • 2016
  • 전산유체역학은 유체현상을 기술하는 미분방정식의 근사해를 컴퓨터를 이용하여 풀고 해석하는 학문으로, 다양한 종류의 변수의 계산을 필요로 한다. 대용량의 유동해석 데이터의 경우, 스토리지의 제약으로 계산한 변수들 중 필수적인 변수만을 저장하고, 데이터 분석 시점에 필요한 2차 변수를 계산하는 경우가 많다. 본 논문에서는 전산유체역학 응용에서 많이 사용하는 2차 변수의 종류를 정리하고, 임의의 초기 변수가 주어졌을 때 최적의 2차 변수 계산 경로를 효율적으로 추정하기 위한 방법으로 2차 변수 종속 그래프를 일반적인 유향 그래프로 변환하는 기법과 이를 이용한 최단 경로 탐색 기법을 소개한다. 또한 제안하는 기법을 실제 데이터 분석 및 가시화 도구에 적용하여 효용성을 측정하였다.

동심 원통형 용기내에서의 플라즈마 회전 (Plasma Rotation in Plasma Centrifuge with an Annular Gap)

  • Hue Yeon Lee;Sang Hee Hong
    • Nuclear Engineering and Technology
    • /
    • 제14권2호
    • /
    • pp.78-85
    • /
    • 1982
  • 정상상태의 플라즈마의 이론적 해석을 통해서 플라즈마의 회전과 동위원소 분리기로서의 적합성을 분석하였다. 이 장치는 두개의 동심원통형 전극과 이들 사이의 원통형 공동으로 구성되었으며, 축 방향으로 외부자장이 걸려 있다. 두 전극사이에 생성되는 전류밀도는 전기방전의 형태로 동위원소 혼합물로부터 플라즈마를 생성하고, 자장과 교차되어 발생하는 Lorentz힘에 의해서 플라즈마를 회전시킨다. 자기 유체역학 방정식을 바탕으로 이 계를 설명하는 두개의 연립편미분방정식을 얻었고, 네 경계조건을 사용하여 Fourier-Bessel로 표현된 이차원적 전류밀도와 속도분포의 해를 얻었다. 실제로 가능한 조건하에서 플라즈마 회전속도는 $10^4$m/sec 정도에 달하고, Hartmann수가 커짐에 따라 플라즈마회전 속도도 커진다. 이 같은 고속의 회전속도를 감안해 볼때 플라즈마 원심분리기는 기계적으로 회전되는 가스원심분리기보다 훨씬 높은 효율을 가지게 될 것이다.

  • 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

대두포장(大豆圃場)에서 수분흡수(水分吸收)에 관(關)한 수치해석학적(數値解析學的) 모형(模型) (Numerical Simulation of Water Uptake of Soybean Field)

  • 정영상
    • 한국토양비료학회지
    • /
    • 제14권1호
    • /
    • pp.8-16
    • /
    • 1981
  • Ohm의 법측(法測)을 이용(利用)한 물 이동(移動) 방정식(方程式)의 해(解)를 수치해석법(數値解析法)으로 구(求)하여 대두(大豆)의 물 소모양상(消耗樣相)을 추정(推定)하여 실측자료(實測資料)와 비교(比較)하였다. 본(本) 모형(模型)은 고속전산기(高速電算器)를 사용(使用)하기 위(爲)하여 FOPTRAN으로 쓰여졌으며 모형(模型)의 검정(檢定)은 미(美) Iowa 주서부농업시험장(洲西部農業試驗場)에서 얻어진 자료(資料)를 사용(使用)하여 시도(試圖)되었다. 모형검정치(模型檢定値)와 실측치간(實測値間)에 잘 맞는 경향(傾向)이었으나 다소(多少)의 오차(誤差)가 시기적(時期的)으로 인정(認定)되었다. 계산(計算)된 토양증발량(土壤蒸發量)은 작물생육초기(作物生育初期)에는 약(約) 2mm/day이었으나 후기(後期)에는 0.4mm/day로 감소(減少)하였으며, 한발기(旱魃期)에 0.1mm/day까지 감소(減少)되었다. 증발량(蒸發量)은 최대(最大) 5mm/day이었다. 토심(土深) 180cm와 210cm 층위(層位)에서 이루어지는 지하배수량(地下排水量)은 약(約) 0.2mm/day이었고, 후기(後期)에는 역전되어 모관상승(毛管上昇)을 보였으며 이 양(量)은 약 0.07mm/day로 근계(根系)의 물 소모량(消耗量)을 충족(充足)시키지는 못하였다. 따라서 심토(深土)의 물 이용(利用)을 위(爲)해서는 근계(根系)의 확장(擴張)이 더 중요(重要)하다.

  • PDF

미분구적법(DQM)을 사용한 비대칭 곡선 보의 신장 좌굴해석 (Extensional Buckling Analysis of Asymmetric Curved Beams Using DQM)

  • 강기준
    • 한국산학기술학회논문지
    • /
    • 제22권4호
    • /
    • pp.594-600
    • /
    • 2021
  • 곡선 보는 철교 그리고 자동차와 갈은 구조물의 구성으로 널리 사용되어왔다. 많은 연구자들의 관심분야인 이러한 구조물의 안정성 거동 해석분야는 괄목할 만한 성과가 있어 왔다. 곡선 보 구조물의 기하학적 구조 및 물성치가 탄성 및 강성에 미치는 영향을 분석하기 위하여 정역학적 동역학적 해석이 필요하다. 그러나 구조물의 복잡성 때문에 어떠한 경계조건에서도 엄밀해를 얻기가 매우 어렵다. 전통적으로 미분방정식의 해법은 유한차분법 혹은 유한요소법으로 해결해왔으나 이러한 방법들은 때론 복잡한 비선형 구조물에는 과도한 컴퓨터 용량사용과 복잡한 알고리즘 프로그램을 요구한다. 이러한 어려움을 해결하기 위해 미분구적법(DQM)이 여러 분야에 사용되어왔다. 본 연구에서는 복잡한 편미분 방정식의 해를 구하기 위하여 미분구적법이 사용되었다. 중면 신장을 고려한 등분포 하중 하에서 선형으로 변하는 비대칭 곡선 보의 내평면 신장 좌굴의 지배방정식을 유도하였고, DQM을 이용하여 지배방정식의 해를 구하였다. 다양한 열림 각, 경계 조건, 그리고 파라미터에 의한 임계하중을 계산하였다. DQM 결과는 비교 가능한 엄밀해와 비교하였고 DQM은 적은 격자점을 사용하고도 정확성을 보여주었다. 예를 들어 열림 각이 180°인 비 신장 고정단 곡선 보의 경우, 엄밀해의 임계하중 값은 8.0이고 DQM의 임계하중 값은 7.98로, 오차가 0.3% 미만 이었다. 곡선 보의 내평면 비 신장 임계하중도 계산하였고, 두 이론을 상호 비교 분석하였다. 아크축의 중면 신장을 고려한 연구는 곡선 보의 임계하중에 중대한 영향을 미치는 것을 보여준다.

대형이산 행렬 시스템의 초대형병렬컴퓨터에서의 해법을 위한 병렬준비 행렬의 비교 (Comparison of Parallel Preconditioners for Solving Large Sparse Linear Systems on a Massively Parallel Machine)

  • 마상백
    • 한국정보처리학회논문지
    • /
    • 제2권4호
    • /
    • pp.535-542
    • /
    • 1995
  • 이 논문에서 우리는 CM-5와 같은 초대형병렬컴퓨터에서 대형 이산선형체제를 풀기 위한 준비행렬로써 두 가지를 소개한다. 대다수의 초대형병렬컴퓨터들은 프로세서간의 통신을 메세지패씽(messagepassing)에 의존하는데 현재의 기술수준하에서는 이 통신속 도가 실수계산속도에 비해 매우 느리므로 종래의 메모리공유컴퓨터에서와는 달리 데이 터통신량을 최소화하는 알고리듬이 요구된다. 블록 SOR에 다중색채기법을 가미한 알고 리듬이 그 한 예로써 우리는 이를 CM-5에서 구현한 결과 N=512x512 행렬에서 프로세서 의 수가 16에서 512의 범위 하에서 50%의 효율을 실현하였다. 반면 종래의 효율적인 병렬 준비행렬로 알려진 AKI알고리듬은 방대한 량의 데이터통신 때문에 매우 열등한 결과를 보여준다.

  • PDF

다양한 적층 배열을 갖는 복합 신소재 판 구조물의 유한차분 비선형 해석 (Finite Difference Nonlinear Analysis of Composite Plate Structures with Various Layer Sequences)

  • 이상범;이상열;이래철
    • 한국구조물진단유지관리공학회 논문집
    • /
    • 제9권4호
    • /
    • pp.159-168
    • /
    • 2005
  • 본 연구는 다양한 적층 배열을 갖는 비등방성을 보이는 첨단 복합 신소재 판구조물의 유한 차분 비선형 해석을 수행한다. 복잡한 편미분 방정식으로 표현되는 역학문제들을 수치해석 하는 경우 본 연구에서 사용한 유한차분법은 유한요소법에 비하여 체눈 생성 및 수치적분 과정을 피할 수 장점을 갖는다. 유한차분법을 이용한 많은 연구들은 단지 에너지 방법을 사용한 고정 혹은 단순 경계조건에 대하여 수행되었다. 그러나 이러한 접근방법은 자유경계에 대하여 불가피하게 발생하는 가상점 문제를 충분히 만족시킬 수 없다. 그러므로 본 연구에서는 임의의 경계조건을 갖는 비등방성 복합 적층한의 비선형 거동 문제를 보다 효과적으로 해결할 수 있는 유한차분식을 정식화 하였다. 적층 배열 변화를 비롯한 다양한 매개변수에 대하여 본 연구에서 제안한 접근방법을 사용하여 적층판의 복잡한 비선형 거동을 분석하였다.

두개의 제어기를 사용한 건물 내부의 온도변화와 에너지소비량을 계산하기 위한 해석적 연구 (A study on the analytical method for calculating the inside air temperature transient and energy consumption load of the building using two different controllers)

  • 한규일
    • 수산해양기술연구
    • /
    • 제48권1호
    • /
    • pp.82-90
    • /
    • 2012
  • Four different buildings having various wall construction are analyzed for the effect of wall mass on the thermal performance and inside building air and wall temperature transient and also for calculating the energy consumption load. This analytical study was motivated by the experimental work of Burch et al. An analytical solution of one-dimensional, linear, partial differential equations is obtained using the Laplace transform method, Bromwich and modified Bromwich contour method. A simple dynamic model using steady state analysis as simplified methods is developed and results of energy consumption loads are compared with results obtained using the analytical solution. Typical Meteorological Year data are processed to yield hourly average monthly values. This study is conducted using weather data from two different locations in Korea: Daegu having severe weather in summer and winter and Jeju having mild weather almost all year round. There is a significant wall mass effect on the thermal performance of a building in mild weather condition. Buildings of heavyweight construction with insulation show the highest comfort level in mild weather condition. A proportional controller provides the higher comfort level in comparison with buildings using on-off controller. The steady state analysis gives an accurate estimate of energy load for all types of construction. Finally, it appears that both mass and wall insulation are important factors in the thermal performance of buildings, but their relative merits should be decided in each building by a strict analysis of the building layout, weather conditions and site condition.

Thermal diffusion and diffusion thermo effects on an unsteady heat and mass transfer magnetohydrodynamic natural convection Couette flow using FEM

  • Raju, R. Srinivasa;Reddy, G. Jithender;Rao, J. Anand;Rashidi, M.M.
    • Journal of Computational Design and Engineering
    • /
    • 제3권4호
    • /
    • pp.349-362
    • /
    • 2016
  • The numerical solutions of unsteady hydromagnetic natural convection Couette flow of a viscous, incompressible and electrically conducting fluid between the two vertical parallel plates in the presence of thermal radiation, thermal diffusion and diffusion thermo are obtained here. The fundamental dimensionless governing coupled linear partial differential equations for impulsive movement and uniformly accelerated movement of the plate were solved by an efficient Finite Element Method. Computations were performed for a wide range of the governing flow parameters, viz., Thermal diffusion (Soret) and Diffusion thermo (Dufour) parameters, Magnetic field parameter, Prandtl number, Thermal radiation and Schmidt number. The effects of these flow parameters on the velocity (u), temperature (${\theta}$) and Concentration (${\phi}$) are shown graphically. Also the effects of these pertinent parameters on the skin-friction, the rate of heat and mass transfer are obtained and discussed numerically through tabular forms. These are in good agreement with earlier reported studies. Analysis indicates that the fluid velocity is an increasing function of Grashof numbers for heat and mass transfer, Soret and Dufour numbers whereas the Magnetic parameter, Thermal radiation parameter, Prandtl number and Schmidt number lead to reduction of the velocity profiles. Also, it is noticed that the rate of heat transfer coefficient and temperature profiles increase with decrease in the thermal radiation parameter and Prandtl number, whereas the reverse effect is observed with increase of Dufour number. Further, the concentration profiles increase with increase in the Soret number whereas reverse effect is seen by increasing the values of the Schmidt number.