• Tribology and Lubricants
• /
• v.15 no.4
• /
• pp.343-353
• /
• 1999

The role of viscosity index improver's(Ⅶ) additives for modem engine lubrication is complex. Under the condition of atmosphere or low shear rate, the characteristics of Ⅶ added lubricant is verified and quoted frequently for mathematical model of lubricant behavior. However, recent research shows that added lubricant has the characteristics of shear thinning at high shear rate condition although it performs well enough over the whole range of working temperature. At high shear rate, they show significant decrease of apparent viscosity irrespective of temperature. Many experimental researches verify that Ⅶ added lubricant shows boundary film layer formation on the solid surface as well as shear thinning effect by its polymeric molecular characteristics. The intend of our research is to verify the effects of Ⅶ from the viewpoint of continuum mechanics, because conventional Reynolds'equation with only pressure-viscosity relation cannot fully predict the lubricant behavior under the Ⅶ added condition. In these aspects, Reynolds'equation of Newtonian fluid model lacks the reflection of real fluid behavior and there is no way to explain the non-linear characteristics of Ⅶ added lubricant. In this research, we mathematically modeled the Ⅶ added lubricant behaviors which are the characteristics of non-Newtonian fluid behavior at high shear rate and boundary film formation on the solid surface. The consideration of elastic deformation in the contact region is also included in our computation and finally the converged film pressure and the film thickness with elastic deformation are obtained. The results are compared with those of Newtonian fluid model.

• Journal of Mechanical Science and Technology
• /
• v.18 no.6
• /
• pp.1026-1035
• /
• 2004

An immersed boundary method for solving the Navier-Stokes and thermal energy equations is developed to compute the heat transfer over or inside the complex geometries in the Cartesian or cylindrical coordinates by introducing the momentum forcing, mass source/sink, and heat source/sink. The present method is based on the finite volume approach on a staggered mesh together with a fractional step method. The method of applying the momentum forcing and mass source/sink to satisfy the no-slip condition on the body surface is explained in detail in Kim, Kim and Choi (2001, Journal of Computational Physics). In this paper, the heat source/sink is introduced on the body surface or inside the body to satisfy the iso-thermal or iso-heat-flux condition on the immersed boundary. The present method is applied to three different problems : forced convection around a circular cylinder, mixed convection around a pair of circular cylinders, and forced convection around a main cylinder with a secondary small cylinder. The results show good agreements with those obtained by previous experiments and numerical simulations, verifying the accuracy of the present method.

• Journal of Institute of Control, Robotics and Systems
• /
• v.6 no.1
• /
• pp.8-14
• /
• 2000

In this paper, three adaptive sliding mode control algorithms, which self-tune both the sliding mode gain and the boundary layer thickness, are proposed. The first algorithm uses a gain adaptation rule is combined with the boundary layer thickness adaptatioin rule to satisfy the sliding condition. In the third algorithm, the computation burden of the second algorithm is reduced further, and therefore no extra cost is required for real-time implementation. Due to the mixed sliding mode gain and the boundary layer thickness adaptation scheme, the tracking error and the chattering of the control input can be reduced greatly.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.4
• /
• pp.775-786
• /
• 2002

A Mixed Volume and Boundary Integral Equation Method is applied for the effective analysis of elastic wave scattering problems and plane elastostatic problems in unbounded solids containing general anisotropic inclusions and voids or isotropic inclusions. It should be noted that this newly developed numerical method does not require the Green's function for anisotropic inclusions to solve this class of problems since only Green's function for the unbounded isotropic matrix is involved in their formulation for the analysis. This new method can also be applied to general two-dimensional elastodynamic and elastostatic problems with arbitrary shapes and number of anisotropic inclusions and voids or isotropic inclusions. In the formulation of this method, the continuity condition at each interface is automatically satisfied, and in contrast to finite element methods, where the full domain needs to be discretized, this method requires discretization of the inclusions only. Finally, this method takes full advantage of the pre- and post-processing capabilities developed in FEM and BIEM. Through the analysis of plane elastostatic problems in unbounded isotropic matrix with orthotropic inclusions and voids or isotropic inclusions, and the analysis of plane wave scattering problems in unbounded isotropic matrix with isotropic inclusions and voids, it will be established that this new method is very accurate and effective for solving plane wave scattering problems and plane elastic problems in unbounded solids containing general anisotropic inclusions and voids/cracks or isotropic inclusions.

• Computational Structural Engineering
• /
• v.4 no.2
• /
• pp.87-94
• /
• 1991

