• Proceedings of the Korea Water Resources Association Conference
• /
• 2004.05b
• /
• pp.660-665
• /
• 2004

본 연구에서는 비선형적인 메커니즘을 갖는 수문현상의 불확실성을 해석하고자 하는 목적으로 새로운 개념의 지배방정식이 유도된다. 제안된 모형의 불확실성은 토양 특성치의 공간적 변동성에 기인하고 있는 것으로 가정하여, 유도된 방정식은 Fokker-Planckl 방정식의 형태를 가지고 있다. 실제 유역단위에서 토양 내 수분 흐름의 연직방향 흐름을 모의하기 위해 미소단위에서 유도된 Richards 방정식은 토양의 공간적 변동성으로 말미암아 불확실한 매개변수를 갖는 비선형 추계학적 편미분방정식의 형태를 갖게 된다. 이는 먼지 수직 방향적분을 통하여 단순화된 비선형 추계학적 상미분방정식으로 전환되고, 이렇게 전환된 비선형 추계학적 상미분방정식은 다시 추계학적 Liouville 방정식을 이용하여 선형 추계학적 편미분방정식으로 전환되어진다. 최종적으로 cumulant 급수방법을 이용하여 상기 방정식을 선형 결정론적 편미분방정식으로 전환시킴으로써, 강우 시 토양 내 수분 침투현상을 모형화할 경우 유역단위에서 토양의 공간적 변동성을 설명할 수 있는 지배방정식을 유도할 수 있다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.4
• /
• pp.427-436
• /
• 2000

This study presents the nonlinear responses of a hinged-clamped beam under broadband random excitation. By using Galerkin's method the governing equation is reduced to a system or nonautonomous nonlinear ordinary differential equations. The Fokker-Planck equation is used to generate a general first-order differential equation in the joint moments of response coordinates. Gaussian and non-Gaussian closure schemes are used to close the infinite coupled moment equations. The closed equations are then solved for response statistics in terms of system and excitation parameters. The case of two mode interaction is considered in order to compare it with the case of three mode interaction. Monte Carlo simulation is used for numerical verification.

• Journal of the KSME
• /
• v.23 no.3
• /
• pp.200-206
• /
• 1983

FAM의 기본적인 구상은 해석 하고자하는 선형 또는 비선형 편미분 방정식을 국부적으로 해석 적인 해를 구하여 이용하자는 것이다. 그러기 위하여 유한차분법(FDM)과 유한변분법(FEM)에 서와 같이 전체유동장을 작은 요소로 나누고 그 요소 내에서 국부해를 구한 다음 이들 요소를 중첩시킴으로써 각 요소의 미지수에 대한 대수식을 얻어서 수치해를 구하자는 것이다. 그러나 FDM에서와 같이 국부요소에서 미분항을 구하지 않고, FEM 에서와 같이 요소에서 형상함수를 도입하지 않는 상태에서 해석적인 해를 구하고 있기 때문에 수치해석에서 얻어지는 미분양들은 비교적 정확하게 구해진다. 따라서 Navier-Stokes 방정식이나 에너지 방정식에서 최고차항이 작은 파라메타, 즉 레이놀즈수나 피크리수의 역수로 곱하여서 있는 경우에도 안정된 해를 구할 수 있다고 알려져 있다. 요소자체의 계수를 구하는 데는 계산시간이 많이 소요되지만 수치해석 상의 안정성이나 수렴성이 좋기 때문에 전체계산시간은 오히려 적게 걸리는 경우도 있다고 한다.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1991.04a
• /
• pp.147-152
• /
• 1991

양단이 고정된 보가 변형할 때에는 중간 평면의 신장을 수반하게 된다. 운동 의 진폭이 증가함에 따라 이 신장이 보의 동적 응답에 미치는 영향은 심각 하게 된다. 이러한 현상은 응력과 변형도와의 관계가 선형적이라 하더라도 변형도와 변위와의 관계식은 비선형이 되며 결국은 보의 비선형 운동방정식 을 낳게된다. 보는 연속계이긴하지만 근사를 위하여 다자유도계로 간주할 수 있다. 비선형 다자유도계에 있어서는 선형화된 계의 고유진동수끼리 적절한 관계를 가질 때 내부공진이 발생할 수 있다. 양단이 고정된 곧은 보의 비선 형 동적응답이 그동안 많이 연구되어 오고 있으며, 집중질량을 가지고 직각 으로 굽은 보의 해석을 위하여 내부공진을 고려한 해석적 혹은 실험적 연구 가 이루어져 왔다. 그중에서도 Nayfeh등은 조화가진 하의 핀과 꺾쇠로 고정 된(hinged-clamped) 보의 정상상태응답을 해석하기 위해 두 모우드 사이의 내부공진을 고려하였다. 이 연구에서는 세 모우드 사이의 내부공진을 고려하 여 강제진행 중인 보의 비선형 해석을 다루고자 한다. 이 문제에 관심을 갖 게 된 동기는 "연속계의 비선형 해석에서 더 많은 모우드를 포함시키면 어 떤 결과를 낳게 될 것인가\ulcorner"라는 질문에서 생겨난 것이다. 갤러킨 법을 이용 하여 비선형 편미분 방정식과 경계 조건으로 표현되는 이 문제를 연립 비선 형 상미분 방정식으로 변환한다. 다중시간법(the method of multiple scales) 을 이용하여 이 상미분 방정식을 정상상태에서의 세 모우드의 진폭과 위상 에 대한 연립비선형 대수방정식으로 변환한다. 이 대수방정식을 수치적으로 풀어서 정상상태 응답을 구하고 Nayfeh등의 결과와 비교한다. 결과와 비교한다. studies, the origin of ${\alpha}$$_1$peak was attributed to the detrapping process form trap with 2.88[eV] deep of injected space charge from the chathode in the crystaline regions. The origin of ${\alpha}$$_2$ peak was regarded as the detrapping process of ions trapped with 0.9[eV] deep originated from impurity-ion remained in the specimen during production process of the material, in the crystalline regions. The origin of ${\beta}$ peak was concluded to be due to the depolarization process of "C=0"dipole with the activation energy of 0.75[eV] in the amorphous regions. The origin of ${\gamma}$ peak was responsible to the process combined with the depolarization of "CH$_3$", chain segment, with the activation energy of carriers from the shallow trap with 0.4[eV], in he amorp

