• 제목/요약/키워드: forward differential equation

검색결과 26건 처리시간 0.022초

선형적 물성변화를 고려하는 유한요소법을 이용한 2.5차원 전자탐사 수치모델링 (2.5 Dimensional Electromagnetic Finite Element Numerical modeling using linear conductivity variation)

  • 고광범;서백수
    • 산업기술연구
    • /
    • 제18권
    • /
    • pp.131-138
    • /
    • 1998
  • Numerical modeling for electromagnetic exploration methods are essential to understand behaviours of electromagnetic fields in complex subsurfaces. In this study, a finite element method was adopted as a numerical scheme for the 2.5-dimensional forward problem. And a finite element equation considering linear conductivity variation was proposed when 2.5-dimensional differential equation to couple eletric and magnetic field was implemented. Model parameters were investigated for near-field with large source effects and far-field with responses dominantly by homogeneous half-space. Numerical responses by this study were compared with analytic solutions in homogeneous half-space and compared with other three dimensional numerical results.

  • PDF

A MASS LUMPING AND DISTRIBUTING FINITE ELEMENT ALGORITHM FOR MODELING FLOW IN VARIABLY SATURATED POROUS MEDIA

  • ISLAM, M.S.
    • Journal of the Korean Society for Industrial and Applied Mathematics
    • /
    • 제20권3호
    • /
    • pp.243-259
    • /
    • 2016
  • The Richards equation for water movement in unsaturated soil is highly nonlinear partial differential equations which are not solvable analytically unless unrealistic and oversimplifying assumptions are made regarding the attributes, dynamics, and properties of the physical systems. Therefore, conventionally, numerical solutions are the only feasible procedures to model flow in partially saturated porous media. The standard Finite element numerical technique is usually coupled with an Euler time discretizations scheme. Except for the fully explicit forward method, any other Euler time-marching algorithm generates nonlinear algebraic equations which should be solved using iterative procedures such as Newton and Picard iterations. In this study, lumped mass and distributed mass in the frame of Picard and Newton iterative techniques were evaluated to determine the most efficient method to solve the Richards equation with finite element model. The accuracy and computational efficiency of the scheme and of the Picard and Newton models are assessed for three test problems simulating one-dimensional flow processes in unsaturated porous media. Results demonstrated that, the conventional mass distributed finite element method suffers from numerical oscillations at the wetting front, especially for very dry initial conditions. Even though small mesh sizes are applied for all the test problems, it is shown that the traditional mass-distributed scheme can still generate an incorrect response due to the highly nonlinear properties of water flow in unsaturated soil and cause numerical oscillation. On the other hand, non oscillatory solutions are obtained and non-physics solutions for these problems are evaded by using the mass-lumped finite element method.

Numerical simulation of advection-diffusion on flow in waste stabilization ponds (1-dimension) with finite difference method forward time central space scheme

  • Putri, Gitta Agnes;Sunarsih, Sunarsih;Hariyanto, Susilo
    • Environmental Engineering Research
    • /
    • 제23권4호
    • /
    • pp.442-448
    • /
    • 2018
  • This paper presents the numerical simulation of advection-diffusion mechanism of BOD concentration which was used as an indicator of waste only in one flow-direction of waste stabilization ponds (1-dimension (1-D)). This model was represented in partial differential equation order 2. The purpose of this paper was to determine the simulation of the model 1-D of wastewater transport phenomena based advection-diffusion mechanism and did validate the model. Numerical methods which was used for the solution of this model is finite difference method with Forward Time Central Space scheme. The simulation results which was obtained would be compared with field observation data as a validation model. Collection of field data was carried out in the Wastewater Treatment Plant Sewon, Bantul, D.I. Yogyakarta. The results of numerical simulations were indicate that the advection-diffusion mechanism takes place continuously over time. Then validation of the model was state that there was a difference between the calculation results with the field data, with a correlation value of 0.998.

