• Title/Summary/Keyword: 대수 미분 방정식

Search Result 50, Processing Time 0.022 seconds

Numerical Simulations for Design of a Liquid Rocket Engine (액체 로켓엔진 설계에서의 유동해석)

  • 김영목;채연석
    • Proceedings of the Korean Society of Propulsion Engineers Conference
    • /
    • 1995.11a
    • /
    • pp.87-96
    • /
    • 1995
  • 인공위성의 궤도진입에 사용되는 액체추진제 로켓엔진의 개발에서 분사기 설계를 적절히 수정, 보완 할 목적으로 수행된 핵심부품별 유동해석의 내용이 기술되었다. 단일 격자계를 구성하기 어려운 복잡한 형상의 분사기 유동장에 대한 격자계 구성을 용이하게 하고, 3차원의 점성 유동해석을 컴퓨터 기억 용량에 제한없이 수행하기 위한 다중블럭 격자기법이 사용되었다. 분사기의 내부유동은 3차원 비압축성 Navier-Stokes 방정식으로 pseudocompressibility 방법을 이용하여 수치모사되었다. 정상상태의 해는 근사 인자분해에 의한 ADI 기법으로 계산되고, 공간미분항에 대해 nonstaggered 격자계에서 2차 중앙차분을 사용하며 수치해의 안정성을 위해 인공점성항을 추가하였다. 난류계산을 위해 Baldwin- Lomax의 대수적 난류모델에 다수의 벽면효과를 고려하였다. 해석결과는 분사기의 성능에 영향을 미칠 수 있는 유동조건에 따라 분석되었다.

  • PDF

Nonlinear Vortical Forced Oscillation of Floating Bodies (부유체의 대진폭 운동에 기인한 동유체력)

  • 이호영;황종흘
    • Journal of the Society of Naval Architects of Korea
    • /
    • v.30 no.2
    • /
    • pp.86-97
    • /
    • 1993
  • A numerical method is developed for the nonlinear motion of two-dimensional wedges and axisymmetric-forced-heaving motion using Semi-Largrangian scheme under assumption of potential flows. In two-dimensional-problem Cauchy's integral theorem is applied to calculate the complex potential and its time derivative along boundary. In three-dimensional-problem Rankine ring sources are used in a Green's theorem boundary integral formulation to salve the field equation. The solution is stepped forward numerically in time by integrating the exact kinematic and dynamic free-surface boundary condition. Numerical computations are made for the entry of a wedge with a constant velocity and for the forced harmonic heaving motion from rest. The problem of the entry of wedge compared with the calculated results of Champan[4] and Kim[11]. By Fourier transform of forces in time domain, added mass coefficient, damping coefficient, second harmonic forces are obtained and compared with Yamashita's experiment[5].

  • PDF

Dynamic Analysis of Constrained Mechanical System Moving on a Flexible Beam Structure(II) : Application (유연한 보 구조물 위를 이동하는 구속 기계계의 동력학 해석(II) : 응용)

  • Park, Chan-Jong;Park, Tae-Won
    • Journal of the Korean Society for Precision Engineering
    • /
    • v.17 no.11
    • /
    • pp.176-184
    • /
    • 2000
  • Recently, it becomes a very important issue to consider the mechanical systems such as high-speed vehicle and railway train moving on a flexible beam structure. Using general approach proposed in the first part of this paper, it tis possible to predict planar motion of constrained mechanical system and elastic structure with various kinds of foundation supporting condition. Combined differential-algebraic equations of motion derived from both multibody dynamics theory and Finite Element Method can be analyzed numerically using generalized coordinate partitioning algorithm. To verify the validity of this approach, results from simply supported elastic beam subjected to a moving load are compared with exact solution from a reference. Finally, parameter study is conducted for a moving vehicle model on a simply supported 3-span bridge.

  • PDF

Dynamic Analysis of Constrained Mechanical System Moving on a Flexible Beam Structure(I) : General Approach (유연한 보 구조물 위를 이동하는 구속 기계계의 동력학 해석(I) : 일반적인 접근법)

  • Park, Chan-Jong;Park, Tae-Won
    • Journal of the Korean Society for Precision Engineering
    • /
    • v.17 no.11
    • /
    • pp.165-175
    • /
    • 2000
  • In recent years, it becomes a very important issue to consider the mechanical systems such as high-speed vehicles and railway trains moving on elastic beam structures. In this paper, a general approach, which can predict the dynamic behavior of constrained mechanical system and elastic beam structure, is proposed. Also, various supporting conditions of a foundation support are considered for the elastic beam structures. The elastic structure is assumed to be a nonuniform and linear Bernoulli-Euler beam with proportional damping effect. Combined Differential-Algebraic Equations of motion are derived using multibody dynamics theory and Finite Element Method. The proposed equations of motion can be solved numerically using generalizd coordinate partitioning method and Predictor-Corrector algorithm, which is an implicit multi-step integration method.

  • PDF

The Influence of Instrumentalization of Computer Algebra System(CAS) on the Sequence of Mathematics Curriculum in the Optimization Problem Solving Activities of CAS (최적화 문제해결 활동에서 "CAS의 도구화"가 교육과정 내용제시 순서에 미치는 영향)

  • Han, Se-Ho
    • Journal of Educational Research in Mathematics
    • /
    • v.20 no.2
    • /
    • pp.185-202
    • /
    • 2010
  • This study was designed to investigate the possibility that the optimization problem solving activities based on the instrumented CAS can have an influence on the sequence of mathematics curriculum in secondary mathematics education. Some optimization problem solving activities based on CAS were constructed and executed to eleventh grade(the penultimate year of Korean high school) 7 students for nine class hours. They have experienced using CAS in mathematics class for three months, but never learned calculus. The data which consists of classroom observations(audio and video taped) and post-unit interviews with students were analyzed. In the analysis, with CAS, students can highly deal with the applied optimization problems made up of calculus, cubic equation, solution of radical equation, and graph analysis which never learned. This result shows CAS may have an influence on the sequence of mathematics curriculum in secondary mathematics education.

  • PDF

