• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.19 no.2
• /
• pp.197-215
• /
• 2015

Modeling water flow in variably saturated, porous media is important in many branches of science and engineering. Highly nonlinear relationships between water content and hydraulic conductivity and soil-water pressure result in very steep wetting fronts causing numerical problems. These include poor efficiency when modeling water infiltration into very dry porous media, and numerical oscillation near a steep wetting front. A one-dimensional finite element formulation is developed for the numerical simulation of variably saturated flow systems. First order backward Euler implicit and second order Crank-Nicolson time discretization schemes are adopted as a solution strategy in this formulation based on Picard and Newton iterative techniques. Five examples are used to investigate the numerical performance of two approaches and the different factors are highlighted that can affect their convergence and efficiency. The first test case deals with sharp moisture front that infiltrates into the soil column. It shows the capability of providing a mass-conservative behavior. Saturated conditions are not developed in the second test case. Involving of dry initial condition and steep wetting front are the main numerical complexity of the third test example. Fourth test case is a rapid infiltration of water from the surface, followed by a period of redistribution of the water due to the dynamic boundary condition. The last one-dimensional test case involves flow into a layered soil with variable initial conditions. The numerical results indicate that the Crank-Nicolson scheme is inefficient compared to fully implicit backward Euler scheme for the layered soil problem but offers same accuracy for the other homogeneous soil cases.

• Journal of Korea Water Resources Association
• /
• v.32 no.6
• /
• pp.607-616
• /
• 1999

A fractional step finite difference model for the longitudinal dispersion of nonconservative contaminants is developed. It is based on splitting the longitudinal dispersion equation into a set of three equations each to be solved over a one-third time step. The fourth-order Holly-Preissmann scheme, an analytic solution, and the Crank-Nicholson scheme are used to solve the equations for the pure advection, the first-order decay, and the diffusion, respectively. To test the model, it is applied to simulate the longitudinal dispersion of continuous source released into a nonuniform flow field as well as the dispersion of an instantaneous source in a uniform flow field. The results are compared with the exact solution and those computed by an existing model. Compared to the existing model which uses Euler method for the first-order decay equation, the present model yield more accurate results as the decay coefficient increases.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.11 no.2
• /
• pp.116-124
• /
• 2008

Current paper investigates the dynamic behavior and stability of an infrared search and track subjected to external disturbance having gimbal structure with three rotating axes keeping constant angular velocity in the azimuth direction. Euler-Lagrange equation is applied to derive the coupled nonlinear dynamic equation of motion of infrared search and track and the characteristics of dynamic coupling are investigated. Two equilibrium points with small variations from the nonlinear coupling system are derived and the specific condition from which a coupled equation can be three independent equations is derived. Finally, to examine the stability of system, Lyapunov direct method was used and system stability and stability boundaries are investigated.

• Proceedings of the KIEE Conference
• /
• 2002.07b
• /
• pp.882-884
• /
• 2002

The transient voltage characteristics of capacitor self-exited induction generator are analyzed by the state equation which is obtained from the d-q axis equivalent circuit of stationary reference frame and torque equation. The d-q equivalent circuit is composed using the condition of stationary reference frame. The mutual inductance is only considered as a function of magnetizing current in the equivalent circuit. The characteristics are analyzed and discussed by the backward Euler method for various load conditions under specified initial conditions and input.

• Bulletin of the Korean Mathematical Society
• /
• v.47 no.6
• /
• pp.1225-1234
• /
• 2010

Numerical solutions for the fractional differential dispersion equations with nonlinear forcing terms are considered. The backward Euler finite difference scheme is applied in order to obtain numerical solutions for the equation. Existence and stability of the approximate solutions are carried out by using the right shifted Grunwald formula for the fractional derivative term in the spatial direction. Error estimate of order $O({\Delta}x+{\Delta}t)$ is obtained in the discrete $L_2$ norm. The method is applied to a linear fractional dispersion equations in order to see the theoretical order of convergence. Numerical results for a nonlinear problem show that the numerical solution approach the solution of classical diffusion equation as fractional order approaches 2.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 2000.11a
• /
• pp.413-416
• /
• 2000