2.5차원 전자탐사 수치모델링 및 역해 (2.5-Dimensional Electromagnetic Numerical Modeling and Inversion)

  • 고광범;서정희;신창수
    • 지구물리와물리탐사
    • /
    • 제2권1호
    • /
    • pp.43-53
    • /
    • 1999
  • 전자탐사 수치모델링 및 역산은 복잡한 지하구조에 대한 전자기장의 정확한 거동양상 규명을 위한 기본적인 해석방법이다. 본 연구에서는 2.5차원 전자탐사 수치모델링 알고리즘을 유한요소법을 이용하여 구하였다. 이때 전기장과 자기장이 결합되어 표현되는 2.5차원 미분방정식을 정식화할 때 선형적 물성변화를 고려하는 유한요소 방정식을 도출, 제의하였다. 송신원의 영향이 큰 근거리장과 균질 반무한공간의 반응이 지배적인 원거리장 각각에 대하여 모델변수를 검토하고 이를 토대로 균질 반무한공간에서의 해석적인 해와 비교하여 봄으로써 알고리즘의 타당성을 검증하였다. 본 연구에서 제시한 수치 모형반응 알고리즘에 부합하는 수정된 블록역산모형을 고안하였으며 고립형 이상체에 대한 분해능을 살펴봄으로써 분산방정식으로 표현되는 전자탐사자료의 역해 가능성 및 한계를 고찰하였다.

  • PDF

비선형 천장 크레인시스템의 위치제어 알고리즘 개발 (Development of the Position Control Algorithm for Nonlinear Overhead Crane Systems)

  • 이종규;이상룡
    • 한국정밀공학회지
    • /
    • 제17권4호
    • /
    • pp.142-147
    • /
    • 2000
  • An overhead crane system which transports an object by girder motion, trolley motion, and hoist motion becomes a nonlinear system because the length of a rope changes. To develope the position control algorithm for the nonlinear crane systems, we apply a nonlinear optimal control method which uses forward and backward difference methods and obtain optimal inputs. This method is suitable for the overhead crane system which is characterized by the differential equation of higher degree and swing motion. From the results of computer simulation, it is founded that the position of the overhead crane system is controlled, and the swing of the object is suppressed.

  • PDF

A Diffusion Model for a System Subject to Random Shocks

  • Lee, Eui-Yong;Song, Mun-Sup;Park, Byung-Gu
    • Journal of the Korean Statistical Society
    • /
    • 제24권1호
    • /
    • pp.141-147
    • /
    • 1995
  • A diffusion model for a system subject to random shocks is introduced. It is assumed that the state of system is modeled by a Brownian motion with negative drift and an absorbing barrier at the origin. It is also assumed that the shocks coming to the system according to a Poisson process decrease the state of the system by a random amount. It is further assumed that a repairman arrives according to another Poisson process and repairs or replaces the system i the system, when he arrives, is in state zero. A forward differential equation is obtained for the distribution function of X(t), the state of the systme at time t, some boundary conditions are discussed, and several interesting characteristics are derived, such as the first passage time to state zero, F(0,t), the probability of the system being in state zero at time t, and F(0), the limit of F(0,t) as t tends to infinity.

  • PDF

규칙파중 시간영역 스트립이론을 이용한 비선형 유탄성 해석 (Nonlinear Hydroelastic Analysis Using a Time-domain Strip Theory m Regular Waves)

  • 조일형;한성곤;권승민
    • 한국해양공학회지
    • /
    • 제19권4호
    • /
    • pp.1-8
    • /
    • 2005
  • A nonlinear time-domain strip theory for vertical wave loads and ship responses is to be investigated. The hydrodynamic memory effect is approximated by a higher order differential equation without convolution. The ship is modeled as a non-uniform Timoshenko beam. Numerical calculations are presented for the S175 Containership translating with the forward speed in regular waves. The approach described in this paper can be used in evaluating ship motions and wave loads in extreme wave conditions and validating nonlinear phenomena in ship design.

