• Journal of the Korean Institute of Telematics and Electronics A
• /
• v.32A no.3
• /
• pp.143-151
• /
• 1995

The electrochemical impulse oscillations of the cathodic currents at the platinum group (Pt, Pd) electrode/(0.05M KHC$_{8}H_{4}O_{4}$) buffer solution interfaces have been studied using voltammetric, chronoamperometric, and electrochemical impedance methods. The periodic impulses of the cathodic currents are the activation controlled currents due to the hydrogen evolution reaction, and depend on the fractional surface coverage of the adsorbed hydrogen intermediate and potential. The oscillatory mechanism of the cathodic current impulses is connected with the unstable steady state of negative differential resistance. The widths and periods of the cathodic current impulses are 4ms or 5ms and 152.5ms or 305ms, respectively. The H$^{+}$ discharge reaction step is 38 or 61 times faster thatn the recombination reaction steps and the H$^{+}$ mass transport processes. The atom-atom recombination reaction step is twice faster thatn the atom-ion recombination reaction step. The two kinds of active sites corresponding to the atom-atom and atom-ion recombination reaction steps exist on the platinum group electrode surfaces.

• Journal of the Korean Mathematical Society
• /
• v.55 no.2
• /
• pp.471-506
• /
• 2018

In this article, some projection methods (or fractional-step methods) are proposed and analyzed for the micropolar Navier-Stokes equations (MNSE). These methods allow us to decouple the MNSE system into two sub-problems at each timestep, one is the linear and angular velocities system, the other is the pressure system. Both first-order and second-order projection methods are considered. For the classical first-order projection scheme, the stability and error estimates for the linear and angular velocities and the pressure are established rigorously. In addition, a modified first-order projection scheme which leads to some improved error estimates is also proposed and analyzed. We also present the second-order projection method which is unconditionally stable. Ample numerical experiments are performed to confirm the theoretical predictions and demonstrate the efficiency of the methods.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.16 no.3
• /
• pp.193-204
• /
• 2012

This paper is comparing numerical schemes for a differential equation with convection and fourth-order diffusion. Our model equation is $h_t+(h^2-h^3)_x=-(h^3h_{xxx})_x$, which arises in the context of thin film flow driven the competing effects of an induced surface tension gradient and gravity. These films arise in thin coating flows and are of great technical and scientific interest. Here we focus on the several numerical methods to apply the model equation and the comparison and analysis of the numerical results. The convection terms are treated with well known WENO methods and the diffusion term is treated implicitly. The diffusion and convection schemes are combined using a fractional step-splitting method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.3
• /
• pp.173-182
• /
• 2014

The mixed convolved action principle provides a new rigorous weak variational formalism for a broad range of initial boundary value problems in mathematical physics and mechanics in terms of mixed formulation, convolution, and fractional calculus. In this paper, its potential in the development of numerical methods for transient problems in various dynamical systems when adopting temporally second order approximation is investigated. For this, the classical single-degree-of-freedom linear elastic dynamical systems are primarily considered to investigate computational characteristics of the developed algorithms. For the undamped system, all the developed algorithms are symplectic with respect to the time step. For the damped system, they are shown to be accurate with good convergence characteristics.

• Wind and Structures
• /
• v.15 no.3
• /
• pp.189-208
• /
• 2012

Lock-in and drag amplification phenomena are studied for a flexible cantilever using a simplified fluid-structure interaction approach. Instead of solving the 3D domain, a simplified setup is devised, in which 2D flow problems are solved on a number of planes parallel to the wind direction and transversal to the structure. On such planes, the incompressible Navier-Stokes equations are solved to estimate the fluid action at different positions of the line-like structure. The fluid flow on each plane is coupled with the structural deformation at the corresponding position, affecting the dynamic behaviour of the system. An Arbitrary Lagrangian-Eulerian (ALE) approach is used to take in account the deformation of the domain, and a fractional-step scheme is used to solve the fluid field. The stabilization of incompressibility and convection is achieved through orthogonal quasi-static subscales, an approach that is believed to provide a first step towards turbulence modelling. In order to model the structural problem, a special one-dimensional element for thin walled cross-section beam is implemented. The standard second-order Bossak method is used for the time integration of the structural problem.