• Publications of The Korean Astronomical Society
• /
• v.5 no.1
• /
• pp.40-64
• /
• 1990

복사전달식, 통계평형식 및 전하 입자 보존식을 동시에 만족시키는 비열평형상태의 대기에 대하여 그 방출원자스펙트럼을 수치계산하였다. 등온, 중압의 비열평형 태양홍염을 대상모델로 선정하였는데 여기에 적응한 제한조건중 복사전달에 관련된 편미분방정식은 3점 근사 차분법에 의해 정리하였고 홍염중심에 대칭성을 가정하여 경제조건을 부여하였으며 대기의 물리상태에 관련된 비선형 연립방정식의 해는 완전선형화 기법을 통한 퓨트리에 소거법을 이용해 구하였다.

• Proceedings of the KIEE Conference
• /
• 1999.07e
• /
• pp.2441-2443
• /
• 1999

광섬유에서 전송되는 신호의 예측을 위하여 편미분방정식인 비선형 슈래딩거 방정식(Nonlinear Schrodinger Equation, NLSE)을 단계분할 유한 요소법(Split-Step Finite Element Method, SS-FEM)을 적용하여 해석하였다. 수치결과를 해석적인 해가 알려진 솔리톤의 해로부터 전송되는 거리의 증가에 따라 각 단계마다 오차를 계산하였으며, 그 결과를 단계분할 푸리에법(Split-Step Fourier Method, SS-FM)에 의한 수치해와도 비교하였다.

• Clean Technology
• /
• v.10 no.1
• /
• pp.9-17
• /
• 2004

Direct methanol fuel cell(DMFC) with proton exchange membrane (PEM) has advantages over the conventional power source (e.g. vehicle). DMFC, however, has a problem to be solved such as methanol crossover, high anodic overpotential and limiting current density, etc. The physicochemical phenomena in DMFC can be described by coupled PDEs (partial differential equations), which can be solved by a PDE solver. In this paper, we utilized a commercial software FEMLAB to solve the PDEs. The FEMLAB is one of the software programs available which are developed as a solver for building physics problems based on PDEs and is designed to simulate systems of coupled PDEs which may be 1D, 2D, 3D, non-liner and time dependent. We performed simulation using the Tafel equation as an electrochemical reaction model to analyze methanol concentration profile in DMFC system. We confirm that the rapid decrease of methanol concentration at anodic catalyst layer with the increase of the current density is a main reason of the low performance in DMFC through simulation results.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.1 no.2
• /
• pp.82-88
• /
• 1977

1장이 축방향으로 자유로이 움직일 수 있는 단선지시보의 진동에 있어서 그 진폭이 크고, 축방향의 관성력과 전단변형을 고려할때 의 진동을 살폈다. 운동방정식을 처짐을 독립변수로 하는 비선형편미분방정식으로 표시하고, model expansion과 Galerkin 방법에 의해서 비정형연립상미분방정식으로 변형한 다음에 Perturbation method of muthod of multiple scale로 유사해를 구하였다. 또한 보의 진동수-진폭의 관계에 대한 일반적 표현을 구하고, 간단한 구체적 해에 대하여 1차근사해와 진동수-진폭관계를 계산하여 이미 이루어진 연구결과와 비교하였다.

• The Transactions of the Korean Institute of Power Electronics
• /
• v.6 no.2
• /
• pp.115-124
• /
• 2001

In this paper, authors propose a new nonlinear rotor flux observer for rotor field oriented control of an induction motor which is designed based on the extended Luenberger Observer theory. Extended Luenberger Observer requires minimal solution of nonlinear partial differential equation on its coordinate transformation and linearization needed on a nonlinear observer design in general. The proposed rotor flux observer is derived from the 2 phase model of induction motor on the orthogonal coordination and it has the reduce gain matrix. Simulation and experimentation were performed under the conventional indirect vector control and direct vector control with the proposed observer at different rotor resistance. Simulation results show that the convergence of the proposed observer is influenced by the chosen eigenvalues. Experimentation results on load operation show the direct vector control with the proposed observer is better than the indirect vector control to maintain the characteristics of the vector control.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.3
• /
• pp.80-86
• /
• 1994

A numerical method for simulations of inviscid incompressible flow fields around a ship advancing on the free surface is developed. A body fitted coordinate system, generated by numerically solving elliptic type partial differential equations is used to conform the ship and free surface configurations. Three dimensional Euler equations transformed to the non-staggered body fitted coordinate system are discretised by finite difference method. Time and spatial derivatives are discretised by forward and centered differencings, respectively, and artificial dissipations are added to discretised convection terms for improvements of numerical stability. At each time steps, free surface elevations are recomputed to satisfy nonlinear free surface conditions. Poisson equations for pressure field are solved iteratively and the velocity field for next time step is extrapolated. To verify the developed numerical method, flow fields around a Wigley model are simulated(Fn=0.250-0.408) and compared with experimental data to show good agreements.