Elastic Behavior of Contact Lense(I) : Effect of Vibration (콘택트 렌즈의 탄성에 관한 연구(I) : 진동에 의한 영향)

  • Kim, Dase-Soo
    • Journal of Korean Ophthalmic Optics Society
    • /
    • v.4 no.2
    • /
    • pp.1-16
    • /
    • 1999
  • Differential equations and its numerical solution program using Turbo-C were formulated to describe the radical distribution and average displacement amplitude of vibrating dehydrated contact lens(HEMA) driven by sinusoidal or rectangular pressure. The natural resonant frequency of the lens diaphram(thickness 0.08mm, diameter 14mm, curvature radius 8mm) was measured to be 5730 Hz from the extrapolation of frequency vs addedmass to the diaphram curve. The Young's modulus of the lens was measured to be $4{\times}10^9$ Pa with altering the original shape. The effect of parameters such as thickness, effective radius, damping coeff., amplitude of driving pressure on the vibration characteristics was illustrated by the computer simulation of the derived program. When the frequency of driving pressure coincides with the integral multiple of fundamental resonance frequency of the lens the wave pattern changes from arc to bell-shape along the radial position of the diaphram. If this happens to the contact lens on the cornea in vivo, it might create the feel of pull of the lens due to the increased rise of central part of the lens.

  • PDF

Analysis for A Partial Distribution Loaded Orthotropic Rectangular Plate with Various Boundary Condition (다양한 경계조건에서 부분 분포 하중을 받는 이방성 사각평판 해석)

  • See, Sangkwang
    • Journal of the Korea institute for structural maintenance and inspection
    • /
    • v.22 no.5
    • /
    • pp.13-22
    • /
    • 2018
  • In this study, a governing differential equation for the bending problem of orthotropic rectangular plate is drived. It's exact solution for various boundary conditions is presented. This solution follows traditional method like Navier's solution or Levy's solution that transforms the governing differential equation into an algebraic equation by using trigonometric series. To obtain a solution by Levy's method, it is required that two opposite edges of the plate be simply supported. And the boundary conditions, for which the Navier's method is applicable, are simply supported edge at all edges. In this study, it overcomes the limitations of the previous Navier's and Levy's methods.This solution is applicable for any combination of boundary conditions with simply supported edge and clamped edge in x, y direction. The plate could be subjected to uniform, partially uniform, and line loads. The advantage of the solution is that it is the exact solution as well as it overcomes the limitations of the previous Navier's and Levy's methods. Calculations are presented for orthotropic plates with nonsymmetric boundary conditions. Comparisons between the result of this paper and the result of Navier, Levy and Szilard solutions are made for the isotropic plates. The deflections were in excellent agreement.

Modeling of a Pervaporation Process for Concentrating Hydrogen Peroxide (과산화수소 농축을 위한 투과증발공정 모델링)

  • Nguyen, Huu Hieu;Lee, Sung Taek;Choi, Soo Hyoung
    • Korean Chemical Engineering Research
    • /
    • v.49 no.5
    • /
    • pp.560-564
    • /
    • 2011
  • The objective of this study is to propose a mathematical model for a pervaporation process for concentrating hydrogen peroxide. The process was developed by NASA, which consists of a shell and membrane tubes, where a liquid hydrogen peroxide solution flows in the shell, and a sweep gas flows in the tubes countercurrent to each other. The liquid retentate is concentrated as more water molecules permeate and evaporate through the membrane than hydrogen peroxide. For this process, a mathematical model has been developed in the form of a system of nonlinear partial differential algebraic equations based on a sorption-diffusion mechanism for permeation, an Arrhenius relationship for the temperature dependency of the permeate flux, and mass and momentum balances for the liquid concentrations and flows in the membrane module. The dynamic behavior of the concentration of hydrogen peroxide in the retentate side has been simulated by solving a simplified version of the proposed model, and the result is compared with the experimental data reported in the NASA patent.

Damped Oscill ations of the (Hard)Contact Lenses Posterior to the Blink (순목 후 콘택트(하드)렌즈의 감쇄 진동)

  • Kim, Dae-Soo
    • Journal of Korean Ophthalmic Optics Society
    • /
    • v.10 no.3
    • /
    • pp.173-184
    • /
    • 2005
  • A capillary action-induced tension develops in the tear layer between the contact lens and cornea, which leads to the restoring force due to difference in the layer thickness between either upper and lower or left and right side of the lens when it is displaced off the equilibrium position as a result of blinking. Suppose the lens was displaced a certain distance from the equilibrium position, lens starts to oscillate toward the equilibrium position with the decreasing amplitude due to the restoring force as well as the velocity dependent viscous damping force in the tear layer. A mathematical model which consists of the differential equations and their numerical solution was proposed to analyze the damped oscillations of lenses. The model predicts the time dependence of lenses after the blink varying the various parameters such as Be, diameters, masses and positions displaced from equilibrium. As the Be and mass of lens increases the rate of amplitude reduction decreases, which requires a more time for the lens to return to the equilibrium position. It seems that varying the lens' displacement and diameters affect the lens' motion very little.

  • PDF