• Journal of computational fluids engineering
• /
• v.10 no.1
• /
• pp.45-50
• /
• 2005

A fundamental understanding of the flow around the wind turbine is important to investigate the performance of new type of wind turbine. This study presents the simulation of three dimensional flow fields around the Darius wind turbine as an example. Incompressible Navier-Stokes equations are used for this simulation. The rotating coordinate system that rotates in the same speed of the turbine is used in order to simplify the boundary condition on the blades. Additionally, the boundary fitted coordinate system is employed in order to express the shape of the blades precisely. Fractional step method is used to solve the basic equations. Third order upwind scheme is chosen for the approximation of the non-linear terms since it can compute the flow field stably even at high Reynolds number without any turbulence models. The flow fields obtained in this study are highly complex due to the three dimensionality and are visualized effectively by using the technique of the computer graphics.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.228-229
• /
• 2003

Complex flow field around the Darius turbine rotating stationally are simulated by solving the three dimensional incompressible Navier-Stokes equation numerically. The rotating coordinate system is employed so that the boundary conditions on the blades of the rotor become simple. In order to impose the boundary condition on the blades precisely, the boundary fitted coordinate system is employed. Fractional step method is used to solve the basic equations. The complex flow fields due to the three dimensionality of the geometry of the turbine and the rotation of the turbine are obtained and they are visualized effectively by using the technique of the computer graphics.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.910-915
• /
• 2003

A Navier-Stokes equation solver for incompressible viscous flows with free surface is developed and tested. This is based upon a fractional time step method and a non-staggered finite volume formulation for unstructured meshes. For time advancement scheme, Adams -Bashforth method for convective term and Crank-Nicolson method for diffusive term are applied. The interface between two fluids with different fluid properties is tracked with Piecewise Linear Interface Calculation(PLIC) Volume-of-Fluid(VOF) methods. Computational results are presented for some test problems: the broken dam, the sloshing in a rectangular tank, the filling of a cylindrical tank.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.6
• /
• pp.753-765
• /
• 2003

A finite element code for the numerical solution of the Navier-Stokes equation is parallelized by vertex-oriented domain decomposition. To accelerate the convergence of iterative solvers like conjugate gradient method, parallel block ILU, iterative block ILU, and distributed ILU methods are tested as parallel preconditioners. The effectiveness of the algorithms has been investigated when P1P1 finite element discretization is used for the parallel solution of the Navier-Stokes equation. Two-dimensional and three-dimensional Laplace equations are calculated to estimate the speedup of the preconditioners. Calculation domain is partitioned by one- and multi-dimensional partitioning methods in structured grid and by METIS library in unstructured grid. For the domain-decomposed parallel computation of the Navier-Stokes equation, we have solved three-dimensional lid-driven cavity and natural convection problems in a cube as benchmark problems using a parallelized fractional 4-step finite element method. The speedup for each parallel preconditioning method is to be compared using upto 64 processors.

• International Journal of Reliability and Applications
• /
• v.7 no.2
• /
• pp.167-175
• /
• 2006

Effective and quick methods that are easy to carry out even by hand, or easy to be programmed by hand-held calculators are needed for assessing the sizes of contrasts of unreplicated $2^P$ factorial designs. Moreover, they have the advantage to use the original numerical measurements which makes the analysis easier to explain. Basically, Lenth (1989) is one of the most familiar of such quick and powerful methods. Later on, Aboukalam (2001) proposes under constant effects an alternative sophisticated method to Lenth's method. The proposed method is the supreme from two considerable powers. The first utmost indicates less inflation deficiency while the other utmost indicates less shrinkage deficiency. Also under constant effects, Al-Shiha (2006) introduces an alternative quick method which is less shrinkage deficiency while the inflation deficiency is the same. If effects are random, Aboukalam (2005) introduces an alternative quick method in which the first power is favored as long as the second power is within a small margin. In the spirit of quickness and fixed effects, this article adds another method which is supreme from the two considerable powers. The method is based on a one step of the scale-part of a suggested M-estimate for location. Explicitly, we suggest adapting the skipped median (ASKM) estimate. Critical values of ASKM-method, for several sample sizes often used, are empirically computed.