• Journal of Industrial Technology
• /
• v.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.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.20 no.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.

• Environmental Engineering Research
• /
• v.23 no.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.

• Geophysics and Geophysical Exploration
• /
• v.2 no.1
• /
• pp.43-53
• /
• 1999

Numerical modeling and inversion for electromagnetic exploration methods are essential to understand behaviour of electromagnetic fields in complex subsurface. 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. Blocky inversion model was modified to be applied to the forward calculation in this study and it was also adopted in the inversion algorithm. Resolution for isolated bodies were investigated to confirm possibility and limitation of inversion for electromagnetic exploration data.

• Journal of the Korean Society for Precision Engineering
• /
• v.17 no.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.

• Journal of the Korean Statistical Society
• /
• v.24 no.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.

• Journal of Ocean Engineering and Technology
• /
• v.19 no.4 s.65
• /
• 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.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2004.11a
• /
• 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.

• Proceedings of the KSEEG Conference
• /
• 2002.04a
• /
• 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.

• Smart Structures and Systems
• /
• v.21 no.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.