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

검색결과 539건 처리시간 0.024초

NUMERICAL COUPLING OF TWO SCALAR CONSERVATION LAWS BY A RKDG METHOD

  • OKHOVATI, NASRIN;IZADI, MOHAMMAD
    • Journal of the Korean Society for Industrial and Applied Mathematics
    • /
    • 제23권3호
    • /
    • pp.211-236
    • /
    • 2019
  • This paper is devoted to the study and investigation of the Runge-Kutta discontinuous Galerkin method for a system of differential equations consisting of two hyperbolic conservation laws. The numerical coupling flux which is used at a given interface (x = 0) is the upwind flux. Moreover, in the linear case, we derive optimal convergence rates in the $L_2$-norm, showing an error estimate of order ${\mathcal{O}}(h^{k+1})$ in domains where the exact solution is smooth; here h is the mesh width and k is the degree of the (orthogonal Legendre) polynomial functions spanning the finite element subspace. The underlying temporal discretization scheme in time is the third-order total variation diminishing Runge-Kutta scheme. We justify the advantages of the Runge-Kutta discontinuous Galerkin method in a series of numerical examples.

Convergence Characteristics of the Crank-Nicolson-Galerkin Scheme for Linear Parabolic Systems

  • Cho, Jin-Rae;Ha, Dae-Yul;Kim, Tae-Jong
    • Journal of Mechanical Science and Technology
    • /
    • 제16권10호
    • /
    • pp.1264-1275
    • /
    • 2002
  • This paper is concerned with the investigation on the stability and convergence characteristics of the Crank-Nicolson-Galerkin scheme that is widely being employed for the numerical approximation of parabolic-type partial differential equations. Here, we present the theoretical analysis on its consistency and convergence, and we carry out the numerical experiments to examine the effect of the time-step size △t on the h- and P-convergence rates for various mesh sizes h and approximation orders P. We observed that the optimal convergence rates are achieved only when △t, h and P are chosen such that the total error is not affected by the oscillation behavior. In such case, △t is in linear relation with DOF, and furthermore its size depends on the singularity intensity of problems.

INTERNAL FEEDBACK CONTROL OF THE BENJAMIN-BONA-MAHONY-BURGERS EQUATION

  • Piao, Guang-Ri;Lee, Hyung-Chen
    • Journal of the Korean Society for Industrial and Applied Mathematics
    • /
    • 제18권3호
    • /
    • pp.269-277
    • /
    • 2014
  • A numerical scheme is proposed to control the BBMB (Benjamin-Bona-Mahony-Burgers) equation, and the scheme consists of three steps. Firstly, BBMB equation is converted to a finite set of nonlinear ordinary differential equations by the quadratic B-spline finite element method in spatial. Secondly, the controller is designed based on the linear quadratic regulator (LQR) theory; Finally, the system of the closed loop compensator obtained on the basis of the previous two steps is solved by the backward Euler method. The controlled numerical solutions are obtained for various values of parameters and different initial conditions. Numerical simulations show that the scheme is efficient and feasible.

Tracking Filter Design for a Maneuvering target Using Jump Processes

  • Lim, Sang-Seok
    • Journal of Electrical Engineering and information Science
    • /
    • 제3권3호
    • /
    • pp.373-384
    • /
    • 1998
  • This paper presents a maneuvering target model with the maneuver dynamics modeled as a jump process of Poisson-type. The jump process represents the deterministic maneuver(or pilot commands) and is described by a stochastic differential equation driven by a Poisson process taking values a set of discrete states. Employing the new maneuver model along with the noisy observations described by linear difference equations, the author has developed a new linear, recursive, unbiased minimum variance filter, which is structurally simple, computationally efficient, and hence real-time implementable. Futhermore, the proposed filter does not involve a computationally burdensome technique to compute the filter gains and corresponding covariance matrices and still be able to track effectively a fast maneuvering target. The performance of the proposed filter is assessed through the numerical results generated from the Monte-Carlo simulation.

  • PDF

GLOBAL ROBUST STABILITY OF TIME-DELAY SYSTEMS WITH DISCONTINUOUS ACTIVATION FUNCTIONS UNDER POLYTOPIC PARAMETER UNCERTAINTIES

  • Wang, Zengyun;Huang, Lihong;Zuo, Yi;Zhang, Lingling
    • 대한수학회보
    • /
    • 제47권1호
    • /
    • pp.89-102
    • /
    • 2010
  • This paper concerns the problem of global robust stability of a time-delay discontinuous system with a positive-defined connection matrix under polytopic-type uncertainty. In order to give the stability condition, we firstly address the existence of solution and equilibrium point based on the properties of M-matrix, Lyapunov-like approach and the theories of differential equations with discontinuous right-hand side as introduced by Filippov. Second, we give the delay-independent and delay-dependent stability condition in terms of linear matrix inequalities (LMIs), and based on Lyapunov function and the properties of the convex sets. One numerical example demonstrate the validity of the proposed criteria.