Dynamic Elastica에 의한 유연매체의 거동해석 (Analysis of Flexible Media Behavior by Dynamic Elastica)

  • 홍성권;지중근;장용훈;박노철;박영필
    • 한국소음진동공학회:학술대회논문집
    • /
    • 한국소음진동공학회 2004년도 추계학술대회논문집
    • /
    • pp.600-605
    • /
    • 2004
  • In many machines handling lightweight and flexible media such as magnetic tape drives, xerographic copiers and sewing machines, the media must transit an open space. It is important to predict the static and dynamic behavior of the sheets with a high degree of reliability. The nonlinear theory of the dynamic elastica has often been used to a nonlinear dynamic deflection model. In this paper, the governing equation is derived and simulated by the finite differential method. The parametric cubic curve is applied for defining the guide shape. The dynamic contact conditions suggested by Klarbring is used to predict the direction of the flexible media according to the initial velocity and the friction coefficient. The analysis is also compared to the conventional model, showing that after contacting a $45^{\circ}$ wall, the directions of flexible media of two models are different.

  • PDF

A Fast Scheme for Inverting Single-Hole Electromagnetic Data

  • Kim Hee Joon;Lee Jung-Mo;Lee Ki Ha
    • 대한자원환경지질학회:학술대회논문집
    • /
    • 대한자원환경지질학회 2002년도 춘계 공동학술발표회
    • /
    • pp.167-169
    • /
    • 2002
  • The extended Born, or localized nonlinear approximation of integral equation (IE) solution has been applied to inverting single-hole electromagnetic (EM) data using a cylindrically symmetric model. The extended Born approximation is less accurate than a full solution but much superior to the simple Born approximation. When applied to the cylindrically symmetric model with a vertical magnetic dipole source, however, the accuracy of the extended Born approximation is greatly improved because the electric field is scalar and continuous everywhere. One of the most important steps in the inversion is the selection of a proper regularization parameter for stability. Occam's inversion (Constable et al., 1987) is an excellent method for obtaining a stable inverse solution. It is extremely slow when combined with a differential equation method because many forward simulations are needed but suitable for the extended Born solution because the Green's functions, the most time consuming part in IE methods, are repeatedly re-usable throughout the inversion. In addition, the If formulation also readily contains a sensitivity matrix, which can be revised at each iteration at little expense. The inversion algorithm developed in this study is quite stable and fast even if the optimum regularization parameter Is sought at each iteration step. Tn this paper we show inversion results using synthetic data obtained from a finite-element method and field data as well.

  • PDF

Nonlocal elasticity approach for free longitudinal vibration of circular truncated nanocones and method of determining the range of nonlocal small scale

  • Li, C.;Sui, S.H.;Chen, L.;Yao, L.Q.
    • Smart Structures and Systems
    • /
    • 제21권3호
    • /
    • pp.279-286
    • /
    • 2018
  • The free longitudinal vibration of a circular truncated nanocone is investigated based on the nonlocal elasticity theory. Exact analytical formulations for tapered nanostructures are derived and the nonlinear differential governing equation of motion is developed. The nonlocal small scale effect unavailable in classical continuum theory is addressed to reveal the long-range interaction of atoms implicated in nonlocal constitutive relation. Unlike most previous studies applying the truncation method to the infinite higher-order differential equation, this paper aims to consider all higher-order terms to show the overall nonlocality. The explicit solution of nonlocal stress for longitudinal deformation is determined and it is an infinite series incorporating the classical stress derived in classical mechanics of materials and the infinite higher-order derivative of longitudinal displacement. Subsequently, the first three modes natural frequencies are calculated numerically and the significant effects of nonlocal small scale and vertex angle on natural frequencies are examined. The coupling phenomenon of natural frequency is observed and it is induced by the combined effects of nonlocal small scale and vertex angle. The critical value of nonlocal small scale is defined, and after that a new proposal for determining the range of nonlocal small scale is put forward since the principle of choosing the nonlocal small scale is still unclear at present. Additionally, two different types of nonlocal effects, namely the nonlocal stiffness weakening and strengthening, reversed phenomena existing in nanostructures are observed and verified. Hence the opposite nonlocal effects are resolved again clearly. The nano-engineers dealing with a circular truncated nanocone-based sensors and oscillators may benefit from the present work.