The spent fuel skeleton is processed in the cutting processing after compacting. If the cutting length is processed in the same interval length. The spent fuel skeleton is stayed on the connection of bottom nozzle and guide tube. In the case, because the compressive stress is loaded along the length, the guide tube is generated the buckling stress and the deforming. But the deformed guide tube interrupted the guide tube inserted through compressive room. therefore, it is experimented for the optimum buckling stress and the preventing of guide deformed. This paper is predicted the all over buckling stress of the spent fuel skeleton by using experiment. The guide of Spent fuel skeleton have buckling characteristics of the medium column. The experiment and analysis is conducted by the comparing among the equation of Euler, Johnson and Engresser. The fittest one of method is Engresser equation.

• Proceedings of the KIPE Conference
• /
• 2003.07b
• /
• pp.659-662
• /
• 2003

This paper analyzes the electro-mechanical characteristics of the spindle motor in a computer hard disk drive due to the trapezoidal and sinusoidal driving methods. The driving circuit equation is modified by considering the switching action of PWM inverter, and is coupled with the Maxwell equation for the analysis of the magnetic field. Mechanical motion of a rotor is calculated by solving Newton-Euler equation. Electro-mechanical excitation and dynamic response are characterized by analyzing the free response of a rotating rotor and Fourier analysis of the excitation force.

• Proceedings of the KIEE Conference
• /
• 1997.07b
• /
• pp.474-477
• /
• 1997

In this paper, an efficient attitude determination algorithm based on the Kalman Filter which combines earth/sun sensor data with gyro data in a mutually compensating manner is presented. Quaternion is used as the attitude state to save computation time and to prevent the gimbal-lock situation associated with Euler angles. Gyro data allows the use of the kinematic equation instead of space vehicle's dynamic equation which is usually based on approximation of the actual dynamics and inaccurate torque information. The gyro data are used to propagate the attitude through kinematic equation and the earth/sun sensor data are used to update the attitude and estimate the gyro bias. Simulation results for the KOMPSAT attitude determination system are presented.

• The Transactions of the Korean Institute of Electrical Engineers B
• /
• v.49 no.11
• /
• pp.742-748
• /
• 2000

This paper presents a numerical method to analyze the electromechanical-coupled field in the spindle motor of a computer hard drive and investigates dynamic response due to the electromechanical excitation, i.e. unbalanced magnetic force and centrifugal force for the rotational asymmetric motor. Magnetic field is calculated from Maxwells equation and voltage equation by introducing nonlinear time-dependent finite element analysis. Mechanical motion of rotor is calculated by solving Newton-Euler equation. Electromechanical excitation and dynamic response are characterized by analyzing the free response of a rotating rotor and Fourier analysis of the excitation force and resulting vibration of a rotor. It shows that centrifugal force produces the unbalanced magnetic force even in the rotational symmetric motor. It also shows that resonance produces quite considerable vibration even when the high excitation frequency with small amplitude matches with the natural frequency of the spindle motor.

• Steel and Composite Structures
• /
• v.46 no.5
• /
• pp.637-647
• /
• 2023

Nonlinear free vibration and stability responses of a carbon nanotube reinforced composite beam under temperature rising are investigated in this paper. The material of the beam is considered as a polymeric matrix by reinforced the single-walled carbon nanotubes according to different distributions with temperature-dependent physical properties. With using the Hamilton's principle, the governing nonlinear partial differential equation is derived based on the Euler-Bernoulli beam theory. In the nonlinear kinematic assumption, the Von Kármán nonlinearity is used. The Galerkin's decomposition technique is utilized to discretize the governing nonlinear partial differential equation to nonlinear ordinary differential equation and then is solved by using of multiple time scale method. The critical buckling temperatures, the nonlinear natural frequencies and the nonlinear free response of the system is obtained. The effect of different patterns of reinforcement on the critical buckling temperature, nonlinear natural frequency, nonlinear free response and phase plane trajectory of the carbon nanotube reinforced composite beam investigated with temperature-dependent physical property.