• Proceedings of the KIEE Conference
• /
• 1979.08a
• /
• pp.119-120
• /
• 1979

A linear state equation of the first order differential form relating the load-frequency dynamic characteristics of interconnected power systems was derived for use in computer simulation. A now solution of the algebraic matrix riccati equation for application in quadratic optimal controllor and least-square state estimator dermination was developed. The program for a dynamic state equation for two interconnected control areas was developed. The optimized load-frequency deviation was analysed and a numerical analysis was tried based on the computer simulation. It was shown that the dynamic response of th loed-frequency could be optimized with weighting factors IR and Q. The result was that the load-frequency and the tie-line deviation were visibly reduced.

• Structural Engineering and Mechanics
• /
• v.14 no.3
• /
• pp.263-285
• /
• 2002

This paper presents a numerical procedure for solving initial-value problems using the special functions which belong to a class of Rvachev's basis functions $R_{bf}$ based on algebraic and trigonometric polynomials. Because of infinite derivability of these functions, derivatives of all orders, required by differential equation of the problem and initial conditions, are used directly in the numerical procedure. The accuracy and stability of the proposed numerical procedure are proved on an example of a single degree of freedom system. Critical time step was also determined. An algorithm for solving multiple degree of freedom systems by the collocation method was developed. Numerical results obtained by $R_{bf}$ functions are compared with exact solutions and results obtained by the most commonly used numerical procedures for solving initial-value problems.

• Journal of applied mathematics & informatics
• /
• v.31 no.1_2
• /
• pp.299-309
• /
• 2013

In this paper, we introduce a new numerical technique which we call fractional Chebyshev finite difference method (FChFD). The algorithm is based on a combination of the useful properties of Chebyshev polynomials approximation and finite difference method. We tested this technique to solve numerically fractional BVPs. The proposed technique is based on using matrix operator expressions which applies to the differential terms. The operational matrix method is derived in our approach in order to approximate the fractional derivatives. This operational matrix method can be regarded as a non-uniform finite difference scheme. The error bound for the fractional derivatives is introduced. The fractional derivatives are presented in terms of Caputo sense. The application of the method to fractional BVPs leads to algebraic systems which can be solved by an appropriate method. Several numerical examples are provided to confirm the accuracy and the effectiveness of the proposed method.

• 제어로봇시스템학회:학술대회논문집
• /
• 2004.08a
• /
• pp.1769-1775
• /
• 2004

In presence of collision between two rigid bodies, they exhibit impulsive behavior to generate physically feasible state. When the frictional impulse is involved, collision resolution can not be easily made based on a simple Newton's law or Poisson's law, mainly due to possible change of collision mode during collision, For example, sliding may change to sticking, and then sliding resumes. We first examine two conventional methods: the method of mode evolution by differential equation, and the other by linear complementarity programming. Then, we propose a new method for mode evolution by solving only algebraic equations defining mode changes. Further, our method attains the original nonlinear impulse cone constraint. The numerical simulation will elucidate the advantage of the proposed method as an alternative to conventional ones.

• International Journal of Automotive Technology
• /
• v.2 no.3
• /
• pp.117-122
• /
• 2001

For a dynamic analysis of a constrained multibody system, it is necessary to have a routine for satisfying kinematic constraints. LU decomposition scheme, which is used to divide coordinates into dependent and independent coordinates, is efficient but has great difficulty near the singular configuration. Other method such as the projection matrix, which is more stable near a singular configuration, takes longer simulation time due to the large amount of calculation for decomposition. In this paper, the row space and the null space of the Jacobian matrix are proposed by using the pseudo-inverse method and the projection matrix. The equations of the motion of a system are replaced with independent acceleration components using the null space of the Jacobian matrix. Also a new hybrid method is proposed, combining the LU decomposition and the projection matrix. The proposed hybrid method has following advantages. (1) The simulation efficiency is preserved by the LU method during the simulation. (2) The accuracy of the solution is also achieved by the projection method near the singular configuration.