In this paper, a mixed variational principle applicable to the linear elasticity of inhomogeneous anisotropic materials is presented. For derivation of the general variational principle, a systematic procedure for the variational formulation of linear coupled boundary value problems developed by Sandhu et al. is employed. Consistency condition of the field operators with the boundary operators results in explicit inclusion of boundary conditions in the governing functional. Extensions of admissible state function spaces and specialization to a certain relation in the general governing functional lead to the desired mixed variational principle. In the physical sense, the present variational principle is analogous to the Reissner's recent formulation obtained by applying Lagrange multiplier technique followed by partial Legendre transform to the classical minimum potential energy principle. However, the present one is more advantageous for the application to the general anisotropic materials since Reissner's principle contains an implicit function which is not easily converted to an explicit form.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.7 no.4 s.19
• /
• pp.107-113
• /
• 2004

In order to model a connection part between a control fin and actuator, the related characteristics of a dynamic stiffness were extracted from experiments. These characteristics include the static stiffness of a control fin and the dynamic stiffness of an actuator, so they are called the mixed dynamic stiffness here. This mixed dynamic stiffness is used as the boundary condition of a control fin connected to an actuator when the flutter characteristics are analyzed. The simulated stiffness of an actuator is corrected from the experiment results and the mixed dynamic is finally formulated in the domain of frequencies.

• Journal of Mechanical Science and Technology
• /
• v.18 no.6
• /
• pp.990-1001
• /
• 2004

The retention of a molten pool vessel cooled by internal vessel reflooding and/or external vessel reactor cavity flooding has been considered as one of severe accident management strategies. The present numerical study investigates the effect of both internal and external vessel mixed cooling on an internally heated molten pool. The molten pool is confined in a hemispherical vessel with reference to the thermal behavior of the vessel wall. In this study, our numerical model used a scaled-down reactor vessel of a KSNP (Korea Standard Nuclear Power) reactor design of 1000 MWe (a Pressurized Water Reactor with a large and dry containment). Well-known temperature-dependent boiling heat transfer curves are applied to the internal and external vessel cooling boundaries. Radiative heat transfer has been considered in the case of dry internal vessel boundary condition. Computational results show that the external cooling vessel boundary conditions have better effectiveness than internal vessel cooling in the retention of the melt pool vessel failure.

• Journal of applied mathematics & informatics
• /
• v.34 no.5_6
• /
• pp.383-395
• /
• 2016

The convergence rate of a numerical procedure based on Schwarz Alternating Method (SAM) for solving elliptic boundary value problems depends on the selection of the interface conditions applied on the interior boundaries of the overlapping subdomains. It has been observed that the Robin condition (mixed interface condition), controlled by a parameter, can optimize SAM's convergence rate. In [8], one formulated the twolayer multi-parameterized SAM and determined the optimal values of the multi-parameters to produce the best convergence rate for one-dimensional elliptic boundary value problems. Two-dimensional implementation was presented in [10]. In this paper, we present an implementation for threedimensional problem.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.1
• /
• pp.89-98
• /
• 2004

In this paper, we develope a dispersion model based on the Generalized Langevin Model. Thomson's well-mixed condition is the well known criterion to determine particle dispersion. But, it has 'non-uniqueness problem'. To resolve this, we adopt a turbulent model which is a new approach in this field of study. Our model was greatly simplified under the self-similarity condition, leaving model only two model constants $C_{0}$ and ${\gamma}$$_{5}$ that control the dispersion and spin which measures rotational property of the Lagrangian particle trajectory. We investigated the sign of spin as well as magnitude by using the Direct Numerical Simulation. Model calculations were performed on the neutrally stable boundary layer flow. We found that spin has weak effect on the particle dispersion but it shows the significant effect on the horizontal flux compared to the zero-spin model.

• Journal of applied mathematics & informatics
• /
• v.33 no.1_2
• /
• pp.33-44
• /
• 2015

The convergence rate of a numerical procedure based on Schwarz Alternating Method(SAM) for solving elliptic boundary value problems depends on the selection of the interface conditions applied on the interior boundaries of the overlapping subdomains. It has been observed that the Robin condition (mixed interface condition), controlled by a parameter, can optimize SAM's convergence rate. In [7], one formulated the multi-parameterized SAM and determined the optimal values of the multi-parameters to produce the best convergence rate for one-dimensional elliptic boundary value problems. Two-dimensional implementation was presented in [8]. In this paper, we present an implementation for three-dimensional problem.