이동 CRANE의 제어에 관한 연구 (A Study on Control of Mobile Cranes)

  • 김상봉;신민생;김환성;정용길
    • 한국정밀공학회지
    • /
    • 제8권2호
    • /
    • pp.47-56
    • /
    • 1991
  • The specifications needed for the mobile cranes are summarized as the following : 1) there may be not occured the oscillation of the cargo at unloading point. 2)the required time from departure point to destination point may be as short as possible. 3) there may be not a collapse of cargo caused by the oscillation in the course that the crago is mobilling. In this paper, the linear fractional transformation method is adopted as a method in order to improve the above mentioned problems. A design method of servo system is developed by modifying Davison's method for the case that the homogeneous differential equations of reference input and disturbance are different types. The real time control of a mobile crane system is implemented by 16bits microcomputer with A/D and D/A converters to illustrate the application of the adopted method. The experimental results for the three types of the design methods; linear fractional transformation method, servo system design method and optimal control method are shown for the comparison.

  • PDF

일정 일반속력으로 구동되는 구속 다물체계의 선형화기법 및 진동해석 (Linearization Method and Vibration Analysis of a Constrained Multibody System Driven by Constant Generalized Speeds)

  • 최동환;박정훈;유홍희
    • 한국소음진동공학회:학술대회논문집
    • /
    • 한국소음진동공학회 2001년도 추계학술대회논문집 II
    • /
    • pp.725-730
    • /
    • 2001
  • This paper presents a vibration analysis method for constrained mechanical systems driven by constant generalized speeds. Equilibrium positions are obtained first and vibration analysis are performed around the positions. The method developed in this paper employs partial velocity matrix to obtain a minimum number of differential equations. To verify the accuracy of the proposed algorithm, linear vibration analyses of two numerical examples are performed and the results are compared with results from a commercial program or previous literature.

  • PDF

PARAMETER DEPENDENCE OF SMOOTH STABLE MANIFOLDS

  • Barreira, Luis;Valls, Claudia
    • 대한수학회지
    • /
    • 제56권3호
    • /
    • pp.825-855
    • /
    • 2019
  • We establish the existence of $C^1$ stable invariant manifolds for differential equations $u^{\prime}=A(t)u+f(t,u,{\lambda})$ obtained from sufficiently small $C^1$ perturbations of a nonuniform exponential dichotomy. Since any linear equation with nonzero Lyapunov exponents has a nonuniform exponential dichotomy, this is a very general assumption. We also establish the $C^1$ dependence of the stable manifolds on the parameter ${\lambda}$. We emphasize that our results are optimal, in the sense that the invariant manifolds are as regular as the vector field. We use the fiber contraction principle to establish the smoothness of the invariant manifolds. In addition, we can also consider linear perturbations, and thus our results can be readily applied to the robustness problem of nonuniform exponential dichotomies.

건물의 3차원 구조체에 대한 전열해석 프로그램 개발 중 서로 다른 열전도율을 갖는 복합재질 3차원 구조의 비정렬 격자에 대한 전산해석 방법 (Numerical heat transfer analysis methodology for multiple materials with different heat transfer coefficient in unstructured grid for development of heat transfer analysis program for 3 dimensional structure of building)

  • 이주희;장진우;이현균;이용준;이규성
    • KIEAE Journal
    • /
    • 제16권1호
    • /
    • pp.81-87
    • /
    • 2016
  • Purpose: Heat transfers phenomena are described by the second order partial differential equation and its boundary conditions. In a three-dimensional structure of a building, the heat transfer phenomena generally include more than one material, and thus, become complicate. The analytic solutions are useful to understand heat transfer phenomena, but they can hardly be applied in engineering or design problems. Engineers and designers have generally been forced to use numerical methods providing reliable results. Finite volume methods with the unstructured grid system is only the suitable means of the analysis for the complex and arbitrary domains. Method: To obtain an numerical solution, a discretization method, which approximates the differential equations, and the interpolation methods for temperature and heat flux between two or more materials are required. The discretization methods are applied to small domains in space and time, and these numerical solutions form the descretized equations provide approximated solutions in both space and time. The accuracy of numerical solutions is dependent on the quality of discretizations and size of cells used. The higher accuracy, the higher numerical resources are required. The balance between the accuracy and difficulty of the numerical methods is critical for the success of the numerical analysis. A simple and easy interpolation methods among multiple materials are developed. The linear equations are solved with the BiCGSTAB being a effective matrix solver. Result: This study provides an overview of discretization methods, boundary interface, and matrix solver for the 3-dimensional numerical heat transfer including two materials.

경계조건 변화에 따른 Pasternak 지반으로 지지된 원호형 띠기초의 휨 자유진동에 관한 변수연구 (Parametric Studies of Flexural Free Vibrations of Circular Strip Foundations with Various End Constraints Resting on Pasternak Soil)

  • 이병구;이광범;강희종;윤희민
    • 한국소음진동공학회논문집
    • /
    • 제17권9호
    • /
    • pp.835-846
    • /
    • 2007
  • 이 논문은 Pasternak 지반으로 지지된 변화폭 원호형 띠기초의 휨 자유진동에 관한 연구이다. 단면폭은 띠기초의 중앙점을 중심으로 대칭 일차 함수를 갖는 변화폭으로 선정하였다. 이 연구에서 지반은 전단층을 갖는 탄성지반인 Pasternak 지반으로 모형화 하였고 회전관성과 전단변형을 고려한 곡선 띠기초의 자유진동을 지배하는 상미분방정식을 유도하였다. 이 상미분방정식을 수치해석하여 고유진동수 및 진동형을 산출하였다. 수치해석 예에서는 회전-회전, 회전-고정 및 고정-고정의 3개의 지점조건을 고려하였다. 띠기초의 무차원 변수들이 고유진동수에 미치는 영향을 분석하였다.