• International Journal of Automotive Technology
• /
• v.8 no.3
• /
• pp.289-298
• /
• 2007

A nonlinear dynamic simulation model from cranking to idle speed is developed to optimize the cold start process of a diesel engine. Physically-based first order nonlinear differential equations and some algebraic equations describing engine dynamics and starter motor dynamics are used to model the performance of cold starting process which is very complex and involves many components including the cold start aiding method. These equations are solved using numerical schemes to describe the starting process of a diesel engine and to study the effects of cold starting parameters. The validity of this model is examined by a cold start test at $-20^{\circ}C$. Using the developed model the effects of the important starting variables on the cold starting processes were investigated. This model can be served as a tool for designing computer aided control systems that improve cold start performance.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.145-151
• /
• 1999

An Interactive grid generation program(KGRID) with graphical user interface(GUI) has been improved. KGRID works on the UNLX environment and GUI has been implemented with OSF/Motif and X Toolkit and the graphics language is Open GL for visualization of the 3D objects. It supports more convenient user environment to generate 2D and 3D multi-block structured grid systems. It provides various useful field grid generation methods, which are the algebraic methods, the elliptic partial differential equations method and the predictor-corrector method. It also supports 3D surface grid generation with NURBS(Non-Uniform Rational B-Spline) and various stretching functions to control grid points distribution on curves and surfaces. And some menus are added to perform flexible management, for the objects. We generated surface and field grid system about full aircraft configuration using KGRID. The performance and stability of the KGRID is verified through the generation of the grid system about a complex shape.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.25 no.2
• /
• pp.109-118
• /
• 2013

The absence of a simple and general analytical model has been a problem in the design and analysis of desiccant-assisted air-conditioning systems. In this study, such an analytical model has been developed based on the approximate integral solution of the coupled transient ordinary differential equations for the heat and mass transfer processes in a desiccant wheel. It turned out that the initial conditions should be determined by the solution of four linear algebraic equations including the heat and mass transfer equations for the air flow as well as the energy and mass conservation equations for the desiccant bed. It is also shown that time-averaged exit air temperature and humidity relations could be given in terms of the heat and mass transfer effectiveness.

• Coupled systems mechanics
• /
• v.8 no.2
• /
• pp.147-168
• /
• 2019

This work aims to simulate three-dimensional heavy oil flow in a reservoir with heater-wells. Mass, momentum and energy balances, as well as correlations for rock and fluid properties, are used to obtain non-linear partial differential equations for the fluid pressure and temperature, and for the rock temperature. Heat transfer is simulated using a two-equation model that is more appropriate when fluid and rock have very different thermal properties, and we also perform comparisons between one- and two-equation models. The governing equations are discretized using the Finite Volume Method. For the numerical solution, we apply a linearization and an operator splitting. As a consequence, three algebraic subsystems of linearized equations are solved using the Conjugate Gradient Method. The results obtained show the suitability of the numerical method and the technical feasibility of heating the reservoir with static equipment.

• Coupled systems mechanics
• /
• v.12 no.5
• /
• pp.419-428
• /
• 2023

We have developed a model for estimating the parameters of viscous materials from indirect tensile tests for asphalt. This is a simple Burger nonlinear rheological two-cell model or standard model. At the same time, we begin to develop a more versatile and complex multi-cell model. The simple model is validated using experimental load-displacement results from laboratory tests: The recorded displacements are used as input values and the measured force data are simulated with the model. The optimal model parameters are estimated using the Levenberg-Marquardt method and a very good agreement between the experimental results and the model calculations is shown. However, not all parts of the model are active in the loading phase of the experiment, so we extended the validation of the model to the simulation of the relaxation behaviour. In this stage, the other model parameters are activated and the simulation results are consistent with the literature. At this stage, we have estimated the parameters only for the two-cell uniaxial model, but further work will include results for the multi-cell model.